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Abstract 

This  report  documents  the  W.H.O.I.  VAX-11  programs  used  to 
calculate  available  potential  energy  and  related  quantities  from  CTD  data 
using  the  technique  described  in  Bray  and  Fofonoff  (1981).  The  report 
includes  examples  of  how  the  programs  may  be  used,  as  well  as  complete 
listings  of  all  the  required  FORTRAN  files. 
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Introduction 

This  report  describes  the  structure  and  usage  of  programs  designed 
for  calculating  and  displaying  available  potential  energy  (APE), 
adiabatical ly  leveled  steric  surfaces,  and  related  variables  from  a  group 
of  CTD  stations.  For  a  general  discussion  of  the  technique  it  is 
strongly  recommended  that  the  reader  refer  to  Bray  and  Fofonoff  (1981). 
The  programs  have  an  inherent  requirement  that  the  input  CTD  data  be  an 
even  series  in  pressure,  although  the  input  pressure  interval  may  be 
specified.  This  report  describes  specifically  the  structure  of  the 
programs  as  used  on  the  W.H.O.I.  VAX-11,  with  input  data  in  the  standard 
CTD78  disc  format  (Millard,  et  al  (1978)).  Other  input  formats  can  be 
accomodated  through  modification  of  the  data  input  subroutine  as 
described  in  section  4. 

The  calculation  and  display  are  divided  into  separate  programs. 
POTential  ENergy  (POTEN)  reads  the  input  data,  calculates  the 
adiabatical ly  leveled  reference  steric  field  (see  Bray  and  Fofonoff, 

1981)  and  variables  related  to  the  leveled  field.  Potential  Energy  PLoT 
(PEPLT)  calculates  variables  derived  from  the  leveled  field  variables  and 
displays  POTEN  output  in  the  form  of  lists  and  plots. 

This  report  is  divided  into  four  sections.  The  first.  General 
Structure,  covers  the  non-FORTRAN  aspects  of  the  programs:  file 
structure,  linkage  and  general  usage.  The  second  and  third  sections 
contain  detailed  documentation  for  POTEN  and  PEPLT.  The  fourth  section 
describes  modifications  to  the  data  read  subroutine  in  POTEN,  to  allow 
input  data  in  other  than  CT078  disc  format.  Documented  examples  of  how 
to  run  the  programs  interactively  and  in  batch  mode  on  the  VAX-11  are 
found  in  Appendix  A.  Listings  of  programs  appear  in  Appendices  B  and  C. 
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1.  General  Structure  of  Programs 

Both  POTEN  and  PEPLT  are  accessed  through  a  short  main  program 
which  performs  initializations  of  parameters  as  requested  by  the  user. 
Control  is  then  transferred  to  one  of  three  major  subroutines,  from  which 
point  the  user  is  free  to  access  different  branches  within  that 
subroutine,  or  request  entrance  into  either  of  the  other  two  major 
subroutines.  The  various  branches  are  described  in  detail  in  the 
following  sections.  Schematics  of  POTEN  and  PEPLT  are  shown  in  Figures 
1,  2  and  3.  The  remainder  of  POTEN  and  PEPLT  consist  of  secondary 
subroutines:  data  read,  physical  properties  of  seawater,  etc.,  which  are 
accessed  as  part  of  the  various  branches  available  to  the  user  in  the 
major  subroutines.  The  file  structure  reflects  the  program  structure 
(Table  1).  POTEN  and  PEPLT  are  linked  by  linking  the  object  files  in 
Table  1.  Accessory  files  are  listed  in  Table  2. 

The  input  data  in  CTD78  disc  format  is  accessed  using  subroutines 
from  CT0ATA/LI8,  and  the  plots  in  PEPLT  are  created  using  the  NCAR  plot 
package.  The  plot  package  creates  a  file  on  logical  unit  8  which  must  be 
read  and  translated  into  plot(s)  by  a  Metacode  translator.  Those 
translators  are  available  both  for  the  high  speed  Calcomp  plotter  and  for 
various  screens,  for  plot  previewing.  The  absolute  plot  dimensions  may 
be  altered  after  the  file  is  created,  and  the  plots  can  be  plotted  as 
many  times  as  desired.  The  use  of  the  translators  is  described  at  the 
end  of  section  3. 

The  multiple  branch  structure  of  the  programs  provides  an  extremely 
powerful  and  flexible  framework  for  computations  which  are  often  not 
routine;  however  useful  documentation  of  such  programs  is  correpondingly 
difficult.  It  is  suggested  that  the  new  user  begin  by  studying  Figs.  1., 
2  and  3.  A  documented  command  file  (ENERGY.COM)  for  a  routine 
computation  and  display  is  found  in  Appendix  A.  This  file  allows  the  new 
user  to  become  familiar  gradually  with  the  options  available  in  the 
programs.  After  studying  and  experimenting  with  the  command  file,  the 
user  may  wish  to  explore  other  options  available  by  referring  to  the 
detailed  branch  descriptions  found  in  sections  2  and  3  of  this  report. 


'Met*' 


STATION  NAME  SHIP, 
CRUISE,  STATION 
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TA8LE  1:  FORTRAN  and  Object  Files 


Main  program 
Major  subroutines 

Data  read  subroutine 
Secondary  subroutines 

Library  subroutines 

System  subroutines  (associated 

with  the  NCAR  plot  package) 


POTEN  Files 

PEPLT  Files 

POTEN 

PEPLT 

PTENS 

PEPLS 

COMPS 

AVRGS 

AVRCP 

TABLE 

DATA 

TDATA 

POTENSUB 

PEPLTSUB 

POTENSUB 

CTDATA/LIB 

AUTOGRAPH 

OASHCHAR 

NCAR 

ft  •»  +-*0**s  r 
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TABLE  2 


Accessory  Fi les 

Common,  dimension  and  equivalence 

required  for  compilation  (FORTRAN) 

Station  lists  (JSHP.PTN  is 

generated  by  COMPS  —  unit  12) 

Input  data 

*.CTD  is  CTD  file 
*.REG  is  generated  by  COMPS 
*.AVG  is  generated  by  AVRCP 
*  is  the  station  identifier 

Common  storage  (binary  file) 

POTEN  and  PEPLT  generate  these 
if  they  do  not  exist  in  the 
directory 

Documentation  files  (formatted) 


Command  file  used  to  set  up 
assignments  and  start  an 
interactive  job 

Command  file  used  to  run  POTEN  and 
PEPLT  sequentially  in  batch  mode, 
supplying  standard  output. 

Command  file  to  create  plots  shown 
in  Fig.  4  using  TABLE  subroutine. 

Command  file  to  compute  dynamic 
height  station  by  station  and 
output  in  map  format. 


POTEN 


COMPOTEN.FOR 


JSHP.PTN 

*  .CTD 
*.REG 


KPTCM.DAT 


POTEN.DOC 


POTEN.COM 


ENERGY.COM 


PEPLT 


COMPEPLT.FOR 


JSHP.PTN 

*.AVG 


KPLCM.DAT 


PEPLS.DOC 

PEPLS3.DOC 

AVRGS.DOC 

PEPLT.COM 


TABLE.COM 


DYNHT.COM 


k.  *»*»***«  '- 
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2.  Documentation  of  POTEN 

In  this  section,  the  major  subroutines  in  POTEN  described  above  are 
documented  in  detail.  They  are  all  structured  around  the  multiple  GO  TO 
statement  of  the  form 

GO  TO  (#1,  #2,  #3 .  n)  BRANCH 

with  #1  through  #n  FORTRAN  statement  numbers,  and  the  BRANCH  #  an  index 
such  that  8RANCH  #  =  n  transfers  control  to  Statement  #n.  The  BRANCH 
is  input  by  the  user  following  an  appropriate  program  prompt.  The  branch 
numbers  are  keyed  to  various  computations  which  may  be  accessed  at  the 
user's  option.  The  keys  are  listed  later  in  this  section,  and  short 
versions  of  them  may  be  obtained  on  the  terminal  any  time  the  program 
prompts  the  user  for  branch  number  input,  by  typing  0/. 

Within  each  branch  there  may  be  options  which  are  accessible  by 
varying  parameters  input  by  the  user  at  the  time  the  branch  number  is 
input.  These  options  are  also  listed  in  the  branch  keys. 

In  addition  to  input  parameter  options,  there  is  an  array  called 
ISSW  with  16  elements  found  in  both  programs.  Within  the  different 
branches,  different  elements  of  ISSW  may  be  tested  for  values  of  -1  or  0, 
and  options  either  accessed  or  skipped  depending  upon  the  value.  In 
general,  ISSW  elements  determine  whether  a  given  type  of  output  is 
generated.  (Historically,  the  ISSW  array  derives  from  the  binary 
switches  available  on  the  shipboard  computer,  the  HP2100  series.)  The 
elements  of  ISSW  may  be  altered  by  accessing  the  appropriate  branch  in 
both  POTEN  and  PEPLT  as  described  below. 

As  described  earlier,  POTEN  is  accessed  as  a  short  main  program 
which  initializes  parameters  if  requested  by  the  user,  or  reads  from  a 
binary  file  KPTCM.DAT  the  most  recently  stored  parameters,  if  no 
initialization  is  requested.  The  main  program  POTEN  then  transfers 
control  to  the  major  subroutine  PTENS  which,  as  shown  in  the  schematic 
Fig.  1,  controls  the  various  branches  available  to  the  user.  PTENS  is 
the  only  component  of  POTEN  in  which  branches  may  be  accessed.  The  two 
remaining  major  subroutines  are  COMPS,  in  which  the  regressions  are 
performed,  and  AVRCP,  in  which  the  horizontal  averaging  is  performed. 
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This  subsection  charts  the  branches  available  to  the  user  in 
detail,  and  describes  briefly  the  working  of  COMPS  and  AVRCP.  Short 
versions  of  the  branch  documentation  are  found  in  the  Appendix,  and  may 
also  be  printed  on  the  screen  while  the  program  is  (interactively)  on 
line  by  typing  0/  whenever  the  program  prompts  for  branch  input. 

2a.  Main  Program: 

The  main  program  queries  the  user  'Initialize  common  (YES  or 
NO)?'.  A  NO  response  causes  the  present  elements  of  KPTCM.DAT,  the 
binary  storage  file,  to  be  read  into  common.  (If  no  file  KPTCM.DAT 
exists  in  the  directory,  the  program  will  create  a  new  file  named 
KPTCM.DAT,  but  if  the  response  to  the  initialization  query  was  NO,  an 
'end  of  file  during  read'  error  will  result.  Therefore,  the  proper 
sequence  of  commands  to  create  a  new  KPTCM.DAT  file  is  to  run  POTEN, 
respond  YES  to  the  initialization,  thereby  creating  a  new  KPTCM.DAT,  but 
not  attempting  to  read  from  it.  Later  in  the  program  (in  branch  3) 
common  may  be  stored  to  the  newly  created  file,  for  use  next  time  the 
program  is  run.)  A  YES  response  initializes  the  data  selection 
parameters  (subroutine  DATA),  and  certain  other  parameters  not  related  to 
the  regressions. 

Following  this  query,  control  is  transferred  to  PTENS,  and  the  user 
is  asked:  'Initialize  regression  parameters  (YES  or  NO)?'.  A  YES 
response  initializes  the  regression  parameters.  A  NO  response  reads  them 
from  KPTCM.DAT.  (Again,  with  a  newly  created  KPTCM.DAT  file,  the  correct 
response  is  YES.)  Finally,  PTENS  asks  for  the  resolution  of  the  input 
data,  before  going  to  branch  mode.  At  this  point  the  user  may  input  up 
to  7  variables,  as  listed  in  the  program  prompt.  The  current  values  of 
the  variables  are  printed  on  the  screen  along  with  the  prompt  list.  The 
variables  are:  KBR,  the  branch  number;  ISW  and  JSW,  which  may  access 
different  options  in  branch  KBR;  KLIST,  usually  the  list  output  logical 
device  number  (reset  to  6  each  time  the  prompt  is  printed);  KOUT  and  KTP, 
the  data  output  and  input  logical  device  numbers  (note  that  the  program 
uses  named  files  for  data  input  and  output  via  OPEN  statements  which  use 
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KTP  and  KOUT  as  unit  numbers  );  KIN,  the  program  input  device  for  screen 
or  command  file.  (Changing  the  value  of  KIN  to  6  part  way  through  a  COM 
file  transfers  control  to  the  screen,  allowing  interactive  mode  —  see 
Appendix  A  for  an  example  command  file,  POTEN.COM.) 

2b.  Branches  (KBR) 

0:  Short  documentation  printed  on  screen.  See  Appendix  B  for  a 
listing  of  this  documentation. 

SUMMARY  -  POTEN.'PTENS:  KBR  =  0 

Function:  List  on  terminal  the  short  documentation  for  PTENS 

ISW,  JSW  options:  None 

Output  device:  unit  KTTX 

Input  device:  None 

ISSW  options:  None 

1:  Calls  COMPS  subroutine,  which  performs  the  following  sequence: 

a.  Calls  subroutine  DATA,  which  opens  the  subindex  directory 
for  the  default  file  specifications  of  the  input  data.  (Those 
specifications  may  be  changed  by  calling  KBR  ■  13,  which  is 
identical  to  branch  1  except  for  allowing  file  specifications  to 
be  changed.)  Then  the  header  for  the  ISWth  sequential  station 
in  that  subindex  file  is  examined  to  see  if  it  meets  data 
selection  criteria.  If  so,  a  file  name  is  written  to  file 
corresponding  to  logical  unit  12.  Throughout  this  report  that 
file  is  called  JSHP.PTN;  an  example  is  given  in  Table  3.  The 
temperature  and  salinity  data  are  transferred  to  array  DATAX, 
using  Millard  subroutine  GETDAT.  Pressure  is  stored  in  the 
zeroth  element  of  DATAX,  which  is  equivalenced  to  array 

PRESS.  PRESS  is  used  throught  COMPS  and  AVRCP.  The  total 
number  of  scans  (NTOT)  is  also  noted.  The  above  occurs  in 
subroutine  DATA,  after  which  control  returns  to  COMPS. 

b.  COMPS  then  sets  up  the  regression  for  the  first  interval 
using  parameters  which  may  be  changed  using  branch  3  and 
continues  the  computation  through  all  the  Intervals  requested. 
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TABLE  3 

Example  of  a  JSHP.PTN  file  generated  by  COMPS 


Consec. 

Station  I.D. 

Weight 

Number 

1 

GY001302 

1.0 

2 

GY001003 

1.0 

3 

GY001004 

1.0 
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or  until  the  end  of  the  data  (determined  by  NTOT)  is  reached. 

For  each  interval  potential  temperature  and  steric  anomaly 
referred  to  p^  (the  level  pressure)  are  calculated  for  each 
data  scan  to  be  used  in  the  regression.  Potential  temperature 
is  calculated  according  to  Fofonoff  (1977),  using  the  polynomial 
formula  of  Bryden  (1973)  for  the  adiabatic  temperature 
gradient.  Steric  anomaly  6  is  calculated  as: 

6  =  105  x  (a(p,e(p,T,S,pf),  5)  -  a(p,0,35) ) 
with  a  the  specific  volume  calculated  according  to  the  SCOR 
Working  Group  51  new  equation  of  state  for  seawater  (Mi llero,  et 
al,  1980),  for  which  an  algorithm  is  given  by  Fofonoff  (1981). 
Within  each  interval  an  editing  process  occurs  in  which  points 
exceeding  three  standard  deviations  of  the  regression  estimate 
at  a  given  steric  anomaly  are  flagged.  Temperature  and  saliity 
are  then  regressed  against  pressure  over  the  interval.  Any 
points  in  T  or  S  which  exceed  three  standard  deviations  are 
replaced  by  the  regression  estimate.  The  regression  of  steric 
anomaly  is  performed  again  and  rechecked.  The  number  of 
standard  deviations  for  both  tests  may  be  changed  —  see 
KBR  =  3.  The  interpolated  scans  are  printed  out  on  unit  KLIST 
and  data  scans  which  are  flagged  but  not  interpolated  are  also 
listed  as  such  on  KLIST  if  ISSW  (3)  is  set  to  -1.  (ISSW  values 
may  be  changed  using  branch  5.)  Pressure  p  and  potential 
temperature  e  (referred  to  the  level  pressure  p^)  are 
regressed  against  steric  volume  anomaly  (also  referred  to 
Pf)  and  the  coefficients  for  both  p  and  e  are  stored  in 
arrays  CP  and  CT  for  each  interval.  Data  output  occurs  if 
ISSW  (13)  =  -1,  and  is  written  into  a  file  with  the  name  *.REG, 
where  *  identifies  the  station,  a  two  character  (alpha)  ship 
name,  a  3  digit  cruise  number  and  a  3  digit  station  number.  The 
format  of  the  ouput  file  is  a  header  of  150  words  equivalenced 
to  an  1*4  array  followed  by  a  variable  number  of  data  records 
(each  46  words,  also  an  1*4  array),  one  record  per  level 
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TABLE  4 

POTEN  Data  Output  Variables 


HEADER  RECORD: 

VARIABLE 

150  WORDS 

NAME 

DESCRIPTION 

LTYPE 

Identifies  record  as  header  record  (LTYPE  =  1) 

MHDR 

Number  of  elements  in  header 

ICON 

Sequential  number  of  station  (in  POTEN  calculation) 

ISHP 

Ship  name  (A2  format) 

KCAST 

Station  number 

I  DAY 

Julian  year  day 

I  PR 

First  pressure 

LPR 

Last  pressure 

XL  AT 

Latitude  of  station 

XLONG 

Longitude  of  station 

WGT 

We i ght 

XLTJD 

Latitude  of  origin  for  distance  computations  in 
kilometers  (negative  for  south  latitude) 

XLG0 

Longitude  of  origin  (negative  for  west) 

LBBL (3) 

Short  station  label  (3A4  formt) 

LBL (13) 

Run  identification  label  (13A4  format) 

NSC (60) 

Regression  parameters  j  $ee  text 

Regression  parameters 

NPR  ( 60 ) 

NSECTION 

Number  of  sections  in  the  water  column 

DATA  RECORD:  46  WORDS 

KTYPE 

Identifies  record  as  data  record  (KTYPE  *  0) 

MBUF 

Number  of  elements  in  data  record 

IREC 

Level  number 

N 

Polynomial  order 

NDP 

Number  of  data  scans  used  in  regression 

KSW 

Not  used 
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TABLE  4  (continued) 


NAME 

LI 

L2 

PF 

T0,S0,DV0 


PI 

THF 

DVI 

DVF 

PM,THM,SM,DVM 

DH 

PE 

XPE 


CP(8) 

Z1 

CT{8) 

Z2 

F1,F2,F3 

XLT0: 

XLG0: 


DESCRIPTION 

Not  used 
Not  used 
Level  pressure 

Temperature,  salinity  and  steric  anomaly  from  input 
data,  averaged  about  PF  ±  POIFF  (see  branch  3 
description  and  Table  5) 

Pressure  of  the  reference  steric  anomaly  (DVF)  in  the 
unleveled  or  initial  field 
Local  potential  temperature  (referred  to  PF)  as 
estimated  by  the  regressions:  e^(P^) 

Steric  anomaly  corresponding  to  PF  in  the  initial  field 
Steric  anomaly  corresponding  to  PF  in  the  leveled  field 
Average  of  pressure,  potential  temperature,  salinity 
and  steric  anomaly  over  the  regression  interval, 
da/dp  based  on  the  averaged  regression  coefficients 
Potential  energy  anomaly  j  Recommended  that 

Horizontally  average  PE  these  not  be  used, 

but  calculated  in 
PEPLT 

Pressure  vs.  steric  anomaly  coefficients 
Standard  deviation  of  regression  pressure  estimate 
(Fofonoff  and  Bryden,  1975) 

Potential  temperature  vs.  steric  anomaly  regression 
coefficients 

Standard  deviation  of  regression  temperature  estimate 
Steric  volume  minimum,  maximum  and  average  over 
regression  interval 

Latitude  of  origin:  default  is  40.0 

Longitude  of  origin:  default  is  -70.0 
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Pf.  The  output  is  in  binary  (unformatted)  files.  The  variables 
output  are  identified  in  Table  4.  Some  information  at  each  level  may 
be  output  to  unit  KLIST  if  ISSW  (12)  *  -1,  for  purposes  of  checking. 
Header  information  is  ouput  to  unit  KLIST  if  ISSW  (11)  =  -1.  The 
input  data  scans  are  output  to  unit  KLIST  if  ISSW  (5)  =  -1  and  the 
regression  coefficients  and  residuals  are  output  to  unit  KLIST  if 
ISSW  (10)  =  -1.  If  ISSW  (6)  =  -1  statistics  of  the  coefficients  are 
printed  on  unit  KLIST.  The  ratio  of  each  coefficient  to  its  standard 
deviation  (see  Fofonoff  and  Bryden,  1975,  Appendix)  is  computed.  For 
an  infinite  number  of  degrees  of  freedom,  at  95%  confidence  that 
ratio  should  equal  or  exceed  1.96.  The  statistic  which  is  listed  is 
(a^  the  coefficients): 

ai 

std  dev  (a.)  ' 


When  stations  with  subindex  reference  number  (sequential  number) 
ISW  through  JSW  have  been  tested  for  data  selection  criteria  and 
either  been  skipped  or  have  gone  through  the  regression 
calculation,  COMPS  returns  control  to  PTENS. 


SUMMARY  -  POTEN: PTENS:  KBR  =  1 

Function:  calls  COMPS  subroutine 

ISW,  JSW  Options:  ISW  to  JSW  are  the  station  reference  numbers 
data  goes  automatically  to  *.REG  file  if 
ISSW  (13)  =  -1;  other  information  output  goes 
to  unit  KLIST,  as  requested  by  elements  of  ISSW 
3  =  -1  Print  out  interpolated  scans 

5  =  -1  Print  out  input  data  scans 

6  =  -1  Print  out  coefficient  statistics 

10  =  -1  Print  out  regression  coefficients  for 
each  scan 

11  =  -1  Print  out  header  information 

12  =  -1  Print  out  selected  data  following 
regression 


Output  device: 


ISSW  options: 
(to  unit  KLIST) 
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13  =  -1  Data  output  to  *.REG 

2:  Initializes  data  selection  parameters  described  in  Table  5 
SUMMARY  -  POTEN:PTENS:  KBR  =  2 

Function:  Initialize  data  selection  parameters 
ISW,  JSW  Options:  None 
Output  device:  None 
ISSW  Options:  None 

3:  Changes  or  lists  regression  and  data  selection  parameters 
described  in  Table  5.  The  data  selection  parameters  are 
straightforward.  For  the  regression  parameters  the  water  column 
is  divided  into  a  maximum  of  nine  sections,  each  of  which  may 
have  a  number  of  levels  whose  regression  parameters  are  the 
same.  The  regression  parameters  consist  of  the  total  number  of 
sections;  in  each  section,  the  interval  between  leveled  surfaces, 
the  interval  over  which  the  regression  is  performed,  the 
polynomial  order,  and  start  and  end  pressures  for  the  section. 

All  of  these  parameters  are  input  using  subroutine  PARAM,  which 
branch  3  calls.  The  prompts  are  (hopefully)  self-explanatory. 
After  parameters  have  been  entered  for  all  sections,  PARAM 
translates  them  into  internal  parameters  which  control  the  way 
the  program  performs  the  regressions.  These  internal  parameters 
are  stored  in  arrays  NPR  and  NSC.  Since  the  arrays  NPR  and  NSC 
are  included  in  common  stored  to  KPTCM.DAT,  the  user  form 
parameters  need  be  entered  only  once,  until  a  change  is 
required.  The  old  parameters  may  be  retrieved  by  responding  'NO' 
to  the  initial  query  in  PTENS  'Initialize  regression 
parameters?'.  Stored  common  is  written  to  KPTCM.DAT  at  the  end 
of  branch  3,  so  any  changes  in  regression  parameters  will 
overwrite  the  most  recent  ones  in  KPTCM.DAT,  provided  branch  3  is 
completed.  It  is  not  possible  to  change  only  a  single  regression 
parameter;  if  a  change  is  required,  all  the  parameters  must  be 
re-entered.  (This  is  because  the  internal  parameters  NPR  and  NSC 
have  elements  whose  value  depends  upon  parameters  for  more  than 
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TABLE  5 

POTEN  Parameters:  Branch  KBR  =  3 


Parameter 

ICON 


KSW 

A2 


A3 


WGT 

PDIFF 

DELP 


Definition 
Consecutive  number 


Not  used 

Number  of  standard  deviations 
allowed  for  a  regression  point 
in  p(a)  before  flagging. 

Numbe  '  of  standard  deviations 
allowed  for  a  regression  point 
in  T(p)  and  S(p)  before 
interpolation 
Weight 

Interval  (db)  about  Pf  for 
averaging  T0,S0,P0 
Pressure  series  interval  for 
input  CTD  data  (db) 


Default  if  Initialized 
1  for  first  station. 
Increments  with  stations 
processed 

1 

3. 


3. 


1. 


6. 


2. 


REGRESSION  Parameters  —  as  described  in  program  prompts 

Data  selection  parameters:  windows  such  that  data  inside  all  windows  is 

accessed;  all  other  data  skipped 


IDAY1 

Minimum  Julian  year  day 

0 

IDAY2 

Maximum  Julian  year  day 

365 

JDO 

Additive  constant  to  actual  day 

0 

XEMN 

Minimum  longitude 

-180.0 

XEMX 

Maximum  longitude 

180.0 

XMN 

Minimum  longitude 

-  90.0 

XNMX 

Maximum  latitude 

90.0 
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one  section.  PARAM  requires  that  parameters  be  input 
sequentially.)  It  is  not  necessary  to  understand  how  NPR  and  NSC 
work  in  order  to  run  the  program  (that  is  the  purpose  of  the 
PARAM  subroutine);  however,  modifications  of  the  program  may 
require  that  the  programmer  know  how  these  arrays  function.  A 
brief  description  is  therefore  presented  here.  The  pressure 
Pf  for  each  level  is  given  by: 

For  IREC  less  than  NPR(section  #) 

PF  =  NPR(section  #  +  total  number  of  sections)  x 
(IREC  -  NPR(section  #  +  2  x  total  number  of  sections)) 

NPR  (section  #  +  total  number  of  sections)  contains  the 

interval  between  pressure  levels;  NPR(section  #  +  2  x  total 
number  of  sections)  contains  an  index  which  allows  the  correct 
Pf  to  be  determined,  while  NPR(section  #  )  contains  the  level 
number  at  which  the  section  commences.  Some  care  should  be  taken 
to  assure  that  the  parameters  input  are  consistent. 

Specifically,  the  first  level  of  a  new  section  must  have  a 
pressure  pf  such  that  p^  is  some  integral  multiple  of  the 
pressure  interval  between  leveled  surfaces  in  that  section.  The 
use  of  the  total  number  of  sections  allows  the  program  to  treat 
NPR  as  a  variable  length  two-dimensional  array,  even  though  it  is 
in  fact  singly  dimensioned.  Subroutine  PARAM  adds  an  additional 
'dummy'  section  below  those  input  by  the  user  to  assure  that 
COMPS  does  not  continue  below  the  desired  depth.  Thus,  the  total 
number  of  sections  (NSECTION)  will  always  be  one  greater  than  the 
number  input  by  the  user. 

Array  NSC  contains  the  remainder  of  the  parameters:  start 
pressure  in  NSC(section#  ),  polynomial  order  in  NSC(NSECTION  + 
section  #),  number  of  data  scans  in  the  regression  interval  in 
NSC(2*NSECTI0N  +  section  #). 

SUMMARY  -  POTEN-.PTENS:  KBR  =  3 

Function:  Change  or  list  regression  and  data  selection  parameters 
ISW,  JSW  Options:  ISW  =  0:  short  list  only 
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ISM  =  1:  full  list 
JSW:  no  options 

Input  device:  unit  KIN 
Output  device:  unit  KLIST 
ISSW  Options:  None 

4:  Call  AVRCP  -  averaging  subroutine.  The  pressure  and  potential 
temperature  coefficients  from  the  regressions  performed  in  COMPS 
are  averaged  horizontally,  level  by  level.  The  average  pressure 
polynomial  at  each  p^  is  set  equal  to  pf 
(corresponding  to  a  mass  conservation  constraint  between  the 
initial  and  leveled  fields)  and  the  resultant  polynomial  is 
inverted  to  obtain  the  reference  steric  anomaly  («^) 
corresponding  to  that  p^.  (See  Bray  and  Fofonoff,  1981  for 
a  more  detailed  discussion.) 

The  averaging  is  actually  done  in  two  'passes'  through  the 
data,  but  a  single  call  to  AVRCP  with  ISSW (7)  *  0  will 
automatically  average  and  output  new  station  data  files  based  on 
the  leveled  field.  (Data  ouput  occurs  if  ISSW(13)  =  -1,  as  in 
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COMPS.  The  new  files  are  called  *.AVG  with  *  as  before  the 
station  identifier.)  Information  about  the  averaged  pressure 
coefficients  is  output  to  unit  KLIST  if  ISSW(12)  =  -1. 

Information  about  the  averaged  steric  field  is  output  to  unit 
KLIST  if  ISSW(ll)  =  -1. 

The  two  averaging  'passes'  may  be  accessed  individually,  and 
separately  from  the  data  ouput  by  setting  ISSW ( 7 )  =  -1  and 
entering  K8R  =  4,  ISW  =  1  for  the  first  pass,  KBR  =  4,  ISW  =  2 
for  the  second  pass  and  KBR  =  4,  ISW  *  3  to  output  the  new 
station  data  files.  However,  since  the  second  pass  must  be 
performed  directly  after  the  first,  and  the  output  directly  after 
the  averaging  it  is  recommended  that  the  automatic  access  be  used 
{ ISSW (7)  =  0).  If  no  output  is  desired,  ISSW(13)  should  be  set 
to  0. 

SUMMARY  -  POT£N:PTENS:  KBR  =  4 

Function:  Call  AVRCP  averaging  subroutine 


ml 
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ISW,  JSW  options:  If  ISSW{ 7)  =  -1 


If  ISSW{7)  =  0: 


ISW  =  1:  First  averaging  pass 
ISW  =  2:  Second  averaging  pass 
ISW  =  3:  Output  of  data  to 
*.AVG  files  if 
ISSW (13)  =  -1 

ISW  =  1:  Averaging  and  output 
performed 
automatical ly. 


Input  files: 
Output  files, 
Output  files, 
ISSW  Options: 


*.REG 

data:  J  .AVG 

lists:  unit  KLIST 

ISSW(7)  =  -1:  individual  access  of  averaging  passes 
ISSW (11)  =  -1:  List  of  averaged  steric  field  on 


unit  KLIST 


ISSW (12)  =  -1:  List  of  averaged  pressure 
coefficients  on  unit  KLIST 
ISSW(I3)  =  -1:  Leveled  field  based  data  output  to 
*.AVG  files 


5:  Not  used 

6:  Print  data  label.  This  label  is  input  by  the  user  in  branch  3, 

and  is  carried  in  both  the  *.REG  and  *.AVG  files  as  an  identifier 
of  the  group  of  stations,  the  version  of  the  POTEN  run,  etc.  Its 
format  is  13  A4  or  a  total  of  52  characters.  Branch  6  lists  this 
label  to  unit  KLIST. 

SUMMARY  -  POTEN :PTENS:  KBR  =  6 


Function:  Write  data  label 
ISW,  JSW  options:  None 
Output  device:  unit  KLIST 

ISSW  options:  None 

7:  Not  used 

8:  Write  header  record  to  unit  KLIST:  Station  label,  position, 

origin,  LTYPE,  MHDR,  ICON,  ISHP,  ICAST,  JDAY,  IPR,  LPR.  This  is 
also  done  automatically  in  subroutine  DATA  when  COMPS  accesses  the 
station,  provided  ISSW (11)  »  -1. 
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SUMMARY  -  POTEN:PTENS:  KBR  =  8 

Function:  Write  station  header  information 
ISW,  JSW  options:  None 
Output  device:  unit  KLIST 

ISSW  options:  None 

9:  Write  *.REG  or  *.AV6  single  data  record  to  unit  KLIST.  Of 

doubtful  usefulness,  this  branch  was  part  of  the  original  program. 
SUMMARY  -  POTEN: PTENS:  KBR  =  9 

Function:  Write  single  output  data  record  to  unit  KLIST 
ISW,  JSW  options:  None 
Output  device:  unit  KLIST 

ISSW  options:  None 

10:  Not  used 

11:  Set  the  values  of  the  ISSW  array.  One  call  allows  up  to  16 

inputs.  Each  input  consists  of  element  number  followed  by  a  comma 
and  the  value  to  assign  to  that  element.  Whenever  input  is 
complete,  if  less  than  16,  the  branch  may  be  terminated  with  a  /. 
SUMMARY  -  POTEN: PTENS:  KBR  =  11 

Function:  Set  ISSW  array 
ISW,  JSW  options:  None 
Output  device:  unit  KTTX 

Input  device:  unit  KIN 

ISSW  options:  None 

12:  Exit  program.  Program  queries  'Exit  program  A  YES  response 
results  in  a  FORTRAN  stop  statement  execution.  A  NO  response 
returns  the  PTENS  branch  prompt. 

SUMMARY  -  POTEN.-PTENS:  KBR  =  12 
Function:  Exit  program 
Input  device:  unit  KIN 

If  a  value  of  KBR  greater  than  12  or  less  than  0  is  entered,  the  short 
documentaton  is  printed  on  the  screen. 


3.  PEPLT  (Documentation 


Like  POTEN,  PEPLT  is  accessed  through  a  short  main  program,  which 
initializes  parameters  as  requested  by  the  user,  and  then  transfers 
control  to  a  major  subroutine,  PEPLS.  From  PEPLS,  the  user  may  call 
subroutine  TABLE,  which  plots  and  lists  station  by  station,  and 
subroutine  AVR6S  which  computes  and  displays  horizontally  averaged 
quantities  as  a  function  of  depth.  Subroutine  AVRGS  has  its  own  set  of 
internal  branches,  one  of  which  returns  program  control  to  PEPLS. 
Subroutine  TABLE  has  no  internal  branches.  As  in  POTEN,  short 
documentation  can  be  displayed  on  the  screen  while  the  program  is  running 
interactively,  by  typing  tl  as  a  response  to  branch  prompts  in  either 
PEPLS  or  AVRGS. 

3a.  Main  Program:  PEPLT 

The  main  program  queries  'Load  in  previously  stored  common?'.  A 
YES  response  causes  the  elements  of  the  binary  array  KPLCM.DAT  to  be  read 
into  common,  beginning  with  the  common  element  KTTX.  a  'NO'  response 
causes  no  action  by  the  program.  Control  is  then  transferred  to 
subroutine  PEPLS. 

3b.  Branches  -  PEPLT 

1:  Calls  subroutine  TABLE.  TABLE  plots  and  lists  station  by 

station.  It  also  outputs  requested  information  in  a  format 
appropriate  as  input  to  objective  mapping  programs.  The  plot 
section  of  TABLE  is  designed  to  permit  a  number  of  stations  to  be 
plotted  on  the  same  frame,  with  the  origin  of  each  station  within 
the  larger  frame.  Examples  are  shown  in  Fig.  4.  In  Fig.  4a  the 
buoyancy  frequency  N  is  plotted  as  a  function  of  geographical 
position  (relative  to  an  origin  at  37*N,  69.65°W),  the 
coordinates  of  the  frame;  and,  for  each  station,  as  a  function  of 
depth,  where  the  station  axes  represent  0  to  3000  db  vertically 
and  -3  to  3  cph  horizontally.  This  is  accomplished  by  scaling 
the  buoyancy  frequency,  and  adding  it  to  the  X-coordinate  (in 
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Figure  4:  Example  plots  from  TABLE.COM 
4a.  Buoyancy  Frequency,  N  in  cph.  Inset  axes  represent  *3  cph 

(horizontal),  0  to  3000  db  in  pressure  (vertical).  The  origin  on  the 
inset  axes  (0,0)  represents  the  station  location. 
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LDE  FS5  VERTICAL  DISPLACEMENTS.  DB 


KM  EAST  0F  ORIGIN 


4b.  Vertical  displacement  *  in  db.  Inset  axes  represent  *100  db 
(horizontal)  and  0  to  3000  db  (vertical). 
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this  example,  kilometers  east  of  the  origin)  and  scaling  the 
pressure  and  subtracting  from  the  Y-coordinate  (here  km  north  of 
the  origin).  The  program  allows  up  to  four  variables  to  be 
included  in  such  a  sum  for  X  and  Y.  The  length  of  the  station 
axes  are  determined  by  the  variables  X2DIM  and  Y2DIM,  which  are 
input  in  user  units  (i.e.  cph  in  Fig.  4a)  for  the  variables  being 
plotted.  In  the  example,  Y2DIM  is  3000  (db)  and  X2DIM  is 
3  (cph).  One  frame  is  created  for  each  call  to  TABLE;  it  will 
encompass  ND  station  plots.  The  plot  parameters  may  be 
initialized  by  calling  TABLE  (KBR  =  1)  with  ISW  >  2.  Control 
then  returns  to  PEPLS.  Plot  parameters  may  be  changed  by  calling 
TABLE  with  ISW  =  1.  Control  again  returns  to  PEPLS.  Plotting 
commences  only  when  TABLE  is  called  with  ISW  =  0.  Figure  4b  is 
the  same  type  of  plot  as  4a,  with  vertical  displacements  plotted 
instead  of  N.  Both  of  these  plots  were  created  using  the 
documented  command  file  TABLE.COM  found  in  Appendix  A. 

A  number  of  variables  relating  to  the  leveled  field,  the  inital 
field,  and  the  location,  and  time  of  each  station  may  be  examined 
using  PEPLT.  A  list  of  these  variables  is  found  in  Table  6;  they 
are  computed  in  function  subroutine  VRBL,  coded  by  number.  Thus 
a  call  to  VRBL  (3)  returns  the  latitude  of  the  station  being 
examined  (variable  XLAT).  Subroutine  TABLE  plots  the  following 
for  x  and  y: 

x  =  A1*VRBL(NX1)  +  A2*VRBL(NX2)  +  A3*VRBL(NX3) 
+  A4*C(IREC,1) 

y  =  B1*VRBL(NY1)  +  B2*VRBL(NY2)  +  B3*VRBL(NY3) 
+  B4*C(IREC,2) 

Here  C(IREC,n)  refers  to  an  array  which  may  be  filled  using  AVRGS 
subroutine  (see  branch  4).  A1  to  A4,  Bl  to  B4,  NX1  to  NX3,  and 
NY1  to  NY3  may  be  changed  by  accessing  branch  3.  The  default 
values  (initialized  by  KBR  =  1,  ISW  =  2)  are: 
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TABLE  6  —  PEPLT  Variables 

(Nomenclature  follows  that  of  Bray  and  Fofonoff,  1981) 


VARIABLE 

NUMBER 

NAME  or 

SYMBOL 

UNITS 

DESCRIPTION 

-1 

1 

None 

Returns  the  number  1  (counts 

number  of  observations  at 

each  level). 

i 

i 

None 

Returns  i 

1 

XPL 

km 

Zonal  distance  from  origin 
(XLT0.XLG0) 

2 

YPL 

km 

Meridional  distance  from 

origin 

3 

XLAT 

degrees 

Signed  decimal  latitude 
(south  negative) 

4 

XL0N6 

degrees 

Signed  decimal  longitude 
(west  negative) 

5 

ICON 

None 

Consecutive  station  number  in 

POTEN  computation 

6 

WGT 

None 

Averaging  weight 

7 

JDAY 

days 

Julian  year  day 

8 

ISHP 

None 

Ship  code 

9 

ICAST 

None 

Station  number 

10 

N 

None 

Polynomial  order 

11 

NOP 

None 

Number  of  data  scans  in 

regression  interval 

12 

PF 

db 

Level  pressure 

13 

4  A 

T0 

1 

T,  S,  6  averaged  over  the 

14 

S0 

ppt  1 

10“^cnr*»gm“^ 

interval  PF  *  PDIFF 

15 

DV0 

16 

PI 

db 

Pressure  of  reference  steric 
anomaly  («^)  in  the 
initial  field 
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TABLE  6  (Continued) 


VARIABLE 

NUMBER 

NAME  or 

SYMBOL 

UNITS 

DESCRIPTION 

17 

THF 

Cc 

Local  potential  temperature 
referred  to  PF.  (See  65.) 

18 

OVI 

10“®cm3*gm-* 

Initial  steric  anomaly 

on  PF 

19 

DVF 

10-^cm^*gm_^ 

Reference  steric  anomaly  (af) 
on  PF 

20 

PM 

db  ' 

\  Pressure,  local  e,  salinity 

21 

THM 

°C 

and  steric  anomaly  referred  to 

22 

SM 

ppt 

1  Pf  averaged  over 

23 

DVM 

10"5cm3*gm_1> 

regression  interval 

24 

a  * 

Pf 

10~5cm3,gm“^db~3 

d«/dp  based  on  averaged 
regression  coefficients 

25 

X 

Potential  energy  anomaly 
calculated  by  POTEN 

26 

APE 

APE  calculated  by  POTEN 

27  to  34 

CP(1)  to  CP(8) 

(db)"1 

Regression  coefficients  for 

pressure 

35 

Z1 

db 

Standard  deviation  of 
pressure  regression  estimate 

36  to  43 

CT(1)  to  CP(8) 

Cc)’1 

Regression  coefficients  for 
local  potential  temperature 

44 

Z2 

'c 

10_5cm3 .gm"1 

Standard  deviation  of  local 

potential  temperature 
regression  estimate 

46,47,48 

F1,F2,F3 

Minimum,  maximum  and  average 
values  of  steric  anomaly  over 
regression  interval 
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TABLE  6  (Continued) 


VARIABLE 

NAME  or 

UNITS 

DESCRIPTION 

NUMBER 

SYMBOL 

48 

»* 

db 

'Boussinesq'  displacements: 

49 

PI* 

db 

"*  =  -(«i  -  «i)/(^)* 

IT*  +  PF 

50 

IT 

db 

Displacement  of  initial  field 

from  reference  field  PI-PF. 

51 

A6 

Positive  implies  downward. 

•i  -  6f 

52 

Vortex  stretching 

db 

„TV  Sin(XLAT) 

PIX  TTHffiw 

53 

APE| 

g 

'Boussinesq'  APE  with  true 

displacements 

1  *  2 
”  ?9  V 

54 

55 

APEg 

g 

gPFaD  -  PE 

Boussinesq  APE  with 

(10-^cm^*gm~^)"^db 

Boussinesq  displacements 

1  *  *2 
-  7g  “p** 

56 

E0 

Inverse  of  local  specific 

volume  gradient  with  pressure 

57 

1/E0 

(lO^cm^gnf^Jdb-* 

Local  specific  volume  gradient 
with  pressure 

58 

N2 

10“®(rad  •  sec-^)2 

Squared  buoyancy  frequency 

59 

6P 

#C  db-1 

Potential  temperature  gradient 

def 

w 

% 

— 

TABLE  6  (Continued) 


VARIABLE 

NAME  or 

UNITS 

DESCRIPTION 

NUMBER 

SYMBOL 

60 

s„ 

ppt  db“* 

Salinity  gradient 

P 

dSf 

dp 

61 

Sf 

ppt 

Salinity  corresponding  to 

0f»  Pf»  6f 

62 

dSf 

de^ 

ppt('C)"1 

Gradient  of  salinity  with 
potential  temperature 

63 

„2 

(db)2 

Squared  displacement  of 
initial  field  from  reference 

field 

64 

N 

cph 

Buoyancy  frequency 

65 

*C 

Potential  temperature  at 

T 

referred  to  zero 

pressure 

66 

a© 

10~3gm-cm"3 

3  -1 

cm  *gm 

Sigma  theta  of  Pf, 

Sf,  9f 

67 

»  4 

af 

Specific  volume  anomaly  in 

reference  field 

1  „  2 

Vertical  gradient  of 

68 

~  7  rk* 

compressibi lity  contribution 

to  GPE 

69 

Not  used 

70 

rk 

10-5cm3-  gm-1*db2 

_  (i£)  (See  Bray  and 

<JjT  'ap'a  ' 

Fofonoff,  1981) 
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TABLE  6  (Continued) 


VARIABLE 

NAME  or 

UNITS 

DESCRIPTION 

NUMBER 

SYMBOL 

72 

&Sf 

ppt 

Salinity  anomaly  from  cubic 
spline  fit  to 
Worthington-Metcalf  and 

Iselin  e-s  curves 

^S(P-f  f  »Sj:) 

(Armi  and  Bray,  1981) 

73 

aSi 

PPt 

Same  as  72  but  using  the 

initial  field 
AS(Pi,e(T(Jf),S((J) 


74 

ef 

r  c)2 

UO"4  J-kg'1)2 

Leveled  field  potential 
temperature  squared.  (Used 
in  calculating  horizontal 
standard  deviation  using 
AVR6S  branch  ISW  =  21.) 

75 

RME 

Random  measurement  error 

(based  on  pressure  error  of 
±5db  magnitude,  temperature 
error  of  ,007*C,  salinity 
error  of  .005  ppt)  for  APE. 
See  Bray  and  Fofonoff,  1981 
Appendix,  for  details  of 

error  calculations. 

V(ap**2/2) 

76 

Not  used 

77 

RFE 

(10"4  J-kg"1)2 

Same  as  78  but  ir 

78 

RFE* 

(lO"4J-kg-1)2 

Random  finestructure  error: 

(based  on  3xZl  as  error  in 
**):  V(op*w*2/2) 
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TABLE  6  (Continued) 


VARIABLE 

NAME  or 

UNITS 

DESCRIPTION 

NUMBER 

SYMBOL 

79 

RME 

(10-4  3-kg'1)2 

(10"4  J-kg-1)2 

Same  as  75  but  pressure  error 
only 

80 

RFEC 

Random  f inestructure  errors  in 

the  vertical  compressibility 
term  (must  be  integrated 
using  AVRGS  branch  ISW  =  17) 

V(  rk*2/2 

81 

RMEC 

(10-4  J-kg'1)2 

Random  measurement  error  in 

the  vertical  compressibility 
term  (pressure  error  only) 

82,83 

Not  used 

10'5cm3-  gm'^db'1 

84 

K 

*v  JC 

Compressibility  (— )a 

85 

Ks 

10'5cm3-  gnf^ppt'1 

Derivative  of  specific  volume 

with  respect  to  salinity; 
temperature  and  pressure  held 

constant: 

P,T 

86 

K  IT 

Contribution  to  GPE  from 

horizontal  gradients  of 
compressibility 

87 

APEB 

10'4  J-kg'1 

‘Boussinesq’  APE  per  unit 
mass  with  true  displacements 

88-90 

Not  used 

10'5cm3  •  gm-1 

91 

«i  -V 

92 

(-©p»)-1 

fC)'1 

93 

Not  used 
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TABLE  6  (Continued) 


j! 


VARIABLE 

NAME  or 

UNITS 

DESCRIPTION 

NUMBER 

SYM80L 

94 

ei2 

Potential  temperature 
corresponding  to  p-  in 

initial  field 

squared 

95 

-(er  ®-j ) /®pir 

None 

(if  ei 

is  in  column  4) 

96 

-(ef-  ef)/epir 

None 

(if  ef 

is  in  column  4) 

97 

ei 

°C 

Local  potential  temperature 
at  5^ 

98 

Not  used 
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A1  =  1 

B1  =  1 

NX1  =  1 

NY1  =  2 

A2  =  2 

B2  =  -.003 

NX2  =  64 

NY2  =  12 

A3  =  0 

B3  =  0 

NX  3  =  0 

NY3  =  0 

£ 

It 

O 

B4  =  0 

These  values  will  cause  the  buoyancy  frequency  in  cph  to  be 
plotted  as  a  function  of  meridional  position  (y-axis),  time 
(x-axis)  and  pressure  (station  axis).  X2DIM  defaults  to  3  (cph) 
and  Y2DIM  to  3000  (db),  resulting  in  station  axes  representing 
*3  cph  for  the  displacements  and  0  to  3000  db  for  the  pressure. 
The  default  number  of  stations  (variable  ND)  is  1  and  may  be 
changed  by  calling  KBR  =  1,  ISW  =  1.  The  plot  information  is 
stored  in  the  file  corresponding  to  unit  8.  It  must  be  read  and 
translated  by  a  Metacode  translator.  PEPLT  may  be  run  on  any 
terminal,  but  the  translators  are  only  i  vail  able  for  graphics 
terminals  and  the  Calcomp  plotter.  See  the  last  part  of  this 
section  for  instructions  on  the  access  of  the  translators.  The 
origin  co-ordinates  may  be  changed  in  PEPLS  branch  2. 

Branch  1  with  ISW  =  0  may  be  used  to  change  PMIN  and  PMAX, 
thereby  selecting  a  range  in  pressure  over  which  data  will  be 
used  (all  other  data  is  excluded),  X2DIM  and  Y2DIM,  described 
above,  JMIN,  the  level  number  corresponding  to  the  pressure  at 
which  the  plot  is  to  start  (this  allows  the  user  to  skip  over 
shallow  points  which  may  have  anomalous  values),  and  various  plot 
parameters.  The  plot  parameters  include  PLABL,  the  overall  plot 
label;  XMIN,  XMAX,  YMIN,  YMAX,  the  axis  limits;  XLABL  and  YLABL 
the  x  and  y-axis  labels,  respectively. 

In  addition  to  plots,  if  ISSW(IO)  =  -1  TABLE  outputs  to  unit 
KOUT  the  following  list  of  variables  in  format  (GF8.3): 

PF,  XPL,  YPL,  ( VRBL ( NV (K ) ),K  =  1,6). 

(See  Table  6  for  descriptions  of  these  variables.) 

If  ISSW(1 2)  =  -1,  a  short  list  of  variables  is  output  to  unit 
KL I  ST:  pressure  (PF),  and  the  variables  x,  y  and  z,  z  given  by 
z  =  C1*VRBL(NZ1)  +  C2*VRBl(NZ2)  +  C3*VRBL(NZ3)  +  C4*C(IREC,3) 


I 


i 


SUMMARY  -  PEPLT:PEPLS:  KBR  =  1 

Function:  Call  subroutine  TABLE  -  multiple  station  plots,  map 
format  output,  lists  by  station. 

ISW,  USW  options:  ISW  =  2  Initialize  plot  parameters 

ISW  *  1  Change  plot,  map 
format  and  list 

parameters 

ISW  =  0  Plot,  list,  map 

format  output 

JSW  No  options 

ISSW  options:  ISSW(5)  =  -1  No  interior  axes  on  plot 
I S SW ( 6 )  =  -1  No  plot 
ISSW ( 10 )  =  -i  List  variables 
I S SW ( 1 2 )  =  -1  List  p,  x,  y,  z. 

2:  Change  data  selection  variables.  Calls  subroutine  to  change  time 

and  space  windows  and  origin  co-ordinates. 

3:  Change  plot  and  list  parameters.  This  branch  prints  a  short 

documentation  on  the  screen  each  time  it  is  called.  Parameters 
which  may  be  changed  and  their  descriptions  are  listed  in 
Table  7.  This  branch  has  internal  branches  1  through  8,  which 
are  prompted  by  '**:  PARAMETERS:  KBR3,  ISW3,  KX,  MV,  MW'.  Only 
KBR3  and  ISW3  have  any  effect  in  this  branch.  KX  is  the  total 
number  of  parameter  input  branches  (5).  To  return  to  PEPLS  from 
branch  3  the  user  must  enter  KBR3  =  1,  ISW3  =  0  followed  by  /. 
This  will  cause  the  new  parameter  values  to  be  written  on  unit 
KLIST,  and  stored  common  to  be  written  to  KPLCM.DAT. 

SUMMARY  -  PEPLT: PEPLS:  KBR  =  3 

Function:  Change  or  list  plot  and  listing  parameters 

ISW,  OSW  options:  None 

Input  device:  KIN 

Output  device:  KLIST 

ISSW  options:  None 


t 

s 
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TABLE  7 

PEPLT:  PEPLS  Branch  3  Parameters 


VARIABLE 


DEFAULT 


■  i 
■> 


NX1 

NX2 

NX3 

NY1 

NY2 


- 1 

NY3 

0 

N21 

25 

■ 

F- 

NZ2 

0 

NZ3 

0  J 

t  . 

>'  • 

Al  B1  Cl 

1.' 

i*  . 

A2-A6,  B2-B6,  C2-C6 

0 

- 

01  to  06 

■  i 

TMIN  to  YT 

None 

fej 

SMIN  to  ST 

None 

DESCRIPTION 


Variable  codes  for  VRBL  used 
in  AVRGS  and  TABLE  computations 


Scaling  factors  used  in 
AVRGS  and  TABLE  computations 


Not  used 
Not  used 
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4:  Calls  AVRGS  subroutine.  This  subroutine  calculates  horizontal 

averages,  allows  operations  such  as  vertical  integration  and 
column  addition,  multiplication,  exponentiation  and  division. 
There  are  23  internal  branches  in  AVRGS,  accessed  with  different 
values  of  ISW{0  to  22).  These  internal  branches  are  described 
below,  with  a  summary  at  the  end  of  each.  As  an  overview,  AVRGS 
reads  the  requested  data  from  *.AVG  files  into  a  two-dimensional 
array  C(100,6).  The  rows  (1  to  100)  correspond  to  the 
pressure  levels  and  the  columns  to  variables  requested  by  the 
user  and  computed  in  function  subroutine  VRBL  (see  PtPLS  branch 
1  for  a  description  of  VRBL).  As  each  successive  station  is 
read,  the  elements  of  C  are  added  to,  forming  sums  of  all  data 
available  at  all  levels.  These  sums  must  then  be  divided  by  the 
total  number  of  observations  at  each  level,  to  obtain  the 
average  values.  For  reasons  of  flexibility,  the  reading/summing 
and  division  are  performed  in  separate  ISW  branches  within 
AVRGS.  Once  the  array  C  is  filled  (one  column  of  which  must  be 
the  number  of  observations)  and  averaged,  then  a  number  of 
operations  can  be  performed  on  the  averages.  The  remaining  ISW 
branches  of  AVRGS  are  devoted  to  these  operations. 

AVRGS  has  its  own  prompt  'AVRGS:KBR, ISW, JSW,KLIST' ,  and 
control  does  not  return  to  PEPLS  unless  KBR  =  4  ISW  =  12  is 
accessed.  Therefore,  only  four  variables  (or  <  4  followed  by  a 
/)  need  be  input  following  the  AVRGS  prompt.  In  order  to  keep 
track  of  the  operations  performed  in  AVRGS,  if  ISSW ( 2 )  =  -1  the 
four  parameters  are  written  to  unit  4  each  time  an  AVRGS  branch 
is  accessed,  along  with  other  pertinent  information.  This  ISSW 
option  will  not  be  noted  in  the  summaries. 

Branches  in  AVRGS:  (ISW) 

ISW  =  0:  Prints  short  documentation  on  unit  KTTX 
ISW  =  1:  Reads  station  data  into  C  array.  Variables  corresponding  to 
NV(JSW)  to  NV(KLIST)  (maximum  of  six)  are  read  into  columns 


JSW  to  KLIST  of  array  C  for  ND  number  of  stations  from  file 
JSHP.PTN  (logical  unit  12),  starting  with  the  first  station  in 
that  file.  All  data  between  PMIN  and  PMAX  is  accessed  for  each 
station.  The  array  C  is  stored  to  KPLCM.DAT  before  returning  to 
the  AVR6S  prompt.  If  ISSW(15)  =  -1,  the  weights  (WT)  from 
JSHP.PTN  file  are  used;  otherwise  a  weight  of  1.  is  used.  Each 
element  of  C  is  a  sum  of 

C( IREC, I )  =  C( IREC, I )  +  D ( I ) *WT* ( AV*VRBL ( NV ( I ) ) ) 

+  { BV+C V*VRBL ( NV ( I ) )*VRBL ( NX ( I ) ) ) 

The  default  parameters  are  set  such  that 
C( IREC, I )  =  C ( IREC, I )  +  WT*VRBL(NV(I) ) 

Some  of  the  parameters  used  by  this  branch  may  be  changed  in 
branch  3  of  AVRGS,  and  some  in  branch  3  of  PEPLS. 

SUMMARY  -  PEPLT :AVRGS:  KBR  =  4:  ISW  =  1 

Function:  Read  and  store  data  to  C  array 
JSW,  KLIST  options:  JSW  is  first  column,  KLIST  last  column 
Output  device:  Array  is  stored  to  KPLCM.DAT  for  emergency 

retrieval.  No  other  output. 

ISSW  options:  None 

ISW  =  2:  Zeros  columns  JSW  to  KLIST  of  array  C 

ISW  =  3:  Changes  or  lists  parameters.  Parameters  involved  are  listed 
in  Table  8.  JSW  =  1  initializes  the  parameters  (defaults 
also  in  Table  9)  before  allowing  changes;  JSW  =  0  retains 
previous  values.  (The  first  access  to  this  branch  must 
initialize. ) 

SUMMARY  -  PEPLT: AVRGS:  KBR  *  4:  ISW  =  3 
Function:  Change  parameters 
JSW,  KLIST  option:  JSW  »  1  initializes 
JSW  =  0  prints  current  values 
Output  device:  KTTX 

ISSW  options:  None 

ISW  =  4:  Average  table:  divide  columns  JSW  to  KLIST  by  column  6, 

which  should  have  the  number  of  observations  at  each  level. 
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TABLE  8 

PEPLT:  AVR6S  Branch  ISW  =  3  Parameters 


VARIABLE  DEFAULT  DESCRIPTION 

( Initialized) 


ND 

1 

Number  of  stations  to  be  processed 

NV(1) 

51 

Variables  to  compute  for  C  array  as 

NV(2) 

68 

VRBL(NV(I ) )  in  column  I. 

NV(3) 

86 

See  Table  6  for  VRBL  codes. 

NV(4) 

87 

NV(5) 

63 

NV{6) 

-1 

JREF 

50 

Number  of  levels  to  be  calculated 

JMAX 

55 

Level  number  corresponding  to  reference 
pressure  for  integrations  over  pressure 

NX (I), I =1,6 

0 

Optional  additive  quantities  in  C  array 
element  calculation  (see  text). 

A1 

1.’ 

A2 

0. 

X  and  Y  scaling  factors  for  plots 

A3 

0. 

► 

B1 

1.1 

Initialized  when  PEPLS  is  called  by 

B2 

0-J 

responding  YES  to  'Initialize  common  ' 

NX1 

12' 

Optional  plot  parameters  (see  text  for 

NX2 

0 

AVRGS  branch  ISW  =  7). 

NY1 

19 

NY2 

Initialized  in  PEPLS  as  above 
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ISW  =  5:  Add  column  JSW  vertically,  starting  from  level  2  and  going 
to  JMAX: 

C( IREC, JSW)  *  C ( I REC  -  l.JSW) 

ISW  *  6:  List  C  array  to  unit  KLIST.  Inlcudes  data  label, 
parameters,  level  number  and  pressure,  and  C  array. 

ISW  =  7:  Plot  one  frame.  Up  to  six  curves  allowed  per  frame.  NCAR 
plot  package  outputs  to  unit  8  a  file  which  must  be  read  and 
translated  into  a  plot  by  a  Metacode  translator.  PEPLT  may 
be  run  on  any  terminal,  but  the  plot  files  may  only  be 
translated  on  graphics  terminals  and  the  Calcomp  plotter. 
Instructions  for  running  the  translators  are  found  at  the 
end  of  this  section  of  the  report.  The  plot  branch  asks  for 
the  number  of  curves  {default  1,  maximum  6),  the  level 
number  for  the  first  point  (default  1),  the  plot  label,  the 
minimum  and  maximum  coordinates  for  x  and  y  (unless  the  user 
opts  to  have  the  NCAR  plot  package  compute  the  scales,  by 
responding  YES  to  the  query  'Use  default  axis  parameters?'), 
x  and  y  axis  labels,  and  the  column  number  to  be  plotted. 

The  program  actually  plots: 

x  =  B1*C( J, JSW)  +  B2*C( J,NX2)  +  B3*PF 
y  =  A1*PF  +  A2*C( J,NV1)  +  A3*C(J,NY2) 

The  default  values  of  the  parameters  plots 
C(J,JSW)  vs  PF  (pressure). 

However,  if  for  example  the  user  wished  to  plot  potential 
temperature  e  vs  salinity  S,  with  &  (VRBL(65))  in  column 
1  and  S(VRBL(61))  in  column  2,  then  the  values  of  the  above 
parameters  should  be  changed  (using  AVRGS  branch  3) 

A1  =  0.  B1  =  1  NY1  =  1 

A2  =  1.  82  =  0 

A3  =0.  B3  *  0. 

The  y-axis  runs  backwards  (maximum  at  the  bottom  to  minimum 
at  the  top)  unless  A1  is  equal  to  0.  An  example  is  given  in 
ENERGY.COM  —  see  Appendix  A. 
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Characters  of  the  user's  choice  which  mark  the  actual  data 
points  may  also  be  plotted  if  ISSW(5)  =»  -1.  Note  should  be 
made  that  these  are  not  centered  characters,  so  that  the 
data  point  actually  occurs  wherever  the  plotter  commences 
drawing  the  character. 

SUMMARY  -  PEPLT:AVRGS:  KBR  =  4,  ISW  =  7 

Function:  Plot  one  frame  containing  up  to  six  curves. 

JSW  option:  JSW  is  the  column  number  to  be  plotted.  It 

may  be  changed  while  in  the  plotting  branch. 
Output  device:  Plot  information  goes  to  Metacode  file,  unit  8. 

ISSW  options:  ISSW(5)  *  -1  plots  character  to  mark  actual 

data  points.  Character  is 
requested  while  in  plot  branch. 

ISW  =  8:  Calculates  gravitational  available  potential  energy  per  unit 
mass  (GPE )  and  per  unit  area  (TGPE),  from  the  horizontal 
averaged  steric  volume  WT  ( VRBL (18) )  in  column  1  and  for 
the  reference  steric  volume  DVF  (VRBL (19) )  in  column  2, 
except  for  a  constant  of  integration.  GPE  and  TGPE  relative 
to  some  reference  pressure  are  calculated  by  subtracting 
from  GPE  and  TGPE  at  each  level  the  value  at  the  level 
corresponding  to  the  desired  reference  pressure  (denoted  by 
level  number  JREF)  in  AVRGS  branch  ISW  =  10.  GPE  is  stored 
in  column  1,  TGPE  in  column  2.  The  units  are  10~^ 

J.kg-*  and  10+^  J.m~^,  respectively. 

SUMMARY  -  PEPLT:AVRGS:  KBR  »  4,  ISW  =  8 

Function:  Calculate  GPE  and  TGPE  except  for  a  constant  of 
integration 
JSW  options:  None 

Output  device:  None  (GPE  and  TGPE  replace  DVI  and  DVF  in 

columns  1  and  2,  respectively,  of  array  C.) 

ISSW  options:  None 

ISW  *  9:  Integrate  over  pressure  columns  JSW  to  KLIST.  This  is  an 
alternate  method  for  calculating  GPE  and  subsequently  TGPE, 
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with  DVT  -  DVF  ( VR 8L ( 5 1 ) )  in  column  JSW.  It  may  also  be 
used  to  compute  the  compressibility  effects  in  the  6PE 
calculation  (see  equation  28  in  Bray  and  Fofonoff,  1981). 

The  integration  is  performed  starting  with  the  first 
element  in  the  column,  and  continuing  to  the  last;  the 
reference  value  must  be  subtracted  in  a  separate  operation, 
using  AVRGS  branch  ISW  «  10. 

SUMMARY  -  PEPLT:AVRGS:  KBR  =  4,  ISW  =  9 

Function:  Integration  over  pressure  (except  for  a  constant)  of 
columns  JSW  to  KLIST 

JSW  options:  Columns  JSW  to  KLIST  are  integrated 

Output  device:  None 

ISSW  options:  None 

ISW  *  10:  Subtract  value  at  reference  pressure  (level  corresponding  to 
JREF )  from  all  other  elements  in  columns  JSW  to  KLIST 
Output  device:  None 

ISSW  option:  None 

ISW  =  11:  Add  up  to  four  scaled  columns,  according  to 
J  =  IREC 

C( J.JCl)  =  CR1*C( J, JC1)  +  CR2*C( J, JC2)  +  CR3*C(J,JC3) 

+  CR4*C( JREF, JC4) 

If  JSW  *  1,  JC1,  CRI  to  JC4,  CR4  are  entered;  no  addition  is 
performed. 

If  JSW  =  0,  addition  is  performed  using  most  recently  input 
parameters. 

SUMMARY  -  PEPLT: AVRGS:  KBR  =  4,  ISW  =  11 

Function:  Add  up  to  four  scaled  columns,  row  by  row 
JSW  Option:  0:  perform  addition 

1:  input  scaling  and  column  parameters 
Output  device:  None 

ISSW  option:  None 

ISW  »  12:  Return  to  PEPLS 
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ISM  =  13:  Multiply  up  to  three  scaled  columns,  row  by  row  according  to 
C(IREC.I)  =  C0N1*C(REC, I )*{ C0N2*C( IREC, J)*[C0N3*C( IREC,K)]} 
If  I  =  -1  no  operation  is  performed. 

If  J  =  -1  then  the  expression  in  {}  is  set  to  1;  if 
K  =  -1,  the  expression  in  [  ]  is  set  to  one,  allowing  one, 
two  or  three  scaled  columns  to  be  multiplied  together.  The 
parameters  may  be  changed  when  the  branch  is  accessed.  The 
default  values  are  I,J,K  =  -1;  CONI,  C0N2,  C0N3  =  1. 

SUMMARY  -  PEPLTrAVRGS:  KBR  =  4,  ISM  =  13 


Function:  Multiply  up  to  three  columns,  row  by  row 

JSW  options: 

None 

Output  device: 

None 

ISSW  options: 

None 

ISM  =  14:  Output  in  map  format  to  unit  KTO.  Branch  requests  output 

file  name  and  level  number  (JREC)  desired.  Variables  output 
are: 

IDSTN  (station  identifier:  ship,  station),  XLAT,  XLONG, 

( VR8L ( NV( K ) ) ,  K  =  1,3),  (C(JREC,K),K  «  4,5) 
in  format  (1H  ,A5,2(F8.2),5F(8.3)). 

SUMMARY  -  PEPLT :AVRGS:  KBR  =  4,  ISM  =  14 
Function:  Output  in  map  format 
OSW  option:  None 

Output  device:  Unit  KTO  (may  be  changed  in  branch;  default 

is  60) 

ISSW  options:  None 

ISW  *  15:  Not  used 

ISM  =  16:  Take  any  single  column  to  any  power,  row  by  row.  Operations 
are  performed  on  the  absolute  value  of  all  elements.  If 
JSW  =>  1,  exponent  and  column  inputs  are  prompted.  If 
JSW  =  0,  exponentiation  is  performed.  The  call  to  JSW  *  0 
should  immediately  follow  that  to  JSW  »  1,  as  the  variables 
used  for  exponent  and  column  number  are  not  unique  to  this 
branch. 
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SUMMARY 

ISW 


ISW 

ISW 

ISW 

ISW 

ISW 


-  PEPLT.'AVRGS:  KBR  .  4,  ISW  =  16 
Function:  Exponentiation  of  a  single  column 
JSW  options:  JSW  =  0:  Operation  performed 

JSW  =  1:  Exponent  and  column  entered 
Output  device:  None 

ISSW  options:  None 

=  17:  Integration  of  error  terms:  interval  pressure  squared  as 
the  integration  variable.  This  is  intended  for  the 
calculation  of  measurement  and  f inestructure  errors  in  GPE 
and  TGPE;  as  such  it  may  be  used  on  columns  containing 
averaged  values  of  VRBL  (75  and  77  through  81)  —  see 
Table  6.  This  branch  uses  the  same  algorithm  as  AVRGS 
branch  ISW  *  9,  with  aP^  instead  of  aP  as  the 
integration  variable.  See  AVRGS  branch  ISW  =  9  for  a 
summary. 

=  18:  Writes  into  column  5  the  difference  in  pressure  between  each 
pair  of  levels,  beginning  at  the  top. 

a  19:  Exchange  columns  JSW  and  KLIST. 

*  20:  Input  a  new  single  element  of  C.  Branch  prompts  for  column 
and  row  of  element  to  be  changed. 

a  21:  Compute  the  standard  deviation  and  store  in  column  1  of  any 
quantity  X  for  which  X  (the  average  value)  is  stored  in 
column  4  and  Xc  in  column  3. 

=  22:  Compute  the  dynamic  height  for  each  station  at  any  range  of 
levels  referred  to  level  JREF  and  output  in  map  format. 
Branch  prompts  for  output  device  (default  is  60),  and  level 
numbers  (JREC1,  JREC2)  for  dynamic  height  calculation. 
Reference  level  JREF  may  be  changed  in  AVRGS  branch 
ISW  =  3.  To  calculate  dynamic  height  NV(1)  must  be  18, 

NV ( 2)  19.  Variables  output  are: 

I DSTN  (station  identifier),  XLAT,  XLONG,  Dynamic  height 
(in  dynamic  centimeters),  (NV(K),K  *  3,6). 

Output  occurs  for  ND  stations,  beginning  with  the  first 
station  in  JSHP.PTN  (unit  12). 
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SUMMARY  -  PEPLT:AVRGS:  KBR  *  4,  ISW  =  22 

Function:  Compute  dynamic  height  relative  to  JREF  for  any  range 
of  pressure,  for  each  of  ND  stations  and  output  in  map 
format.  Four  optional  variables  are  also  output,  for 
the  same  range  of  pressure.  An  example  command  file, 
DYNHT.COM  is  found  in  Appendix  A. 

ISW  options:  None 

Output  device:  Unit  KTO  (default  60;  may  be  changed  by  the 

user  when  the  branch  is  accessed). 

ISSW  options:  None 


PEPLT  Branches  (KBR),  continued 

5:  Set  values  of  elements  in  the  ISSW  array.  Up  to  16  inputs  are 

allowed,  each  consisting  of  the  element  number  followed  by  the 
element  value  (-1  or  0).  Terminate  before  16  by  typing  /. 

6:  Restart  main  program. 

7:  Exit  program:  a  YES  response  to  the  branch  query  ’EXIT  PROGRAM  ' 

results  in  the  execution  of  a  FORTRAN  stop.  A  NO  response 
returns  the  PEPLS  prompt. 


Metacode  Translators 

The  translators  for  the  plot  files  (written  to  unit  8)  created  in 
AVRGS  (branch  ISW  *  7)  and  TABLE  (PEPLS  branch  KBR  *  1)  are  device 
specific.  That  is,  each  graphics  terminal  has  its  own  version.  The 
CALCOMP  high  speed  plotter  has  two  versions:  one  with  default  plotting 
parameters,  and  one  which  allows  the  user  to  enlarge  or  stretch  the 
plots,  alter  their  distribution  on  the  plotter  paper,  etc.  The  IMLAC  and 
Tektronix  terminals  also  have  versions  of  the  translator  to  allow  plot 
previewing. 

For  all  translators: 

If  the  plot  file  was  written  to  any  other  file  than  that  named 
F0R008.DAT  (via  an  ASSIGN  statement  before  running  PEPLT)  then  you  must 
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assign  that  ouput  file  name  to  unit  8  before  running  the  translators. 

For  example,  if  your  plot  file  is  named  PLOT.PPT,  you  must  make  the 
following  assignment: 

ASSIGN  PLOT.PPT  F0R008. 

For  the  CALCOMP  (both  versions)  you  must  also  assign  terminal 
TTA4:  to  FOR 061: 

ASSIGN  TTA4:  F0R061 

Then 

RUN  MCTRNPLOT  (for  MetaCode  TRaNslator  PLOT) 
plots  with  default  parameters,  and 
RUN  MCTRNPLT2 

prompts  the  user  for  changes  in  the  plotting  parameters  before  executing 
the  plots.  MCTRNPLT2  asks  three  questions:  first,  how  many  plots  in  the 
y-direction  (across  plotter)?  The  default  is  1,  and  is  retained  if  a  / 
is  entered.  Second,  what  size  shall  the  plots  be?  The  default  is  10  by 
10  inches.  The  new  dimensions  are  entered  in  inches,  and  need  not  be 
equal  for  x  and  y.  Again  a  /  retains  the  default  values.  Finally,  the 
program  asks  for  the  distance  between  plots,  in  inches.  The  default  is  2 
inches  in  both  x  and  y.  All  plots  in  the  file  assigned  to  unit  8  are 
plotted,  sequentially. 

For  the  Tektronix  (or  the  IMLAC  in  Tektronix  mode): 

RUN  MCTRNTEK 

starts  the  plot  previewer.  If  there  is  more  than  one  plot,  the  program 
prompts  for  continuing  to  the  next  plot  by  asking  'Option  ?'  to  which  the 
user  should  respond  C  for  continue,  until  all  plots  in  the  file  assigned 
to  unit  8  have  been  plotted. 

For  the  IMLAC  (recommended  over  the  IMLAC  in  Tektronix  mode,  since 
it  is  simpler,  and  uses  more  of  the  screen): 

RUN  MCTRNDYN1 

starts  the  plot  previewer.  This  program  also  prompts  for  continuation  if 
there  is  more  than  one  plot. 

This  translator  information  is  accurate  as  of  December  1980.  If 
you  encounter  difficulties  you  should  refer  to  the  current  VAX  manual. 


4.  Modification  of  POTEN  to  accept  input  CTD  data  in  other  than  CTD78 
disc  format. 

This  section  is  intended  as  a  guide  to  assist  users  who  wish  to  use 
POTEN  on  CTO  data  with  formats  other  than  that  read  by  the  standard 
version.  In  this  section  the  header  information  required  by  POTEN  is 
described  in  detail,  and  the  procedure  for  reading  data  is  explained. 

The  only  subroutine  which  must  be  changed  is  DATA,  providing  that  the 
input  data  is  an  even  series  in  pressure  with  no  gaps. 


DATA  requires  the  following  header  information  for  each  station: 


Description 

Variable  Name 

Format 

Ship  Name 

ISHP 

A2 

Cruise 

ICRUIS 

A3 

Station 

I  STAS 

13 

Decimal  Latitude  (south  negative) 

XLAT 

F 

Decimal  Longitude  (west  negative) 

XLONG 

F 

Day 

IDA 

12 

Month 

I  MO 

12 

Year  (last  two  digits) 

IYR 

12 

Time  (24  hour  clock) 

ISTME 

14 

Station  Label 

LBBL(3) 

3A4 

Minimum  Pressure 

PM IN  or  IPR 

F  or  I 

Maximum  Pressure 

LPR 

I 

The  CTDATA  library  subroutines  not 

disc  version  of  CTD78  are: 

needed  for  formats  different 

from  the 

PVER 

CRUISE 

Header  Information 

STATION 

0ATI0X 

Data  Retrieval 

GET DAT 

Also,  the  common  file  IOXREC.DIM  should  not  be  included  in  DATA  — 
see  the  statement  INCLUDE  'IOXREC.DIM'.  The  variable  LLREC  is  the  total 
number  of  stations  in  the  subindex  directory;  all  statements  in  DATA  and 
COMPS  which  refer  to  LLREC  may  be  deleted.  The  data  are  stored  in  arrays 
PRESS  and  DATAX. 


Pressure  is  stored  in  PRESS, (  #),  temperature  in  DATAX(1,  #), 
salinity  in  DATAX(2,  #)  with  #  the  data  scan  number.  Subroutine  DATA 
must  fill  DATAX  and  PRESS  (all  scans)  when  it  is  called  for  each 
station.  Finally,  DATA  must  return  to  COMPS  the  total  number  of  data 
scans,  JRMAX. 

Stations  are  selected  by  the  call  to  DATA  in  COMPS.  The  call  is 
CALL  DATA  (KST.l) 

In  COMPS,  KST  is  the  sequential  number  in  the  DO  loop  from  ISW  to 
JSW  in  branch  1  (or  13).  If  the  input  data  is  on  magnetic  tape,  the  user 
may  wish  to  change  the  DO  loop  in  COMPS  to  go  from  1  to  JSW:  that  is, 
start  at  the  beginning  of  the  tape  and  read  through  ISW  stations. 

The  section  of  DATA  in  which  the  ship  and  cruise  specification  may 
be  changed  (NSW  =  2)  can  be  readily  modified  to  accept  similar 
information  (in  branch  13)  pertinent  to  the  user's  input  data. 

The  header  information  should  be  read  in  following  statement  #5, 
replacing  the  statements  between  #5  and  #54.  The  data  should  be  read 
in  in  statements  which  replace  the  calls  to  DATIDX  and  6ETDAT. 
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Appendix  A 

Example  command  files  for  different  routine 
calculations  using  POTEN  and  PEPLT. 

In  this  appendix  documented  command  files  which  perform  various 
routine  calculations  are  listed.  In  the  order  in  which  they  appear  they 
are:  ENERGY.COM,  POTEN.COM,  PEPLT.COM,  TABLE.COM,  and  DYNHT.COM.  Brief 
descriptions  of  these  files  are  also  found  in  Table  2.  The  files 
themselves  contain  detailed  documentation.  Example  plots  from  ENERGY.COM 
are  also  included:  see  Figures  5  to  12.  Example  plots  from  TABLE.COM 
are  found  in  Figure  4. 
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5  1!  *************  FNERGY.COM  *♦♦♦***'*♦*#*** 

10  i!COM«ANO  FILE  TO  COMPUTE  APE  FROM  CTD78  FORMAT  DATA.  CREATES 

— i-g - TTTnremnm  printout  amp  plots." nm novo ~vs  AWinriN  learning 

26  »M0  USF  THE  VAX-11  PROGRAMS  POTEN  AMO  PEPLT. 

34  »!*  JULY  Ml.  NAN  BRAY 

42 

SO  iSFT  OFFAULT  Oft  A  2 :  <  H6316  .LDEV 

100  iASSrCN  JSHPE5.PTN  F0R012  ’FILE  TO  BE  CREATED  BY  POTEN  CONTAINING 
T*6  _  STATION  IDENTIFIERS. 

200  tASSFON  PRINT. PTN  FORGO*  ! F ILE  FOR  LINEPRINTER  OUTPUT  FROM  POTEN 

300  RUN/MODFB  POTEN 

400  YES  INITIALIZE  COMMON 

500  NO  DO  NOT  INITIALIZE  REGRESSION  PARAMETERS 
600  2.5  I  PRESSURE  INTERVAL  FOR  INPUT  DATA 

'  TOO - TTTTT7R7 - MET  PARAMETERS 'FOR  REGRESSION - 

000  *•«••?./  MET  PO IFF  TO  20B — INTERVAL  OVER  WHICH  TO,SO,DVO  AVFRAGFO 

000  /  !Dn  NOT  CHANGE  OELP— PRES  SURE  INTERVAL  FOR  INPUT  DATA 

1000  NO  DO  NOT  CHANGE  REGRESSION  PARAMETERS 
TTOO  /  !LTAVF  time  WINDOW  AT  DEFAULT  0-7650AY5 
1200  /  ! LEAVE  FAST-WEST  SPACE  WINDOW  AT  DEFAULT:  -1*0,180  DEG 

TTOO - 7 TCriVE  ~  WOTTEff =STTOTH  S  P  ATE B  ENDOW  AT  WE11JUTT'^D0  ,^0  TTfC 

1400  YFS  CHANGE  DATA  LABEL 

1500  LDF  ES5  TEST— STANDARD  VFRSION,  NEW  EOT--??  JULY  81 

1600  11/  MET  ISSW  ARRAY 

TBOO  11,-1/  !L  1ST  STATION  TNFOIMATTlW^TTrWTTnTLTTT  IErLF  FRTNT.TTM) 

1 850  11/ 

t700 - 17,0/ - ron  NOT  ri  ST  RFGRFSSrON "SUHNART'RT  EACH  DEPTH  rOR'~E~ACTT  STATTDN 

1950  11/ 

7000  17,-1/  *  CREATF  *. R EG, *.TVC  FILES 

2050  11/ 

7100  5,0/  TOO  NOT  LIST  INDlVinTJAL  INPUT  DATA  SCANS’ 

7150  11/ 


1YJ  ,0/  Tun  NTJ  T  L  111  REGRESSION  CUFF  El CTt NTS  EUR  EACH  LEVEL 
2500  1, 229,236, */!USlNG  OFF  CRUISE  SPECS,  REGRESS  STA  REF  »  1SW  TO  J5W 

7660  17,701,708,*/  ’CHANCE  THE  CRUISE  SPECSS  THEM  AS  IN  PREVIOUS  COMMAND 

2675  W  ! SUBDIRECTORY  VERSION 

7690  T SCO TOOT  ’SHIP,  CRUISE,  PROJECT  NUMBER - 

7900  11/  ’RESET  ISSW  ARRAY 

3000 - r?v-t / — ’LIST  AVERAGED  REGRTSSTON  COEETTCTENTS - 

3100  *,1,0*4/  'PROCEED  THROUGH  ENTIRE  AVERAGING  PROCESS.  LISTS  TO  PRINT. PTN 

7700  t 7/  7FXIT  PROGRAM? 

3300  YFS 

"3*00  S ASSIGN  RRINT.RRT  FOR 004  ITT STING  FTLE  FOR  MERIT  OUTPUT 
3500  5ASST0N  PLOT.PPT  EOROOB  !NCAR  PLOT  FILE  FOR  PFPLT  PLOTS 

1600 - RDN/NnOEB  PFPLT - 

3700  NO  DO  NOT  PFAO  IN  PREVIOUSLY  STQREO  COMMON 
7000  YES  INITIALIZE  DATA  SELECTION  PARAMETERS 

3900  5/  MET  TSSW  ARRAY 

4000  7,-1/  'LIST  OPERATIONS  PERFORMED*  TN  PRINTS PPT 

*050  5/ 

4100 - OtM-7 - ’PLOT  'CHARACTERS  ON  PLOTS  FOR  t  DENT  TOt  CAT  ION - 

*200  *,2,1,6/  *  ZERO  C  ARRAY  IN  AVRGS.  PROGRAM  CONTROL  NOW  IN  AVRGS. 

*700  *33,  TT  MET  VARIABLE  SELECTION  EARA  METER  S' 

**00  1000, V,,,, ,50/  ’CHANGE  1000  TO  «  DE  STATIONS  IF  NO  <  ALL 

*500  /  *00  NOT  CHANGE  AV  THROUGH  NX  ft  1  — 

*600  /  ! DO  NOT  CHANGE  Al  THROUGH  B3 

4 TOD - 7 - TOO  NOT  CHANCF  NT1  YKRCOGH' NY7 - 

*800  *,1,W6  'READ  VARIABLES  AS  SELECTED  INTO  C  ARRAY. 

*700  *,*,r,4  ’OTVTOE  COLUMNS  1-5  BT  FT  NUMBER  ~0F~  7TATT0NS3  TO  AVER  AGE 

5000  *,6,t,*  ! MR ITE  C  ARRAY  TO  PRINT. PPT 

MOO  *,7,1*3  ’INTEGRATE  COLUMNS  t-7  WTTH  RESPFCT  TO  PRESSURE ~  -  - 


T7U0 

5325 

6400 

5500 

•>(,00 

3700 

5800 

3900 

6TJ00 

6100 

6200 

6300 

6600 

6300 

"W 

6700 

6600 

6900 

7000 

7100 


7900 

8000 

8100 

3200 

8300 


8300 

8600 

8  700 
8300 
3900 


9100 

9700 

9300 

9600 

9300 


9700 
9800 
99 00 
10000 
10100 


IK/  e  u  v 

10300 

ro6oo 

10300 

10600 

10700 


1 1 /nxjK/ 

10900 

irooo 

moo 

11200 


4.10.1.3  'SUBTRACT  FRO*  ALL  LFVELS  THF  VALUE  AT  LEVEL  JRFF 

4,6, 1*4  'WRITE  C  ARRAY  TO  PRINT.PPT  _ 

4,6, 1*4  'WRITE  C  ARRAY  TO  PRINT.PPT* 

4,11*1/  'RESET  AODIWE  CONSTANTS 

1,1. »?*!., 3,1. ,4,0.  'REPLACE  COMPRESSIBILITY  TER8S 
4 , 1 1 *07 
4,16*1  / 

".■6TS"-  TTART  C0CTT8N  8  TTrTHF"  PTTWFR  T5 -  - - 

4,16,0/ 

4,6, 1*4  'WRITE  C  ARRAY  TO  PRINT. PPT  - 
4,7/  fCALL  PLOT  BRANCH 

4.3  «4  PLOTS  IN  THIS  FRAME,  STARTING'  AT  LFVEL  3  ON  EACH 

YES  INPUT  NEW  PLOT  LABEL  _ 

WO  00  NO  USE  OFEAULT  AXIS  PARAMETERS 
-20,200,0,3000  !  XW IN yXMAX ,YMlN,YMAX 

YES  CHANGE  X-AXIS  LABEL 
APE  10M/SFC) 4*2 
YES  CHANGE  Y-AXIS  LABEL' 


•PLOT  COL  1 

!  *  IS  PLOT  CHARACTER  IDENTIFIER  INOT  CENTERED'! 
!  PI. TIT  COL  7 

•PLOT  COL  1 


4  'PLOT  COL  4 
X 

4,7/  'CALL  PLOT  BRANCH  FOR  NEXT  PLOT 

1,1/  TF  PLOT  IN  THIS  FRAMf;  STARTING  AT  LFVfl  l 

NO  no  NOT  CHANGE  PLOT  LABFL 


0,100*0, 3000 

Y7S  OWIWGF  X-AXIS  LABEL 

RMS  OTSPLACFMENTS,  OB 

NO  OTT  WOT  CHANGF  Y-AXIS  LABEL 

3  'PLOT  COL  3 


♦ 

4,2, 1*3  '/FRO  COLUMNS  1-3  OF  C  ARRAY 

4,3,0/  fRESET  SELECTED  VARIABLE  EARAWETEPS  ,  LEIWINC  OTHERS  AS  BEFORE 
,30,61*63,64,19/ 

/ 


4,1, 

4,4, 

4,4, 

4,12 

3/ 


7,')» 

4,7/ 

1,1/ 

NO 

NO 

34.8 


1,3  'READ  VARIABLES  INTO  COLUMNS  t-5? START  AT  TOP  OF  JSHP.PTN  LIS 
1,5  IOTVIOF  COLUMNS  T-3  BY  NUMBER  OF  STATIONS  -  .  -j 

1*4  'WRITF  C  ARRAY  TO  PRINT. PPT 

/  TRETURN  CONTROL  TO  PTPL3 - 

'SET  TSSW  ARRAY 


:  •,>1  (.numv,  it"  j  i  <i  rirn  i  iri  riui  j 

fCALL  PLnT  BRANCH  FOR  NEXT  FRAME?  SFE  EARtIFR  DESCRIPTION 


,36.8,0,3000 


SALINITY,  PPT 
NO 
2 

4,7/ 


WIU'MBHVIU 


11600  2,22,0-1000 
II TOO  YF 

11S00  POTENTIAL  TEMPER  ATURE *  OEG  C 
VI 900  NO 
12000  1 


100  4, 

12200  1*1 
T7100  NO 
12600  NO 
17500  0,5,0,1000 

12600  YE 


;  12800 
17000 

i  nooo 
i;  moo 

S  11200 


11400 
11500 
11600 
1 1700 
11800 


14000 

14100 

14200 

14100 

14400 


40,220,0,1000 

YF 

DELTA-E,  IE-6  CM**3/GN 
W 


.iiaiuvi. 


! NOW  GOING  TO  PLOT  5  <  THETA  I  RATHER  TURN  SIP),  CHANGE  Al,A? 


14600  NO 
1 4  TOO  NO 

14800  34.R -16.8,2,22 

14000  YF 

16000  SALINITY,  PPT 
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18200  POTENTIAL  TEMPERATURE,  DEG  C 

15100  2  - 

18400  4,7/  ’CALL  PLOT  8®ANCW  FOR  FINAL  FRAME?  DEEP  TNETA-S 

T550TJ  1,17  -  - - - - 

18600  NO 


16  700 

18800  14.01,18.01,2.6,4.8 

15900  NO  -  - 

16000  NO 

T6100  2  - 

16200  4,12/  ‘RETURN  CONTROL  TO  PFPLS 


16400 

16500 


SPRINT/OFL  PRrNT.PTN.PTrTTrr.PFT 


4 


6 


18 


20 


22 


8  10  12  14  16 

POTENTIAL  TEMPERATURE,  OEG  C 


Figure  8:  Example  plot  from  ENER6Y.COM 
Potential  temperature  in  *C  averaged  along  adiabatical ly 
leveled  surfaces. 


.5 


1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5  5.0 

N.  CPH 

Figure  9:  Example  plot  from  ENERGY.COM 
Buoyancy  frequency  N  in  cph  averaged  along  adiabatically 
leveled  surfaces. 
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1 


LDE  FS5  TEST  0F  STRNDRRD  VERSI0N 


PATENT I AL  TEMPERATURE,  DEG  C 


LDE  FS5  TEST  0F  STANDARD  VERSI0N 


Figure  12:  Example  plot  from  ENERGY.COM 
Deep  potential  temperature  vs  salinity. 
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50  !•  POTFN.Cn*1!  #♦♦*♦♦♦*♦♦♦♦♦♦ 

100  ItTASK  IS  TO  SFT  UP  P ARARETfRS  FOR  AN  INTERACTIVE  RUN  OF  POTEN 
110  S • F INAL  INSTRUCTION  IN  THIS  CON  FILE  CHANGES  KIN  TO  ft*  THEREBY 
120  SIPETORNING  CONTROL  TO  THE  TERMINAL. 

130  t! 

200  SBSSTGN  PRINT. PTN  EOROOA 

“3D6  tASSIGN  JSHPFS5, PTN  POROl?  I  CHANGE  JS  HP  E  S  5  TO  CORRECT  STATION  LIST 

500  RUN/NOOFP  POTEN 

“TOO  YYT~TNTTmT7F  “TTONOOK - 

700  NO  00  NOT  INITIALISE  REGRESS  ION  PARAMETERS 
716  2.5  TPRFSSUPF  I NTFRVAL  FOR  INPUT  DATA 

1200  0  ««.*■«•  6/  !  CONTROL  NtlW  RETURNS  TO  TERMINAL  FOR  INTERACTIVE  SESSION 
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5  l»  *************  PEPLT.COM  *♦**♦♦♦♦**♦♦♦* 

10  4  !  T  A  9K  IS  TO  SFT  UP  P5RAMFTERS  FOP  AN  INTERACTIVE  PUN  0E  PFPLT 


iunnmun  n  iii  ■■  utiii  ■inmni  i  ivm  n  ■»  ■tni  n  n  i 


tv  *  .  un  j  i  jt  ^  i  l  n  l  '»  1  t  n  i  n  i  J  t,  u  n  rite  t  n  unvju  j  i  n  i  u  w  9  1  ntwen  1 

30  4'RETURNING  CONTROL  TO  THE  TERMINAL. 

ATI  1! 

100  4ASSIGN  JSHPF5.PTN  F0P012  ! CHANGE  JSHPF5  TO  CORRECT  STATION  LIST 

“125  VASSTGN  PIT1T.PP7  TOROOP  ! P1T0T.TPT' T5 “M E T I CUITE  FILE 

ISO  4ASS1CN  PRINT. PPT  FQR00A  IPRINT.PP'T  IS  LIST  FILE 


11IUUA  J 


300  NO  DO  NOT  READ  IN  PREVIOUSLY  STORED  COMMON 
500  YES  INITIALIZE  DATA  SELECTION  PARAMETERS 

600  4,?, 1,6, ..6/  'CONTROL  NOW  RETURNS  TO  TERM! NIL  FOR  INTERACTIVE  SESSION 


69 


50  S!  *•*»♦**♦♦*♦♦♦  TABLE.COM  ♦*♦♦♦♦♦♦♦♦♦♦♦♦ 

100  $ ! T  A  SK  IS  TO  CRFATf  TWO  PLOTS  CORRESPONDING  TO  FIGURES 

— 7uu — rr*wr  5  vr  wmrm-r  blue'  tovfv  tfeptirt;  thfttrst 

300  \'F1GU»F  IS  BUOYANCY  FR’EQUFNC  Y  NASA  FUNCTION  OF  DEPTH 
~  500  i'ANB  POSITION.  THE  SECOND  IS  VERTICAL  0! SPL A<TEWE NT. 

500  S'POTH  ARE  CRFATEO  USING  TABLE  SUBROUTINE  OF  PFPLT; 

ROD  TITHE  FIRST  FIGURE  USES  THE  OEEAULT  PLOT  SPFTT  FT  CAT  IONS . 
700  $  A S S I ON  PLOT.PPT  FOR008 

- BUD - STVSTUN'JSHPFH.PTN  FUR01? - 

900  RUN/NODFB  P FPLT 

1000  NO  00  NOT  »EAO  IN  PREVTUOSt?  STORED  COMMON 
1100  YFS  TNTTIALI7F  DATA  SELECTION  PARAMETERS 
T700  1,?/  ITNITIALI7F  PLOT  P1RWETERSTN  TABLE 


1300  16/  IFNTFR  NUMBER  OF  STATIONS  IN  PLOT 

— 1750 — 37 - rrnrm  pst  xrm- potttto - 

1 AOO  YES  INPUT  NEW  PLOT  LABFL 

1500  LOF  FS5:  n,CPH 

1600  NO  no  NOT  USE  DEFAULT  AXIS  PARAMETERS 

1700  /  !U5E  THESE  MTN  A NO  MAX 

l BOO  YES  CHANGE  X-AXIS  LABEL 

- T«rOO - RW  FAST  OF . 0RTT7TN - 

2000  YES  CHANGE  Y-AXIS  LABFL 

7100  K*  NORTH  OF  ORIGIN 
7700  1,0/  *PLDT 

7300  1,1/  TFHANGE  PrRAWFTERS  FDR~  SECOND  PLOT 

2 AOO  16, ,,100,3000/  1100  IS  DISPLACEMENT  AXIS  IN  OB 

- - jy - 

7500  YES  CHANGE  PLOT  LABEL 

7600  LOF  F75:  VFPTICAL  DISPLACEMENTS,  OB 

7633  NO  00  NOT  USE  OEFAULT  AXIS  PARAMETERS 

7666  /  !USE  THESE  MIN,HAX  VALUES 

2700  NO  DO  NOT  CHANGE  X-AX35  LABFL 

- TirtTO - NIT  00 'NOT  CHANGE  T=TTT*S  C1TBFC - 

7B16  5,3,0/  'CHANGE  PLOT  PARAMETERS — AVRGS 

ZB37  / 


7858  / 

7P65  ,.05/  IRFSCALF  DISPLACEMENTS 

2880  ,50/  ! PLOT  DISPLACEMENTS  RATHER  THAN  N 

7W9TJ - T7T77 IRTTUR  N  TO"  FEPLS - 

7900  1,0/  'PLOT 

300U  7/  rEXTT  PROGRAM  .  -  . 

3100  YFS 


tf  ♦♦##*#**♦*****  dynht. com  *♦**♦•#♦•****♦ 

t ! TASK  IS  TO  CREATE  GPCP  COMPATIBLE  OUTPUT  FRO»  R.AVG  FILES  AT 


tiiiiuiii 


TOO 

WOO 

<>00 

IOOO 

1100 


1250 
1 100 
1800 
3  wo 
3500 


$!XLONC,OYN  HT  REF  TO  PF  AT  LEVEL  507T0,S0,REE  SPECIFIC  VOL.'TIWE 
srrjUCTfN  DAYS  FRTJM  1  JANWFCTMAL  HOTJRSTJFTJRMITTS 
t'lH  VAS,IA,2Ff0.?),3F8.3»  28  FAY  81  NAN  BRAY 

♦IWUWRFR  OF  STATIONS  OUTPUT,  OUTPUT  VARTITLFSY  AND  REFERENCE  PRESSURE 
t ! FOR  DYNAMIC  HEIGHT  FAY  BF  CHANGFO  IN  AVRGS  BRANCH  3. 


l  ;  u  i  i  u  iil  i  un  i  uui  ryi  i  i  u  inpnn,  LCll  inc  (CR}| 

t! 

♦ASSIGN  JSHPFS.PTN  FOR012  ! CHANGE  JS1^5~T0  APPROPRIATE  STATION  LIST 
RUN/NWOEB  PEPLT 

Nn  no  MOT  READ  IN  INITTALT7ED  COMMON 
YF  INITIALIZE  OATA  SELECTION  PARAMETERS 


Y  f  T  B  If Of  ;  tCBU 

8*3,17  ! SFT  OATA  VARIABLES 

1000718,19,13, 18, 10, 7, SO, 55/ 

/  _ 

/ 


YES  INPUT  NEW  FILE  NAME 
TEST. 'DAT  " 

IS* l 7  (RANGE  OF  LEVELS  FOR  WHICH  DH  WILL  BE  OUTPUT 

8,177  TRFTWN  CONTROL  TO  REELS - “ - 

7/  ’EXIT  PROGRAM? 


3700  SSORT/KFY ■(POSITION:?, "S I ZFsftl  TEST. OAT  TESTP.URT  (SORTS  BY  PRESSURE 
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50  POTEN*  rrfNSs  SHORT  DOCUMENTATION 


100  KB#  ISW  JSW  KLIST  DESCRIPTION 

~T50 - 0  -  - - - - SHORT  50Cu«ENT*tt0N - 

200  1  «  -  -  COMPUTE  REGRESSIONS  E-Ot  ISW 

100  SEQUENTIAL  ST  ATT  TONS. 

A 00  2  INITSALIZF  OUT*  SELECTTON 

100  PARAMETERS. 

*00  3  0  SET  PARAHFTFRS*  SHORT  LTST. 

TOO  I  - - = - S  ET-  FIMIWETFRYT  FULU  TTSTT 

SOO  STORES  COMMON  TN  FILE  RPTC* 

000  4  l  -  LU  AVRCP  AVERAGING'  SUBROUTT NF. 

1000  5  _NOT^JSeD 

1100  6  -  -  LIST  LABEL 

1300  7  NOT  USED 

rtrm  *  - - - - - - 1  i\i  wcrrrny  t  tub  matt  wb - 


1500  B 
1600  10 
2000  11 
17 


LIST  DATA  RECORD1. 
NDT  USED 

SET  tSSW  (SWITCH* 
T5TTT  PROGRAM. 


ARRAY. 


73 

93  c  CONPOTEN.EOR  FILEi  01MENSI ON, COMMON  AND  EQUIVALENCE  FOR  POTFN 
96  C  AVAILABLE  POTENTIAL  ENERGY  PROGRAMS.  N.8RAY 

- ^ - c - 

100  PARAMETER  KC*»-  235 

700  BYTE  LBtPROVFR 

400  C 

433  C  01  MANSION 
466  C 

“TUU - DIMENSION  KH0CT1WJ  ,KBUF  ( 46 1  ?VR IT) - 

600  DIMENSION  KPTCNIKCM) 

650  DIMENSION  PRESS  I  33001 


700  C 

716  C  BLINK  COMMON  ' 

737  C 

750 - COWNTTNTTTN' - 

775  C  BEGINNING  OF  STORED  COMMON 

500  COMMON  KTTX,KLtST, KOUT , KTP  , TSW ,JSH  , K ft R 


525  C 

550  C  BEGINNING  OF  HEADER 


875  C 

too - nTMMowrTipnE  ;5rnrwn  con  ,  t  SHFixe astvda  — 

1000  COMMON  XLAT,XLQNO,WOT,XLTO,XLGO 

UOO  COMMON  L  B  8  L  f  3 )  ».L'BL'f  1 31, NSC (  60TVNRR  (60)  «P#SE C t  ION 

1125  C 

TI50  C  8FGTNNTNG  OF  OATA  BUFFER  X OOF  . 

1175  C 

7200 - common  rrypEvMmjF-,  iRFCtw,Npp^xsiifTi»L2 - 

1300  COMMON  PF  ,  TO,  SO,‘OVO 

1X00  COM-ON  P  T  , THF»SF**OVF  . . 

1500  COMMON  PM,  TMM,SW,OVM 

TOOTT  COMMON  DH,TFfTFF 


1700  COMMON  CP<8),7l,CTf 5),Z2,F1*F2,E3 

T7T7 - C - 

1724  C  FND  OF  KBUF 

1716  C  -  ' 


1750 

1500 

1550 

COMMON 

COMMON 

COMMON 

DELP , DP 

A1 , A7* A3,N1 ,W7VN1  ' 
C  16  )  ,  ISSNM6) 

1556 

ICRUI  St  IMffJi^tOVER 

1959 

rOH"ON 

J*AX 

1562 

C 

—  — 

1874  C  FNO  OF  STORED  COMMON 
1 556  C 

2000  COMMON  P I  3300 1  * T( 3300) *S  1 3300) ,0V f  33001 

7TU0 - COM  MDN  T H I  3  3 00  )  ,  P TT7 3 0  0 )  ,  TT  n 100 ) - 

2200  COMMON  Rm,BP<8),BTm,BA(8),C0<36),MRt8) 

2100  COMMON  EYI6), JEYI6)  -  — 

2400  COMMON  WTI600), JSHP1600) 

2475  COMMON  DATlYf 3300,01?)  ~ 

2450  COMMON  JSTN, JRMAX,M1,M2 

7475 - COMMON  tLRgC»«KTT - 

2500  C 

2533  C  EOirPTALENCE  - 


2566  C 

7600 

2700 


FOUTVALENCF  I KF|!TO,l.TYPTT, TKBUF iJRTYPF) 
EQUIVALENCE  <P01PF,A1) ,!VR,PF)/tKTTX,XPTCM) 
TOUTVACTNCE  I PRESS  ,OAT  AX) - 


2800  C 


1  100 
200 

c 

c 
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INDEX  RECORD  FIELDS  DECOMPOSED 

USES  LABELLED  COMNOM  X  USE  INCLUDE  STATEMENT  TO  MERGE  INTO  PRODS. 

TOO” 

c 

R0M  MAR  27  1980  ~~ 

500 

c 

ARRAY  IOXREC  CONTAINS  THE  INDEX  RECORD 

500 

c 

ALWAYS  THE  FIRST  RECORD  OF  A  DATA  FILE 

-  * - 

■4  600 

C 

IOXREC  HAS  THE  SANE  STRUCTURE  IN  THE  SUBINDEX  FILE 

700 

800 

C 

COMMON/TNOX/  I0XRECI256) 

• 

1000 

1100 

1200 

TWO 

1*00 


1500 

1600 


1800 
“  "I960 
2000 


2200 

2300 

2*00 

2500 

2600 


2800 

TWO 

3000 

TOO 

3200 


J  J  vu 

3*00 
“  1500 
3600 
~  7700 
3800 


*000 

*100 

*200 

5700 

**00 


C  LPGVER  IS  AN  ACSII  DESCRIPTION  OF  LOAOING  PROGRAM  VERSION 

INTEGER  LPGVFRf 6),LSTREC 

C ♦♦*♦♦•♦**♦♦♦*♦♦ *********************************** ************** 

C  IMW  AND  LTtNE  CONTAIN  ASCII  OATEfT f MF~FTLC'  CREATED 
INTEGER  L0ATE(3)VLTTHE(2I 


VICE  AND  FILSPEC  FORM  A  COMPLETE  FILE  SPECIFIER 
INTEGFR  FILSPEC(8)*0EVICE 

C  IFWEO  ARRAY  HAS  TWF  SAME  STRUCTURE  AS  CT0T8  VAX  FORMAT 
INTEGER  IFHEOCW)  4ICMNTI35) 

C  CNTRL  HAS  THE  LOCATIONS  OF  Tiff  BEGINNING  WORD  OF 
C  INFORMATION  FIELDS  0*  THE  INDEX  RECORD 


*600 

*700 

*800 

*000 

5000 


'  IUV 

5200 

5700 

5*00 

5500 

5600 


5800 

5500 

6000 


INTEGER  CNTRL ( 6  )■ 

VARDES  ARRAYS  CONTAIN  A8RE V I ATEO  VAR TABLE  0E5CR I PTORS 

MIN/MAX  VALUES  IVARDES  CONTAINS  MNEMONIC  IDENTIFIERS 
Of  MENS  ION  VAROES*  *  1 16T*  f  VARDFST  4*T6)~ 

EQUIVALENCE  ( KEYWO* I DXREC 1 1)1 


LSVRFC  IS  THE  NEXT  AVAILABLE  RECORD  IN  SUB INDEX  FILE  1ST  RFC.  ONLY 

f out valence i lstrec* idxrec  t9i  i  -  - 

EQUIVALENCE! TFNEOI II, IOXREC (13)1 

FOU I VAL ENC E  ( VARDFS < 1*0,  fDXRF C f  1T5TT,  IfVARDF $11 4 llVlDXRCC  111 5 1 1 
FQUI VALENCE ( KSCAN, IOXREC (105)1 


EOU 1  VALENCE  ( NSCA^*  IOXREC ( 106 11VINSCANS **IOXR EC  11121 1 
TOUT  VALENCE  *  PMTN/TOYRE  CITIOT V, *  P WS INT,TtTTR  E  C (111  IT 
EQUIVALENCE INTOT, IOXgEC< 108 1 1 

EQUIVALENCE 


2,'!IFHED(7),TYR),!  IEHEO I  8 1 ,IMO> ,1 rEHED! 5) ,’TDA 1 

3^fIFHED*ITltILTS0I,*TEHE0TT71,1ET5in - 

4/(IFHE0(l3)*ILNS0>,!IFHED(l5),tLNSM> 

5*1  IFWED  ri5)  *TN1P»77  * TTETIEOIIST  *7WRlt  Ml  FWEUTT7 
6y< IFHE0(19I* ILTED),(IFHED«20),rLTEM) 


If il  aananrj  in  iiaatifia  ianaii  r/j  n  i  ■  •  i  aj  ■ 


8V!IF HE 0(231*  TETNE) *<IFHED(10)*ISTNEI 

94  f  IFHED*  38  KICA’ST),  f  TEHFD  f 277  ,JDATT,  MPNFDf 281  *TW3  T) 

X*(IF HE 0*55), ICMNT) 

C  *♦  «  *  *  *#  *♦  *  *  •  *  **  *  *  *  *  ******  **  **  *  *  *  *  vmmTTVr*  C**  *  MM  M*  M  ****** 
EQUIVALENCE 


1701  »  «  iriLirruinn^^eu 

2  *  (LDATF flit  TOXR’ECI  195  1 1  *  (LUMP!  II*  IOXREC*  198 1 ) 

7 ** LPGVERTT1  *1  DXREC  1 99  Vf~ - - 

C  RMAX  IS  THE  LAST  RECORD  OF  A  DATA  FILE 

C  TUXT7TJC  IS  TWE~  REEURfO  IT  OF  THE  INDEX  FILE  «  7  FOR' SINGLE  ST  ATI  0? 
*, (RMAX* IOXREC ( 1091 1  *( I DXLOC* IOXREC! 20811 


C *************** *****************  ************************************ 

C  END  UABELLEO  CTJWMON  FTJR  TNPEY  RECORD -  -  - 
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'  1  I0O  C  PDTFN  PROG  ♦♦♦♦*♦♦♦♦♦*♦*¥♦♦***♦♦**♦*♦*♦♦**♦*♦♦♦♦♦*♦*♦*♦ 

200  PROGRAM  POTEN 

- TUO - c 

400  c 

"  - 500  C  PROGRAM  TO  COMPUTE  RFFtftEWtE  SOffFAOFS  SPLAT fVF  TO  PF  FO* 


600  C  CALCULATION  OF  AVAILABLE  POTENTIAL  ENERGY,  REGRESSION 

- TOO  ~  C'FTTY  ARF  "AOF  TO  PTTFSTORF  ANT  P'OTETETTAL  TEMPERATURE  AT  ' 

BOO  C  ^UNCTIONS  OF  SPECIFIC  VOLUME  ANOMALY  AT  PF* 


- "TOO — 

~r 

' 

1000 

c 

JUNE 

28  1076  N  FOFONOFF 

TTOO 

c 

1200 

TTOO 

c 

TNCLUOE  *  CON POTEN* FOR* 

1400 

c 

OPEN 

BINARY  FILE  FOR  STORAGE  OF  COMMON 

- TTOO — 

1600 

c 

OPEN! UNI T ■ 10  »NAMg«*  K  PTCM. OAT *V ACCESS** 01 RECT**TYPF«*OLO** 

~  ITOO 

V 

RFCOROTYPP-’ FfYFO'*  *PEC0PDS77FSHCM*ERR*TT00) 

1800 

ITOO 

c 

l 

continue 

-- 

2000 

KIN  -  5 

_ J 

- TTOO — 

TCTTST  -~5 

2200 

KLIST  -  4 

7700 

KOUT  •  9 

2400 

KTP  -  l 

- 7500 

30 

MR ITE  f KTTX» lOOOJ 

2600 

1000 

F0RMATC1H  * *POTENt  POTENTIAL  ENERGY  PROGRAM*! 

^  tv u 
2800 
7000 
TOGO 

non 

3200 
3300  " 
3400 
7500 


v,  all.  ntm 
CO  TO  50 
C 

C  CREATE  NEW  BINARY  FrtF  FOR  STORAGE  OF  COMMON  IF  NO  OLD  ONE  EXISTS 

c  '  .  . . . 

1100  OPEN(  UNI  T-10 .NAME- »KPTCM. OAT  »  VACC FSS «** OTR ECT*  ,TYPF«*NEN** 

- * — RECnROTYPF.'FnrFD'tRECOinJSTIF-KCMtERR-mTOI 


GO  TO  1 
50  FNO 


tl.M 


400 

500 

600 

TOO 

800 


woo 

1000 

1100 

1200 

TTOO 

1400 


€  PTE  NS  SUBPROG  POTEN  PTENS.FOR  PILE  ♦♦♦*♦ 

"SUBROUTINE  PTENS 


riTtTTTITTITTII TTTTI TTTTIIITT TTTTTT ▼ TTTT TT T T TT T Til If I1IIIT] 


c  PROGRAM  to  CONFUTE  REEtERENCE  SURFACES  RELATIVE  TO  PF  FOR 
C  CALCULATION  OF  AVAILABLE  POTENTIAL  ENERGY.  REGRESSION 
C  rm'MF  MADE  TO  PRESSURE  ANOPOTENTTAL  TEMPERATURE  AS 
C  FUNCTIONS  OF  SPECIFIC  VOLUMF  ANOMALY  AT  PF.  _ _ 


C 

C  JUNE  28  1976  N  FOFONOFT 
C 

c  MODIFIED  TO  ACCEPT  CT078  VAX  DISC  OATA  AS  INPUT  15DEC80  N.8RAY. 

c .  . 

OINFNSION  Df5),DOCI10) 


i  - T500 

1600 
:  TTOO 

1800 

C 

-■  - 

INCLUDE  * COMP OTEN. FOR* 

CHAR ACTER*8  OOC 

- .  . 

a  T9O0- 

2000 

C 

KIN  -  5 

"T  - TTOO - 

KTTT  *  6 

:  2160 

WR ITE 1 KTTX ,40) 

7270 

40 

F0PMATI1H  ,*INTTnrtT7F  COMMON  (YES  OR  NO)  V  I 

1  2280 

IF  I  NOYES  f  KTN,KTTX ) «NE • 1 ) GO  TO  14 

- 7740 

C 

.  . 

2400 

C 

INITIALIZE  OATA  SELECTION  PARAMETERS 

TT 

2520 

CALL  DATA(KTPy-l) 

r  7500 

GO  TO  30 

. 

2640 

14 

RFADIIOUIKPTCN 

7700 

18 

WRTTE (KTTX, 2 01 

2800 

20 

FORMAT ( 1H  , 'INITIALIZE  REGRESSION  PARAMETERS 

ITES  OR  NO)?*) 

,  - 2500 - 

. 

3000 

RFADI10* 1,END*10) KPT  CM 

•  ‘  3700 

5 

WRITE  IXTTT4  25) 

"  i  7300 

25 

FORMAT! 'OWHAT  TS  THF  RESOLUTION  OF  THE  TNPUT 

DATA,'  IN  OB? *1 

3400 

RFADIRTN, *)OEUR  - 

3450 

1  -  H  ■■  —  -■  ■  ■ 

KLIST  -  6 

- - - - - 

1 


tl*l 


1  3600 

4  ~  7WU0 
*  3900 

“4000 
M  4100 


HUM  1  *  ,  1  Rt  U  I  r» 

1005  FORMAT  ( IH  ,  •  POTENt  KB  R  •  I SW,  JS  W*RL  I  ST*  MOOT  ,  KTP  yK  IN  *  *7,  714  ) 

R E AU IlCIN » *  1 KUR" t’T SIT*  J3W*  K L  I  S  »  yRUU)  » X  rPVKTN 
IFfKBR.GT. 12IKBR-13 

tf  cRBirnTOOTiiooviT  —  - ~ -  “ 

IT  GO  T0( 100 *200 *300, 400, 500, 600 *703, 800* 900* 1000,1100,1200, 100* 


■ 

H 


4200 
4  TOO 
4400 
4500 
4600 


4800 

“4000 

5000 

TTOO 

5200 


5400 

5500 

5600 

TTOO 


C  ♦»»>»♦»♦»»» *TNIT I ALT7ATION 
15  KTYPE  •  0 
MHOR  •  150 
MPUF  •  46 


:  wf  an  ■  fti  I.  ri] 


NPR III  «  4 
NPRT7T  *  TT 
NPR  f  3 )  -  17 
NPR  1 41  -»•  24 
NPRC5)  •  50 


NPR (71  «  100 
NPPT8T  •  700 
NPR (9 1  •  500 
NPR  HOT  *  500” 


5900 

— BTTTTO - 

5100 
"  6?50 
6900 
6900 
6500 

~rms 

6  700 
6000 
6900 
7000 
7100 

- 77TT0 - 

7300 

7400 

7500 

7600 

7700 

— 77T0 - 

7720 

7730 

7740 

7750 

7800 

— T955 - 

8000 
ITl  00 
8200 
8705 
8400 

— *500 - 

8600 

8700 

8800 

8900 

8925 

— 8950 - 

8975 

9000 

9100 

9700 

9300 

9400 - 

9500 

9600 

9700 

9800 

9900 

toooo — 
10100 
10200 
10300 
10405 
10500 

10600 - 

10700 

10800 

10900 

11000 


NPUdll  -o 
NPR  f 1 2  >  ■  2 

- upvamr-T - - 

NPR (14)  •  13 
NPRI15)  -  13 
Nsctn  *  o 

NSrm  .  400  - 

NSC (31  •  1000 

wsc(4»  •  mo? 

NSC  ( 5 )  «  3500 
NSC (6 )  «  6 
NSCC7)  -  5 
NSC<8!  •  4 
NSC (9)  -  3 

- Nsriim  -I - 

NSCIll)  -  20 
NSC ( 1 2 1  «  30 
NSC  C 1 3 1  -  40 
NSC (14)  -  50 
NSC (15)  -  *0 

pr$rf\rpir~-  "nr* - 

ISHP  -  'GY* 

rc»urs  -  i 

IPPOJ  -  3 
GO  TO  5 
30  DFLP-2. 

16  VR  C  J  >  -  0.0 

POIFF  -  6.0  - - 

A 2  •  3.0 

43  -  3.0  - 

LTYPF  -  1 

- TTW«"|j - 

N  -  2 

NOP  -  TO  .  . 

KSW  -  1 

WGT  «  1.0  . 

00  17  J-1,16 

- nsTctiinr - 

17  CONTINUE 

GO  TO  18 

C  *♦♦*♦♦♦**♦♦♦  SELECT  DATA  AND  COMPUTE  «1  ♦*•♦♦♦♦♦♦*•**♦♦* 

too  call  corps'  '  - - -  -  -  --- 

GO  TO  10 

T - 

C  IN ITTAL  T  7E  OATA  SELECTION  PARAMETERS  92  ♦♦*♦**♦♦♦*♦*♦♦♦♦♦♦ 

C  •  •  - . . - 

200  CALL  DAT  A ( K  TP  V-l 1 

C  p»p******»pp5E'r  PARAMETERS  «3  r*r*4*4»»****r*9*r*»***r »♦ 

300  MRITFfKLIST,30001!C0NtKSW,A2,A3,HGT,P0fFF 

- read  » i  fcmn  r  s  wy  a?t«3  ?  wgt/pdi  pf - 

W»ITEfKLIST,302010FLP 

RFAO(KTNy*10FUr 

IF<ISW)10,10,>I5 

c 

C  SUBROUTINE  TO  ACCEPT  REGRESSION  PAR ANETERS  IN  ENGLISH  ANO 


315  CALL  PARAH 

WRTTE<*LIST,320» 

C  PRTNT  OUT  POTFN  FORMAT  PHriPETEHr 


11100  320  FORMAT* *OSFCTTON  LEVEL  lfvel  index  start  «  OF  «  OF*,/ 

11200  X  1H  ,  *  NUMBER  NUHBFR  INTFRVAL  PRESSURE  TERNS  CYCLFS*) 

75  nsE-nsfction 

11400  00  330  1*1* NSF 

m60  11*1  ♦  NSF 

11600  12*1 ♦PRNSE 

TI70O  M# !TE ( KL 1ST  , 335 ) T, NPR I  I ) , NPR I f 1 1 , NPR I 1 2 ) *NSCl  T )VNSC« II ) ,NSC<I 2  > 

11800  335  F0RNATU3*6I  *> 


12000 

12100 

12200 
T2300 
12400 
~00 
12600 
12700 
12800 
12000 
13000 


13200 

13250 

13300 

13350 

13400 


C  CHANGE  OR  LIST  DATA  SPLCCTION  PARAMETERS 
CALL  0  A  T  A  I X  T  P  ,0  ) 

C 

C  CHANGE  OR  LIST  DATA  LAB Ft  I  IDENTIFIES  THE  DATA  SOURCE  1 
C 

MR ItEfKTTX* 30151 ILBC II), I *1,13) 

IF  I  NOYES  I K IN, XTTX) ,EQ • 1 ) THEN 
MRtTEIKLIST, 3-010) 

RFADf KIN, 6005)1 LBLIX ),K«1, 13) 

FNOIF  ~  ' 


niLi: 


IMJ 


14200 

14300 

14400 

T4SUU 

14600 


fUU 

14800 

14000 

15000 

T5T00 

15200 


IF  (  JSW. NE  .  2  )  THEN 
WRITE! 10* 1)KPTCN 
FNOIF 
GO  TO  10 


FORMATUH 
FORMATUH 
FORMAT  1 1H 
FORMATUH 

format  I  ih 


,  '  ICON  ,.K5W,S0P»SDT  «/WGT  »P0IPF**/»2I4*3F6.?,F7.0) 
,* NSC SP,N, NOP*,/, 51 5,1013) 

, *  PRESSURE  CONSTANTS*,/, I5f 4) 

,* INSERT  LABEL  <27  CHAR.*) 

,  *  INPUT  NEW  LABEL?  OLD  LABEL  IS:  %/,?H  ,13A4) 


C  ♦♦♦♦♦♦♦♦♦♦♦♦AVER AGING  SOBROUTTNF  r4"  **#*4*****#*4»4.*4*4*4#*44 
400  CALL  AVRCP 
GO  TO  10 

c  ♦♦♦♦•♦♦♦♦♦♦♦  NOt  PRESENTLY  USED  ♦♦♦♦♦♦♦♦♦♦*♦♦♦ 


C 

C 

C 

C 

C  501 


Mil 


NCR  -  5 
XTNP  -  5 
JM AX  -  23 
ROUT  -  I 
no  505  M-9,13 


15400  C  505  CR IX , M)  -  0.0 

V55U0  C  507  MRTTFIXTTX,5D1D)X0UT,NGR, JNRTfXTNP 

15600  C5010  FORMATUH  ,  »  AVDVF  t XOUT ,NGR*  JNAX,K INP*  ,7*4 1 4 ) 

T5700  C  512  RFADrXTN,#)XWr^Cr,  JETAr^nwr 

15800  C  DO  530  J*1,NGR 


16000  C  REAO(KINP,*H,NST,XPR,  rDIK),K-I,5) 

I6TOO  C  1IFRIT1  -  XFX 

16200  C  no  520  N-9,13 

T67DO  C  520  CRf T ,NI  -  CPU, Nl  ♦  OfN-B? 

16400  C  530  CONTINUE 


16600  C  00  540  I«1,JMAX 

T6TO0  C  540  CRII,J)  •  CR  II  ,  31  /FLOAT!  NCR  1 

16800  C  545  DO  550  I-1,JMAX 

16000  C  MR  TTF  TXTJUT  ,  5  BOOT  I", "NS  T  ,  I IPR  mVTCRTT,  XT*X*9 , 137 


79 


17000 

17100 
r 72U0 
17300 
17*00 
17800 

17300 

17900 


18100 

18700 

18300 

18*00 

18800 


18700 

18000 

18900 

19000 

19100 

"T970U" 

19300 

19*00 

19800 

T9600 

19800 


70100 


c 

C 

XT 


880  CONTINUE 

_ KOUT  -  t  ]  bt>VKUUT 


C  GO  TO  (801,80?V8*8,10)I$W 

C8000  rO^HAT ( 2 1  *  »  T8V8P10. *1 

C  ♦♦♦♦♦♦♦♦♦♦♦♦LIST  DATA  LARFL  *6  ♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦*♦♦**♦♦♦♦*♦**♦♦♦< 

800  HR  ITE  f  KL  I  ST* 6008  MLBL  (  K )  *K»1  *  13  ) 

GO  TO  10 

~8W8 — TTorir/rrmr 

************  y 7 

703  MFOF  -  0 

707  CALL  RFAO(KQUT,KBUF»  MBUF, IFOE 1 

IF(IF0F)720,7O*,70* 

70*  MFOF  -  0 


Tram - 

NOT  PRESENTLY 


USED  ************************ 


-Ttrs — rrryTYPFT90tn-90<n7io 


710 

718 


770 

778 


00  718  ".i,khO* 

KHDG(M)  «  KRUFfM) 

GO  TO  800 

TF I MF  OF ) 1 0  ,775*10  - 

MFOF  -  1 

- T0”10 - 

703  GO  TO  10 

♦♦♦♦♦♦♦♦♦♦♦♦LIST  HEAOE*  RECORD  *8  *********************** 

800  WPITFfKLIST,  80001  I L« 8L f K I , K- 1 , 3  I , XL A T , XLONG, XLTO, XLGO 

HR  TTE I  KL  1ST,  8008  JLTY PE  *WHDRTTCUN,  1  SHE  *3X015  T ,  JOTIT ,TFR  ,LTR 
GO  TO  10 


8000 — rotmrTt 3a*  i;*t-b  .  n - 

8008  FORMAT!/, ♦  Typ^  MOHR  ICON  SHIP  CAST  JOAV  IPR  LPR*,/,8I5 
C  ♦♦♦♦♦♦♦♦♦♦♦♦LIST  0A7A  RECORD  *9"  ************************* 

900  WRITF1KL I  ST, 90001 IRFC ,PF, TO, SO,DVO, PMVTHM, SM*DVH, 71 ,72 

903  IFTISSWC8))  905, 'TO,  10 

908  WRITFIKLIST ,9008) !CP!K),K-l,N) 

JDrrTrnarrr,  wrim  ct  to  ,k  «i ,  n  > - 

GO  TO  10 

,  1 2 ,7 IF 7.  T, F7. 3 *F7 »3 ,F7 .71, P6.2,Fft. *  1 
,3HCP  »  6F  l 1  •  * ) 

, 3NCT  ,6FI1.*) 

TAPF  FUNCTIONS  »10  ♦♦♦♦♦♦♦•♦*♦♦♦♦♦♦♦♦♦♦♦• 
FTAPETTtWf JSK,KLTSTI - 


70700 

70300 

70*00 

70800 

70600 

70700 

70800 

70900 

71000 

71100 

71700 


71300 


71*00 

71800 

71600 

71700 

71800 

7T900 


77000 

77100 

77700 

77300 

77*00 

77300 


77600 

77700 

77800 

27900 

73000 

73300 

73700 

73300 

73*00 

73300 


) 


9000 

9008 

9010 


FORMAT! IH 
FORMAT  C 1H 
FORMAT! 1H 

C  •♦♦♦♦♦♦♦♦♦♦♦MAG 

croon 


crctr 


c  klist  -  ft 

1000  GO  TO  10  - 

C8  ♦♦♦♦♦♦  SFT  ISSN  SNITCHES  ill  ♦♦♦♦♦♦ 

1100  HR  ITEIKTTY,rT801  TKfK»l ,  X6TTT  TSTWIK  )  »K-T»X6T 

1180  F0RMAT(2(1H  ,7,161*,/),*  ENTER  K, l SSH! K ) •» ) 

- yy  ATH  K  I N  ,-RTT  KVTS  S  WEK  r,~M  ♦1,16) - 

GO  TO  10 


C 


C  ♦♦♦♦♦♦♦♦♦♦  EXIT  PROGRAM  *12  ♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

1700  HR  T  TF  !  K TTX , 1 71 0 ) 

rF!NOYFS(KIN,KTTX).NE.l)G0  TO  10 
- - - STOP- - 


1710  FORMAT  !  1H  ,  *EXIT  PROGRAM?*) 

C*»  POTEN!  SHORT  nOClJ9ET«rr*TTON— 8ITAWOT  0  ************** 


1300  0PFN(UNIT»80,NAME«*P0TFN.00C**TYPF  *  *010* , READONLY ) 


00  1380  N- 1,708  1 

RFAD<SO,132S,PND«1312)IOOC!I ) VI -1,8)  | 

1380 

1317 

1378 

1330 

l  '  HR  1  I  *  ,  l  HUM U'JvJ  1  1  1,1*1,81 

CONTINUE 

CLOSE ! UNIT-801 

FORMAT  1 8  A  8 ) 

FORMAT !  1H  ,8*3)  '  '  - 

— 

Ihhl 


COMPS  SURPROC  POTPn  ********************************4**** 
********************************************************* 


k-iM.l:hllk*ft.'l: 


lil-lii 


********************************************************** 


TO  COHPUTF  RFGRFSS  ION  CDFPFICIFNTS  AT  SPECIFIED  DEPTHS. 


JUNF  28  1976  N  FOFOMOfF 


1*1*1 


1000 

1100 

1200 

1700 

1*00 


1600 

T700 

1800 

IflOO 

2200 


(alt 


N.8RAY 


UiM.BiiLl.lll 


l»I' 


INCLUDE  'COMPOTEN.FOR* 


lllfl 


i*ni 


hUkiJiki 


NUMBER  OF  STATIONS 


100  CONTINUE 


INITIALIZE  AND  ACCESS  INOEX  FILF  A NO  CRUISE  INFORMATION 


2*00 
2500 
2600 
2  "TOO 
2800 


UJltlil 


3000 

77U0 

3200 

3300 

3*00 


CALL  0ATA<KTP,7» 


RFGIN  COMPUTATION  FOR  ISW  TOTAL  STATIONS 


IF  1  TSW.GT.LLRECHSN-LLREC 


IU' 


RF/TO  STATION  HFADER  FRO*  UNIT  XTP  UNO  CHECK  IF  IT  MFFTS 
SELECTION  CRITERIA 

READ  TFHPEPATURF  AND  SALINITY  DAT*  TNTO  DATA  ARRAY. 


*000 

*TOU 

*200 

*700 

**00 


GO  TO  200 
106  CONTINUF 

RFTUPN  TO  PTFN5 


*600 

*700 

*800 

*900 

8000 


COMPUTE  REGRESSION  VFR5T0N  OF  DATA 
IF  TSSWI131— 1  OUTPUT  TO  FILE  * . REG 

PTSCFLLANFOUS  ITDORDATTON  MAY  "BF  RFOUESTFD  TO  B"F  PRINTED 
TO  FILE  PRINT. PTN  IKLIST-*)  BY  SETTING  ISSW  VALUES. 


8200 

8300 

8*00 

8800 

8600 

"8700 

8800 

7900 

6000 

6100 

6200 


KTYPF  DISTINGUISHES  BFTRFEN  MEADE*  AND  DATA  RECORDS! 

0-OATA,  1-HEADER. 

KE.KT,K*»  APE  INDICES 
N  IS  POLYNOMIAL  ORDER 

NT7P  I~S~T  OF  DATA  CYCLES  TJ VER'WHICH  REGRESSION  IS  Pf» 
KERR  COUNTS  THE  *  OF  RFPL ACFMENTS  MADE  BY  SOSR  EDIT 
TN  EACH  RFGRFSS ION  INTERVAL 
IPR  ANO  LJP  KFFP  TRACK  OF  PRESSURE  AS  AN  INDEX 
IRFC  INDEXES  THP  LEVFLS  PF 


6*00 

6700 

6600 

6700 


KF  -  2 

N  -  NSCI NSFCTTON^l 1 

NOP  -  NSCI ?*NSECT ION* 1 1 
XNDP  •  NDP 


MOO 

7000 


KFRR  -  3 


82 


7i  t ro - ptfc  =~i - 

7400  C 

7800  C  CO*PuTF  REGRESSIONS' 

7600  C 

7673  C  JRAX  IS  THE  TOTAL  WT78HFR  of  lfvfct. 

7666  C 

-7EW2 - nrmT*r  wvrrrcwix - 

7714  P(«)*PRFSS(f*t 

7710  Tm-nATAXIN,.!  I 

7746  S(H)-DATAXf 8*2) 

7770  2188  CONTINUF 

7800  ?10  00  270  J-1,J8AX 

“WO — c - 

7900  C  SUBR  JPR  COMPUTES  CORRECT  PRESSURE  PF  GIVEN  SECTION  ANT) 
8000  C  INTERVAL  INF  OR  NATION 


8100 

8700 

8100 


10600 

roToo 

10733 
I0T66 
10800 
10900 
11000 
TT 070 
11040 
11060 
11080 
1 1 100 
13400 
17500 
13600 
11700 
13800 
13900 
14000 
14800 
14900 
15000 
18100 
~TT?oF 
18300 
1S4THJ 
18800 
15600 
18700 
~  1 8800- 
18900 
16000 


C 

IP  *  JPR(  TREC.NPHiNSECTrONl  “ 

If  I IP . FO. NSC IKF) I  THEN 

7157 - TT~g  NS  CT7F7  NYFTT  IT1W1 - 

NOP  -  NSC<KF»2*NSECTI0N> 

XNOP  -  NOP 
KE  -  KF  ♦  l 

FN01F  - 

INO  -  ( I P -PRESS! I 1 l/OELP  ♦  1 

- wt  „  tNO -NOP 77 - 

82  -  IN0»N0P/:2-l 

TFI8I.LT.11THE1T  - 

81-1  \ 

R2-N0P  '  . . .  . " . .  _  " 

ENDIF 

- IFI82.GI. JR8AXJG0  TO  280 - 

C  IF  PSSWI5I--1  WRITE  OUT  SCAN  i,  SCALED  PRESSURE*  TFNP, 
C  SALTNTTT.  -  - 

C 

TFITSSWI81 1 216*717,217 

216  WRTTF(KLlST,?160)(K,PtK)«T<KI*S<K),K**t,n2) 

7TFU — ^FnRNATirH  ,~1 4Y3F9.  31 - 

217  CONTINUE 
C 

C  PERFORM  REGRESSIONS  OVER  INTERVAL  CORRESPONDING  TO  PF 

C  FIRST,  FIND  WEANS  OF  P,S. 

C  ' 

220  PF  -  IP 

P8  -  0.0  . 

S  8  .  0.0 

THR  -  0.0  - 

0V8  -  0.0 

- T0-a  0>0 - 

SO  -  0.0 

DVO  •  0.0  "  - 


16100  XN  -  0.0 

16700  231  00  230  8-81,82 


16300 


16800 

T6600 


P8  •  PR  ♦  PIN) 

7135 - SR~-  S8  ♦  SIR) 

P8  •  pw/XNOP 
S8  -  S8/XN0P 


16700  238  00  280  8  -81,8? 

16500  C 


16900 

17000 


l  f  100 

17200 

17300 

17600 

IT500 

17600 

"T7700 

17800 

17900 

18000 

19100 

18200 

"moo 

18600 

18500 

18600 

18700 

18800 


4  O  71/U 

19000 
19600 
19500 
10600 
19700 
~T TflOO' 
19900 
7D000 
20100 
70700 
20300 


20500 

70600 

20700 

20000 

20900 


C  CALCULATE  POTENTIAL  TF*P  AND  SPECIFIC  VOLUME  ANOMALY 
C  RFEFRRFO  TO  PE. 

C 

2350  TH  f  8 )  .  THETA(P(M),r(M),S(M)yPE) 

OV(M)  «  0VAIPF,TH<M> ,SIM)) 

7352  thm  -  thm  ♦  rmn) 

OVM  -  OVM  ♦  OVtM) 

P  T  <  M 1  -  P!M| 

TTTwT”-”mnn - 

DVX  -  OVtM) 

C 

C  Fl,E2  APE  MIN  AND  max  SPECIFIC  VOLUME  ANOMALY’  WITHIN 

C  tnf  regression  interval. 

C 

TTTM-M 1  236,7767717 - 

736  FI  -  DVX 
F?  -  OVX 

737  IF (0VX-F1  1  7372,738,738 
737 T  FI  -  DVX 

738  IMF2-DVX»  2382,'239,239 


t  c  r  f  — 

739  CONTINUE 

TF( ABSI PI Ml-PFI-POIFF  1 760, 260,750 

C 

C  AVTPACF  T , S , 0 V  OVER  PF  ♦-  PDTEF - 

C 

740  TO"  '<i~TTT"*~TT  81 - 

SO-SO  ♦  SIMI 

OVO  -  OVO  ♦  071 M) 

XN  -  XN  ♦  1.0 

750  CONTINUE  - 

THM  -  THM/XNOP 


OVF  -  PVM 
TO  -  T 0 / XN 
SO  -  SO/XN 
DVO  -  OVO/XN 


21100 

71700 

21300 

71600 

21500 


21700 

21800 

21900 

77000 

22100 

"77700 

27300 

77600 

22500 

22600 

22700 


22900 

21000 

23100 

23133 


2503  CALL  LSFT 
C 

C  IE  TSSWI10J— 1  PRTirr  OUT  RFGRESSI08I  CUEFFTCTETTTTTOR  THIS  LEVEL 
C 


751  00  7510  M-N1,M* 

DVI  *  0VAIPIM1 , TIM), SIMI ) 

PTD  -  PTIM)  -  PM 
TTO  «  TT(MJ  THM 

7510  WRITE  (KLIST,?5lllM,p  f  m1,THIMI  xSIMI  ,OVI,’OVIMUPTD,TTt) 

- X - ,-71777 - 

7511  FORMAT! 1H  , 16 , E7. 1 ,2 F 7. 3, 3F7. 2 ,F7 . 3 , X , F6. 3 ,F 7. 6 , 2F3.0 1 
WRITE!XLTST,2515HrPIM3,Wil,Tn,PM 
WRITE!KLI5T,2515) (CTIMI ,M-1,N»,0VM 

7515  FORMATIIH  ,6Eri.51 


C  A2  TIMFS  THF  STO  OEV  71  I OEFAULT  IS  3>,  AND  EXCLUDE.  RE-EDIT. 

C  - 

753  rF IKFRR.GT.3 I  THEN 

WRITEIXLIST»253001  - 


25166 

23177 


25300 


FORHATllH  ,'KER*  IS  GREATER  tH«l  3-COIT  GIVE?  UR*  I 
GO  TO  7536 


23194 

7T7U0 

23300 

71400 

23900 


2  3700 
~7T750 
23775 
'71000 
23825 


IF (KERR  12  536,7532,253 2 

253?  00  7535  M-H1,*7  . 

IF  ( A8S( PTC  5 )-PM|-A247 1)2535,2534, 2534 

c  . . . . 

C  CALL  EDITING  SUBROUTINE 


2534  CALL  FOITIKERR) 

TF  I KERR 1 2  53 7 *<720 *2  20  . 

2537  TF(  ISSWm.EO.-llTHEN 

OVT -OVA I PF  ,TWFTOT  FT M>  ,T (Ml  , STM)  ,~PF )  ,  ?! Ml  ) 
WR  ITE (KL 1ST, 25370  IP TM),DVI 


23875  X  INTERPOLATION  OF  T  OR  S') 

T71F3  ENT5TF 

23891  GO  TO  220 

"23900  7535  WNTTNUF  - 

24000  2536  KERR  -  0 


24200  C  IF  OUTPUT  IS  REQUESTED  WRITE  DATA  SUFFER  KBUF  TO  FILE  ♦•REG 

75100  C  -  -  - 

24400  IF(ISSW(13) 1255*260,260 

75500  251  MRTTETKOUTirSTIF —  -  -  - 

24600  C 


24800 
- 24100 

25000 

"75100 

25200 


260  TFlTSSWCrZl  I765V267*767 

265  WR  ITE  (KL  rST,  2I650HREC*  PF  ,  T0*S?0*DV0  *DVN  VS*H,  THN,'71  *  7  7*N*NDP 
7650  FORMAT (IH  4‘*F7«TV7F7.  3,3F  7.7 VF7.1,X, 4*6.3, F7.4 ^7141 

267  CONTINUE 


25400  770  CONTINUE 

75500  280  TFT  TTSWf  I T1 1 255Y795V795 

25600  285  CLOSE (UNIT»K0UT) 

25700  WTO  106  - 

25800  FNO 


300 

400 

500 

“woo 

TOO 

— rfrr 
«oo 
<roo 
1000 
1100 
1200 
"T300 
1400 
1500 
1600 
l  TOO 
1800 


AVRCP  SUPPROC  POTEN  **♦♦*  A VfttfP  ******+♦*****♦*♦****♦****** 

♦•♦♦♦♦♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦I* 

- stjbpouti nf  mr? - 

*♦•■♦•♦♦***♦♦*♦♦♦*♦***•♦  4**  ***♦*****•*♦**♦*♦*♦♦  ♦*♦♦♦♦♦***** 

SUBROUTINE  TO  AVPRSCF  SPFCIFIC  VOLOWE  ANO  COETF ICIEMTS. 

JUNE  28  1976  N  EOFONOFF 

IFTTTUF  TITT  D!  TC  "CTO  28  FOtfWTT  18DECB0  N.‘  BRAY - 

INCLUDE  •COHPOTfN.FOP*  ' 

DIMENSION  - 

dtwfn-stun  tr  nwrn tttwctt iwTTCTOTnrnTmrwr - 

or MENS  ION  APF(100),$VI (1001 *EOB( 1001 
DIMENSION  VWlN(iOO>,VNAX(iOOi 

CHARACTER 


2000 

2100 

2200 

71013 

2600 

- 2500 

2600 

7800 

2900 

3000 

3100 


CHAR  AC TFR *  l  IVlVTV2fTV3 
FQUrVALFNCE 

FOUI VALFNCE  ( CR , IHDC 1 . ( PF ,VR l 

FXnj'T V ACE NC E ~ TCR  T  ITV 1 T5TTGT l«  l  CR  f  1  * 103  ,  WMVrCinTV  II  I  ^5VI  1 - 

FOUIVALFNCF  t  C*  ( 1  •  1 2 1  *  APE  1,  (  VM IN*CR  <  l  » 1 2  H  •  ( VM  AX  ,  CRI1  *  13  11 

READ  IN  STATION  »«S  TO  BE  AVEPAGEO.  ENCOOF  INTO  CORRESPONDING 
ETLF  NAMES. 


3300  RFWINO  12 

3500  OO  61  K-l,IOOO 

3600  RFA0(1?*  6 10»EN0«62  1B»FNAME(K1  *'WT  ( K  I 

3700  610  FOPWATUA*  A12**F5.71 

3900  61  CONTTNUF 


3932 

3050 

4000 

4300 

4200 


JSTN-K-1 

GNAME(9: IPI-t-AVG’ 

ENOIF 


AVPRAGING 


4400 

4800 

4600 

4700 

4800 

9900 

5000 

5100 

5200 

5300 

5400 

9500 

5600 

5700 

5800 

5^00 


100  ISW2  •  ISW  -  2 

rF(TSW2110 1.113  *113 

BRAfTCH  1  —  INTTIA1  AVER  AGING- —  BEG  TNT  HER  F 

101  00  110  J-l,iOO 

- 00-t  TO  ~T»1~*T3 - 

no  rRtj,n  -  o.o 

112  IRWX  -0  -  -  - 

BRANCHES  2  ANO  3  BEOTN  HERE  — 


01*PN  APPROPRIATE  FtlTF*  RTAO  HFT0E1T 
00  1200  KK«1*J9TN  ~~  - 


1 


% 

'4 

'I 

TK 


6000 

6050 


6100 

6200 

6100 

6600 

6500 

6600 


OTWC  tf  NT  T-KOUT  VNAMF-FNAME  (KKIVRTAWNLY.TY  PE-mD  *  V 
*  EORM-^UNFOPMATTEO*) 


H 

NSW-5 


C  CHECK  IF  OAT*  S FLECTION  PARAMETERS  ARE  SATISFIED.  TF  NOT, 
C  ~  NSW  IS  RE TOftNTTT  PRffM  OITA  AS  6,  AW5TATtON  ! $  SKIPPED. 


rF(NSW.NE.5)G0  TO  1200 
NST»NSt*l 
OELB-O.O 
P  PR -0.0 
82-0.0 


B3-0.0 


8000 
-i  DTOO 
B200 
8100 
8600 


.MiM 


8600 

8700 

8800 

8000 

5000 


5200 

5600 

550IT 

5600 


10000 
TO  TOO 
10200 
lOTUO 
10600 


IU7UU 

10600 


C  IF  TSSWI151--1  WETOHTS  ARE  TAKEN  FROM  JSHF.PTN" FILE?  OTHERWISE i 
C  THEY  ARE  SET  TO  l.' _ _ 

TF I tSSW 115)11350, 1357, 1357 


nu  i  —  R  i 

GO  TO  1370 

1357  WGT  -  1.0  " 

1370  TF  < ISW2  )  1 20, 120,138 _ _  _ 

C  BRANCHES  1  AND  2  CONTINUE  HERE  FROM  STATEMENT  #1370 


120  READ! KOUT ,ENO«T60 IK6UF 

GO  TO  160  *"  '  “  "  . 

C 

C  BRANCH  3  (  WRITE  OUT  AVERAGEO  FILESI  tONTINOFS  HERE  FROM  STATEMENT  #1370 • 
C 


I  iiiivv&i  kiiiKv/i  rimriii 


135  CONTINUE 

C  OPEN  NFW  FILE  NAMEO  *. AVG  CORRESPONDING  TO  INPUT  R.REC*  FOR  OUTPUT 
C  ~  TIN  WTHT.~  —  - - 


LrL.ur.4i 


10800 

T050IT 

11000 


I  t»ii] 


11200 

TITOO 

11600 

TI5O0 

11600 


OPEN! UNI T-ll ,NAME -GNAMF ,TTPE- *NEN  * , FORM*  * UNFORMATTED* ) 

C  WRITE  HEADER  TO  *• AVG _ 

WRITEU1IKHDC 


X  BRANCHES  1  AND  ?  CONTINUE  HI 
C 

160  TFT  TREC-TRNXTI5B  VT  67V 1 6  ^ 
162  IRMX  -  IREC 


X  TTRANCH  2  XTJNTTNUET  TfE5fE~FR0M“ 
C 

T6T  TTONTINUF  -  - 

1675  OEL  -  SVA(IREC)  -  OVM 


11800 

CPM(l)  •  CPMU>  ♦  PM  -PE 

TT50U 

12000 

00  150  J -I , N 

150  CR(IRFCtJ)  -  CRMRECiJI  ♦  WGTPCPMIJI 

T7T0XF  X 

1 

127U0  C  3RANCH  2  RETURNS  TO  RE  AO  NEXT  OAT*  tECQRD  EROW "*«W¥t 
12300  C 


12500 

176W 

12700 

"mroo 

12900 

nnnnr 

13100 

"13700“ 

13250 

13300 

13500 


13600 
13700 
13800 
13900“ 
16000 
_TMUO". 
16200 
16300 
16600 
1550 0 
16600 


C  BffTNCH  1  CONTINUES  EDO*  S  TATENENT  #166 
C 

155  CONTINUE  —  -  “  - - "  '  ‘  1 

158  SVAIIRFC)  -  9VAIIRFC)  ♦  WGT60VN 

WTN7TRTCT  -  ^NTTRED  ~^CT*E1 - - 

V«AX<IRFC>  -  VNAXITREC)  ♦  WGT6F2 

SVKIRFC)  -  SVIHREC)  ♦  WGT5DV2R0I  PF*DTN*PN*CP»N*NnP*FlVf  2,'  j 

♦  ISHP,KCASTtICON,OELP )  j 

SWCT(TRFC)  -  SWGTITRECI  ♦  WCT -  ] 


CO  TO  120 

160  CLOSE  I  UN  I T •*  OUT  1 
CLOSE(UNIT-Il) 

_ IE (KK.LT. JSTN1G0  TO  1200 

TTTTSV7 


C  BRANCH  1  CONTINUES  PRO*  PREVI01TS“  STSTTHENT 
C 

161  on  1605  J-IVIRH*  ~  "•  . 

VNINIJI  -  VNINI JI/SWGTI J> 


1 


16800 
T6900 
150 00 
15100 
15200 


JUU 

15600 
15500 
15600 
15700 
15800 
15900 
16000 
16T00 
16200 
1 6  TO  0 
16600 
“T6507T 
16600 
16700 
16800 
15900 
17000 


SVKJl  «  SVM  Jl/SNGTf  J) 

1605  SVAIJ)  «  SVAIJ) /5VGTIJ)  - 

TF(  I5W?H0O,18O»T0O 
C 

C  BRANCH  2  CONTINUES  f RON  S  T  ATENFNT  PRECEEDIW?  »161 


iJ JTrfj 


1610  OFLR  -  0.0 
PPP  •  0.0 
KF  -  ? 

N  «  NSCl NSFCTION+l 1 


IEIlSSWfll).FQ.-l)WPITE(KLlST*16230) 

no  161  j«i , rwMX  - - 

^  IFUPRI  JfNPP,NSPCTI0N)-NSCfKF))I612«16Il»1612 

NDP  -  NSC  I KF ♦2*NSFCT ION) 

- ^—a  T - 

1812  DO  162  I » 1  *  N 

162  CPNIT1  •  CR  I  J»  T  ITSNGTf  J1  -  - 

PF  -  JPRI J,NPR,NSECTION» 
rFiissNUinifimri576*i625 
1621  WRITEIKLIST*16?3)PF»  (CPfHI)  t  I«l«N) 


t 


17200  16210  F0BMATI1H  ,'AVERAGFD  REGRESSION  COEFFICIENTS*  *■*/* 

17300  *  ,PPESSORE^»2XVrrCPITl  %5T,-rCFr27Tr5Tt^CPTl3T%5TV»m6)%  “  - 

17600  ♦  5X»*CPI5)»,5K,**CP<6)*) 

17500  1626  OVT  -  SVAIJ) 

17600  OVE  •  DVZRniPF,DVI,PF,CPH*N*NOP*VNINf 7V,VNAXf J»*15HP,KCAST*t 


l rruu  UtUB  ■  iTIIJf  -  uv^ 

17800  SVAIJ)  -  OVE 

17900  0FLP1  «  PF  -  PPR  -  ' 

17933  EOBIJ)  -  r)POV10VFfOVI*CPN*N*VNlN(J)*VNAX(JI) 

17965  EOBIJ)  -  l./SOBIJ)  -  - 


74000 

24100 


TYi-^y* 

rv2-»r* 


89 


— rrrtro — 

— rry«T» - — — 

24300 

N1  •  l 

“TRW 

N?  •  7 

24300 

N3  •  3 

24600 

ITT"*  6 

24700 

NV?  -  4 

NV3"-_"8 

24900 

WRITE(KTTX,3l7D)Nl,N2»N3,KT0*IYR,TT8,m-,IV7,lV3 

74000 

R  (EIOC  KlN,  ♦  INI',N2,N3,KT0,IYR,IT" 

29100 

REA0CKIN,3175)1V1, TV2,  IV3 

74700 

WR  TTTf  RTTX,  3 176  VNV1 ,NV2 ,NV3 

24300 

READCKtN,*)Nfl,'NV?,NV3 

— 74400 — 

T179 

— Fmmrrmm 

29300 

3176 

F0R"ATI1H  ,,NV1,NV?,NV3',3!5) 

75600 

ENOfF 

29700 

JSW  «  0 

2580XT 

C 

29900 

C 

N1 

TO  N3  ARE  EFFECTIVELY  IGNORED*  UNLESS  ISS3M81— 1.  SEE  STATEMENT  ■490. 

T 

26100 

320 

IFCTREC—N 1)321 V4 90 ,321 

r«?oo 

321 

IF  CIREC-N 2)327 *<4  70 *322 

26300 

322 

IF ( IRE C-N3 ) 120, 324*120 

'  76400 

774 

OUNTTNUE 

26900 

C 

— 2'6B00 

~ c 

— TFT 

VALUES  FRU8  ""AF  FURPA1  VARTAflLTS 

26700 

c 

26800 

3243 

VR1  •  VR(NVl) 

26900 

VR 2  •  VRCNV21 

27000 

V»3  •  VRCNV31 

27100 

c 

77700 

~z 

— nr- 

27300 

c 

OMFCK  PURPOSES  only). 

27400 

c 

2  7300 

329 

IFCISSWf 101  1 376, 327, 3 27 

77600 

376 

WRTTEIKLIST*3T50nCDN,!RFC,!SMP,TCAST,»F,PI,T0,VRl,'VR7,VR3, 

27700 

X 

OH, PE 

— mnro — 

27900 

c 

c 

IF 

9AP  FOR"  AT  NOT  RFOUFSTEO,  RFTURN  TO  ST*TF"ENT  *120  TD  READ  NFXT  DATA 

28000 

c 

RFCORO. 

28100 

c 

28700 

c 

888 

F0R8AT  WRITTFN  TO  UNIT  KTO 

28300 

c 

"377 

28900 

330 

LFTLE  •  0 

28600 

I PF  -  PF 

28700 

XLG  -  XLONG 

28800 

339 

WRTTE(8T0,3300)CL«9L  !"),"•  1,21, TSHP,TC AST, LF ILF, TPF,1 

28900 

X 

XL  AT,  XLG,  I YR  *J0 A Y , IT" , I V 1 , VR 1 ,1 V2 , VR 2 , t V3 ♦ VR 3 

29000 

IFCLF1lc»390, 12 0,990 

29100 

3190 

F0R8ATUH  ,I3,X,2I7,  I4,X,2F5.0*2F7.3,2F8.3,F6.3»F7.2) 

29700 

3300 

FOR"* M 2*4, 11  2T4,X,Il,T5,F7.2,F8.2,X,T?Vt4,T2*7,3C»T,F9.4n 

29300 

c 

29400 

c 

ITTVm—  1  ALLOWS  SPECIAL  FUNCTTOWT  TO  8F  CO8PUTF0 — 3U8TPACTTH0  VALUES 

29900 

r 

AT 

ONF  LEVFL  FRO"  ANOTHER  BEFORE  OUTPUTINC  TN  "AP  FOR"AT. 

29600 

XT 

29700 

490 

IFC  ISSWC8)  )455*.3?4,3?4 

29800 

435 

V»l  -  VRINVt  l. 

29900 

VR?  -  VRINV2) 

30000 

V«3  -  VRCNV3) 

70100 

30200 


90 


- GO  TO  120 

470  TF  f  ISSN<8>  )4?5V324i324 

TO  TOO - 4T5 - VRT  =-TTCTimT  ~  V*1 - 

30400  VS?  -  VR(NV?K  -  VR? 

7P0700  VR3  -  VRTNV3I  -  "VR3 

30600  GO  TO  325 

70700  C 

30800  C  BRANCHES  1  AND  2  CONTINUE  FROM  183  OR  180 
30000  ~~  T  -  . . 

31000  C  IF  1SSHI7)-0  CONTfNUE  THROUGH  BRAHCH  3  AUTOHAT  I  CAttT. 

31 1 00  T  "  . .  "  '  . 

31200  400  IFIISSWI  7n550f’410|410 

moo  mo  i th  *  rsir  ♦  r  ~ 

31400  IF<ISW-3>100. 100,420 

”71500 - C - 

31600  C  BRANCH  3  CONTINUES1  HERE  F RON  STATEHFNT  PRECEEOING  R161 

7T700  C  . .  ” 

31800  C  IF  NAP  OUTPUT  REQUESTED  SET  LFILE  TO  1  TO  INDICATE  EOF  T"N  HAP  FORHAT. 

71000  -  C . .  ~  . .  '  ■'  "  . 


32000  420  lF(ISSW<14n422V560,560 

I7ITH? - 477 - LFILE'  ■  T - 

32200  GO  TO  335 

32300  550  CONTINUE 

32400  1200  CONTINUE 

37450  ”  I F ( NS  H .  Nf  . 5 . A  NO . T5V. E  O . 1 160  T0  16 1~ 

32475  IFCNSW.NE.5.AN0.ISM.EQ.2IC0  TO  1610 

32487  rF(NSN«NE.5« AND. tSN.E0«3)G0  TO  420 

32500  NftO  RETURN 

—moo -  ENO 


700 

800 

900 

TOOO 

1100 

TTOTT 

1300 

T600 

1500 

1550 

1600 


1800 

1900 

2100 

T20O 

2300 

"7W 

2500 

2600 

2700 

2500 

2900 

-3000 

3100 

3200 

3300 

3500 

3500 

1600" 

3700 

3500 

3900 

6000 

6100 

6200 

6300 

6600 

6500 

6600 

6700 


C  DAT f  SORR  PATEN  *****  OUT!  **«**«»*»* 

- SUBROUTINE  0 mrnrOITY NS tfl - - 

c  **»*♦**♦♦*♦***♦*♦*♦*♦♦*♦*♦*#♦*♦**•♦*** 

c 

C  TO  SELECT  ANO  ACCESS  CT078  FORMAT  OATA  FROM  VAX  0I5C  FORMAT 
C  '  ATtFSSFS  VARIOUS  WTTLXRD  SUBROOTTTfFT  EOUHO  TTT  CT0ITA7tT5 
C 


C  JAN  6  1976  N  FOFONOFT 

C  ROUTE TEO  FOR  CT078  FORMAT  INPUT  TU  OEC  50. 

C 

INCLUDE  *CO*EOTEU.FUr* 

C 


N  BRAY 


VIVIM  ilia 


IIII.U.  koaii. 


kif  tia.iai 


TNCLUDE  MOXREC.'DT'M’ 


C  CHARACTER 

C 


PROGRAM 


IF (NSW. FQ. 5  I  GO  TO  RO  '  ~ 

IFIN5W.EO.2I  GO  TO  30 

- -TEtirsurrvro'.u - 

NSW  LESS  THAN  ZERO'*  INI T I ALIZE  SELECTION  PARAMETERS 


CONTINUE 


JOO  - 
TJAYI- 
0AY2 
XEMN 
XFMX 
XNMN 
XNMX 
—  XtTO 
XL  CO 
JSTN 
RETURN 


0.' 

365. 

-1B0.U 

180.0 

-90.0 

90.0 

ro~;oo~~ 

70.00 

1 


NSW-O:  LIST  OR  CHANGE’  SELECTION  PARAMETERS 


6900  20  CONTTNUF 

5000  IT?  WRITE  naiST.173IDATI  fDAYrrJU0 

5100  173  F0R"AT(1H  .5H0AY l  : F8 . 3 * X, 5H0AY2? . E 6. 3,'X*6H  JOO?  *1  6 » 

5200  REBDfKTNt* 1 0ATT*0RY2»  JOO 

5300  176  WRITFIKLIST*175JXEMN,XEMX,XNNNVXN*'X 

- 3500 - TT3 PUR  MATT  tlT"  f7WT*N  L I M  »  6E  7  .T1 - 

5500  REA0CKrN,*lXEMN»XEMX,XNMN,XNNX 

5600  WPTTE(KLIST»17TtXLTO,XLCO - 

5700  17 7  FORMAT ( 1 H  *8H0R TGIN*  «2(X*F8.3)) 

5500  RFAO(KINf*>XtTOf-XLGO 

5900  RETURN 

- 5000 - C -  -  — 

6100  C  NSW  -  ? :  READ  FROM  FILE  STATIONS. PTN  INFORMATION  TO  IOFNT TFY 
6200  C  STATIONS 
6300  C 

6600  30  CONTINUE 


6450 

6600 

irscnr 

6607 

6614 

6650 

6667 

6674 


TV 

& 


4 


*1 


y  ((iij 


7150 

7700' 

7300 

7175 

7337 


7500 

7600 

7700 

T¥00 

7900 


11300 

'¥000 

9200 

“7300 

9316 


■mj 

9366 

'9400 

9500 

9600 

9700 


9900 
TOO  00 
10100 
10700 
10300 


10400 

T0500 

10600 

roToo 

10600 


L’J'tUU 
11000 
I  IlOO 
11200 
11700 
1 1400 


:t 

g 


11600 

77709 

11750 

TTWOO 


If  TK8R .EQ. 1 3) THEN 
WRITEfKTTX,3l01 


MRlTEfKTTX,320> 

REA0<KIN,330mHP,rcRUlS,TE*Drj 
300  FORMATfA) 

310  F  Of  HAT  (  I  H  ,  *  ENTER  SUBDIRECTORY"  7F9TI  »•  1 

320  F  OR  MAT  f  1 H  ,  *  ENTER  SHIP  C00E,CRtnSE  *,PRTJJ  «*! 


- 6¥¥F  ' 

- TTO - FORMT?  14  2,2131 - 

6693 

ENDIF 

6700 

C 

6800 

6900 

C  "ILLARD  HEADER  RELATED  SUBROUTINES 

c  "  —  ”  - 

- . .  -  -  - 

- -  - 

7000 

CALL  PVFRIPROVERl 

Milil.fi  liiiP  il 


CALL  STATION* OVO*KTP J 
CALL  INDEX  fill 
LRFC  -  IOXP.fCf-91 
LLREC  ■  LRFC 
KKST  -  0 


n-v 

RETURN 

C  NSN  *  IX  RE  AO  STATION  HEADER*  CHECK  AGAINST  DATA  SETE'CTIUN  CRITERIA, 
C  AND  RFAO  TFMPF 7ATDRE  ANDS  AL INTTY TNTDDIYIARRAY. 

C 


iii.'kAl.'Ki: 


rFIKUN.GT.LLRECl  GO  TO  620 
CALL  RECTDXIKON) 
XLAT-SLATU 
XLONG-SLNG ( ) 

IPR  -  PM  I N 


L  PR  -  X  N  •  P  R  S  I  NT  VP  M  TN 
LBBLlll-IEHEOf 31 
ENCODE  I  4,53, LBBLf  211 IFHEDf  4) 
ENCODE  f  4  ,9 4  ,L¥RUm  TTEHE0T5T 
FORMAT* I  3, •-*  » 


C  CONPUTE  JULIAN  YEA*  DAY 
C  CHEEK  AGAINST  SELECTION  PARAMETERS 

c  ■  —  '  - - - 

I OAY-KDAY< TOA/TMO, IYR l-KDAYf 31,12, IYR-1 A 


ItMIBIIf  i 


TPfOAY-OAY1162Di602,602 
607  IF  IUAY-OAY?  J6O4  ,f604,’6  2O 
604  CONTINUE 


C  CHECK  LAT  AND  LONG  AGIANST  SELECTION  PARAMETERS. 


L 

TF  fXL0NG-XEMNI620,606,606 

406  IE  f  XLOWG-7EWX»60D,6DB  ,  670 - “ - 

608  TE ! XL AT -XNMN 1620,610, 610 

610  TEf  XLAT-7NMTI6T6T6T¥i  670  - * 

616  L  TYPE  -  1 


l  v.un  - 

DAY  -  DAY  ♦  TDD 

KCAST  •  I7TA7  - 

JPM AX  •  NTOT 


I 


4 


rt9oo 
12000 
•'intro 
12200 
T7T00 
12400 
17500 
12600 
'1770  O’ 
12000 
17000 
11000 
1 1100 
11150 
1T70D 
11100 
11400 
13500 
11500 
13608 
17616 
13617 
11618 
11620 
13627 
11624 
V1676 
11628 
15612 
11640 
11648 
11664 
~T76W7T 
1 3690 
11700 
13800 
11500 
14000 
T4T00 
14200 
14100 

14400 

14500 

1  4600 

1470  0' 

14800 

14900 

15000 

15100 

15200 

rsimr 

15400 

15401 
15407 
15410 

15414 

Twr 

15428 

15415 
15438 
15442 


93 

C  IF  T5SUI111  — 1  WRITE  OUT  HEAOER  INFORMATION  ON  UNtt  KLTST 
C 

- C1S5VK  1  TTT6T6  O  ?  6'7Tr,T>  70 - 

6160  WRITEIKL  1ST  ,  8000  )  I LBBUK  )  .«• 1  v3 1  «  XL  *T  VXLONG,  XITO  ,XLCO 

WRITE  (KLIST.MO'HLTYPE.HHDIUTS'HP.TSTAS.  rCON.OAY,  TP* 

C 

C  IF  YSSW (12  3 «-l ,  W*T‘TP  OUT  HEADINGS'  FOR  OUTPUT  WHICH  IS  WRITTEN  IN  CUN»S. 
C 

- TFTTOVI  1 2  1 7ETJY- 1  HJff nTTKC  1ST* 8010 ) -  ~ 

8010  FOB*ATI  •01REC  PE  TO  SO  DVO  DVH  S"  TH*» 

V  7 1  Z  2  N  NDP  •  ) 

8000  F  OP  HAT I  /  «  1H  *3A4»4F8.3) 

8005  EOPHATdH  ,  *  TYPF  *DHR  SHIP  TSTN  CAST  DAY  ?PR»,/. 

♦  215, 2**A2, 27,215, F8.3»2I5» 

— NOTTS'" — F"PB  W  A  T  ITS  ',T  %  IT  V,  T7,I  A.T.TS.rS  ,V7.7  ,Fg;  ?  ,7  ,  T  2  V  ITTET.l ) - 

C 

C  MTLLARO  SUflPO UTINE  TQ  FILL  DATA  APPAY  WITH  TEMP  AND  SAL  TNI TY 
C  DATA  FOR  ALL  OBSERVATIONS. 

C 

WGT-1. 

- rrrrTTSTTTT.nr:  -7  jtwfn - 

*.*♦1 

JTCR-lCRUlS 
JIST-ISTAS 
CALL  L2IJTCR1 
CALL  LZ( JtSTI 

- - — - ENCODE  X 1  2,5 2TGNTWF  T T  SKP  V  TT CR  7  JTS  T - 

5?  FORMAT! A2.A1, 41, ‘.RFC* 1 

WRITE ( 12.80?0)n,GNAME,WGT 
8020  E0RMATI14.A12.F5.2I 

OPEN! UNI T "KOUT ,NAHF»GNAHE»TYPE«,NFW,«FORH»  • UNFORM a TTEO* 
WRITEIKOUTIKNOG 

- EWTE - - — - - - 

CALL  DATI0XUCUN1 
CALL  GFTDAT(FTP,DATAX»3300,2f 
6 ?0  RETURN 
C 

C  NSW-5:  CHECK  ONLY  LAT  ANO  LONG  OF  HFAOFR  ALREADY  RF AD  ACAJNS 

T - SFCTC7  JON  "PAR  S*FTPTS  .  - '  '  ' 

C 

80  TF1VL0NG-XFnNM5,87,87 

87  IF  1 YLONG-XENX 1 89, 89, 8  5 
89  TF f XLAT-XN8N)85v83,fl3 
83  rFCXLAT-XN**X)8?,82,85 

8?  CONtl nuf  "  "  ■'■■■  .  ■' 

RETURN 
C 

C  NSW«4  IMPLIES  SELECTION  CRITERIA  A«DVE  NOT  NET. 

C 

85  NS  W»4 

-  R'FTURN  ‘  ”  ‘ - -  " 

END 

C  SLAT  function 

REAL  FUNCTION  SLAT 

INCLUDE  MDXREC.OIM* 

T  TUNCTPON  RE TUFN5  DECIMAL  - UTTTRFE  VALUF "STONFO  -  FUR  SOUTH  F  WF5Y 
XLAT-RTSO 
XLATN-TLTSN 

XL  AT*. XL ATn/6000. 

SLAT.XLAT»STXNTTLAT»,XLAT1 
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VAX-11  PROGRAMS  FOR  COMPUTING  AVAILABLE  POTENTIAL  ENERGY  FROM  C— ETC(U) 
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UNCLASSIFIED  WHOI-81-70  NL 


*1 


94 

15449  RETURN 

15456  ENTRY  SLNC 

15463  XLAT«ILNS0 

15470  XLAT*«ILN$M 

I54T3  XLATK- XL  ATM/6'000. 


15477  $LNG*XL  ASSIGN  EXEAT'S  XL  AT  I 

- 1 54*  4  - RTTUEW  .  ~~  - 

15491  F«m 

- 15500 - C"***-***********  SUBKWJT IWF  LZITA1  WWWTWTTVVT***VVT*T*T*Trv*TTrr- 

15503  SUBROUTINE  L 7  f  T  A 1 

- 15504  C  *  *  *  *  ♦  * *  *  ♦*♦*♦*>**♦44*  WWW** ********  *' *444*  W¥9WfW*f**W9  *  44“ 

15506  TNTEGER  TAin.IWCl) 

155 12  TFTT1T1)  .GEYTOUY~GU  TUTCO  - 

15515  TFtlAm.GE.10)  GO  TO  10 

- rr5rg - YFTTTTT  JTCF .  0  J  770  TO  I - 

15521  RETURN 

15524  lOG  CONTINUE  —  "  “ 

15527  ENC00E<3«2*  TMlin  IAC1* 

"15550  - TATlT-TVflT -  -  - - 

15533  2  PORNATIim 

- 15,3b - RFTTTRN - 

15539  10  CONTINUE 

15542  ENC00F(3,l.T«(tn  IAfir  — 

15545  TA(1I«IW(1I 

“175548  TTWTNTnTHO  *  T  21 

15551  RETURN 

- T5554 - rTTTRTTNUE - 

15557  ENCOPE (3«4«IN(1))  lAfl) 

T5560  4  FORNATI2MOOVTH 

15563  IAm«IW(l) 

155*6  RETURN  ~ 

15584  END 


- 1 


J 


95 


12000 

12100 

12600 

12700 

tttoo 

12900 


c 

c 

*r 

c 

c 


POTFXJSUB.FOR  FILFs  SUBROuriNRl "ft»  FOTEN,PFPl.T  POTENTIAL  fNFRGY 
PPOORAHS.  VAX  VFRStON.  N.BRAY . 

SUBROUTINE  S*INVf  l,N,P  ,0,«,!F»IL  ) 

*#  ***********  *****************************************  ♦**♦♦♦ 


T3000 

c 

TO  T5VFRT  SYW0FT9TC  41TRTX  FOR~  T  R I  AffGUl  AR  SFCTTON  IRtTINOFO 

13100 

c 

TN  A  LTNFAR  ARRAY  MJ). 

13200 

t 

PRO*  5YHTNV?,, ,CAC*  *Y50  BY  RUT15H10SER  VTA  j,‘  NTtTTTTTT 

13300 

c 

IT40  0 

c 

AFRTL  27  1975  N.  FOFOVOFF 

13500 

c 

nnsoo 

01 0FNS  I ON  Ain  *iPTD  *  Of  1  V*MR  ( 11  — 

13700 

c 

13000 — 

- tmt-r-r 

13900 

00  10  I • 1 »  N 

- I4TTOIT 

rcr  0ritj  «o 

14100 

c 

5F/TROH  FOR  PIVOT 

1 4700 

TO  100  T-l ,N 

14300 

8IGAJ  -  0.0 

73  «  ^ 

14500 

00  20  J-1,N 

14600 

JJ  -  JJ4N-J*2  - 

14700 

B  -  ABS(A( JJI ) 

14000 

TFTWRf  JI120»T7,7O  - 

14900 

12  rF<B-BIGAJI20,20Vl4 

mi  1 1 1 1 1 1 1 1 

1 4  B  I G  A  J  “  IT 

15100 

K  «  J 

- 157110 

K*  -  JJ 

‘ 

15300 

70  CONTINUF 

15400 

IFIRTGAJI16, 15*16 

15500 

15  TPAIL  •  1 

1  U^tf 

15700 

c 

PREPARATION  OF  FLINIYfATION 

T50OO 

16  NRfK)  -  1  - 

15900 

QtK)  -  l./AfKK) 

16000 

KJKI  -  1.0 

16100 

AfKK)  «  0.0 

TflT 

Tt^T  *  K-l  -  - 

16300 

TFfK«lU5,  19,160 

16400 

160  JK  -  K  -  N  ~  ‘  " 

16500 

00  30  J» 1 , K N 1 

16400 

JK  »  JK4N-J41  - 

‘  * 

16700 

P«  J)  «  Af  JK  ) 

m%  i  H1  m 

!Ff*Rf  J1  )18,I7,V0 

16900 

17  01  J)  -  -Af  JK)*Qf«) 

17000 

00  TO  30 

17100 

10  OfJ)  -  A(JK)POfR) 

17700 

3TT  Af  JK)  -  0.0 

17300 

19  KP1  *  K*i 

17400 

ic  j  ^  icic 

m 

17500 

IFf  KP1— N) 71, 71. 41 

17600 

71  OO  40  J»KPjfV  '  ' 

. . i 

17700 

KJ  -  KJ  ♦  1 

17000 

TFf*RTJT)34, 37/34 

17900 

3?  *MJ)  -  A ( K J ) 

mono 

on  nr  75  ~  • 

10100 

34  PIJ)  -  -AfKJ) 

10700 

35  Of  J)  -  -AfKJ)*OfK)  “ 

* 

10300 

70400 


40  AfKJ)  - 
ELtWTWATTON 


0.0 

PROPFR 


ist*l 


4 


'! 

$ 


11700 

18600 


*1  JK  -  0 

00  80  J-l.N 


till'] 


10800 

“T1900 

19000 

T9T00 

19200 


bUB.lt] 


JK  -  JK  ♦  1 
50  M  JK 1  -  JK) 
100  CONTINUE 
1 5TJ  FFTWN 
END 


♦  PfJ>*0(KJ 


10400  C  POTENT I AL  ENERGY  PROCWAN 

~T9Wo  c  "  - 

19600  C ******* ******* ***************************************** 

19700  SUBROUTINE  PAR  AW 

19800  C*** ********************  *****************  *********  ****** 


20000 
"70TO  0  C 

20200 

“2rnw  C 

20400  C 


20600 

20700 

20800 

20900 

21000 


INCLUDE  *  CONPOTEN. FOR  * 


INTEGER  PFOELTA ,PDELT  A 


icim  i  *  *  ; 

FORHATI ’OINPUT  NFM  PARANETERST* » 
T  F  f  NOY  ESI  K  IN  ,XT  TX  J  .Tit .  UR  FT  OWN 
HR  ITE(KTTX,20» 

FOPHATf’OENTFR  THE  NUHBER  OFSEC’ 
THE  SAME  PARANPTERS) * J 


OF  LE'VFLS  WITH 


21200 
2T300 
21400 
— 2 1500 
21600 


C  l  fVV 
21800 
21900 

22000 
“77TI TO 
22200 


22300 
22400 
7770  IT 
22600 
27700 
22000 


C  C'fUU 

23000 

73100 

23200 

23300 

23400 


C  31UU 

23600 

“27700 

23800 

27900 

24000 


C7  1WU 

24200 

77700 

24400 

74500 


TF  f  NSECTION. GT.19 ) THEN 
MR  TTEIXTTX *  22  )• 

FOR  NAT  (  *0NAXINUH  ALLOWED  IS  19*1 

CO  TO  IS  - - - - 

END  IF 


FORNATI’OENTER  THE  PRESSURE  FOW  THE  FIRST  LEVEL** > 

R  E  AO  (K  IN*  *11  NT  TTAL  P -  - 

NLEVP-i 

NPRFV--1  “ 

NSE«NSECTION*l 


00  1000  I »l*NSE 
I2-I*NSE 
T3“I*7*N3E 
HRITE(KTTX,iOO>r 

F OP PTATl'rO FUR  TECTT ON* *T  4  ENTE*  THE  INTERVAL  IN  1 
f  LEVELEO  SURFACES*  *  I 


jn*  *  i  r  p  urT.  1  m 

WRITE (KTTX*120) 

FORMAT! ’CENTER  THE  TNTEHVAL  SITE  IN  1>B  FOR  THF9EI 
READIKIN** ) P DELTA 

HR  TTE  IXTTX* 1401  -  - 

FORNAT I  *  OENTER  THE  FIRST  PRESSURE  IN  THE  NEXT  SECTION’ 


01  OFTMFEir 


NLEVEL* I IP2-INITI ALP  ) /PFDELT A 

WRITE  (KT3Tr*160T  - - - - 

160  FORNATI ’OENTER  THE  NUN8ER  OF  TERNS  IN  THE  REGRESSf ON* • , 7* 

T  MN-2  IWFLTES  A  LTNFIR  FTTTXAXINUN  TTTT  RT’l - 

READIKIN, *)N 


ir»  n«n(  «o 

c******  COMPUTE  NPR(I)  ♦•♦♦♦* 

NPR  C  T1 • NL  F VEE ♦NUFVR  '  -  - 

C******  CONPUTF  NPR 1121 

NPRfI?r«7F0ELT»  “  ~  -  “ 


l 


'Z45TOO 

24700 


C.  lOWU 

24900 

'73000 

25100 

752UO 

25300 


25600 

25700 

25800 

25900 

26000 


t  O  1UU 
26200 
76300 
26400 
76450 
26500 


Z6600 
26700 
2  6  BOO 
26900 
27000 
27100 
— 27200 
27300 
77400 

27500 

77500 

27700 
TTfOTJ 
2  7900 
20000 
28100 
70  200 
28300 


c 

28500 
20500 
28700 
2  8800 
28900 


29100 

29200 

29300 

79400 

29500 


C  TOUU 

29700 

79800 

29900 

30000 

30100 

30700 

30300 

30400 

30500 

30600 


C«****t  COMPUTE  NSCIfl  *♦♦♦♦♦ 
TFCI.FQ.il THEN 


J  V<  l  1  I  ■  1  Ti  t  1  l  A 

IPREV-INITI AtP 

VNDTF  '  - 

NSC  Cl  *1 ) ■ IPREV»NLE VEL*PFDELTA 

cr*r»#9  CU8PUTT  NPRT TTT - 

NPRf 131-NLFVP-f TPREV/PFOELTA) 


nrnn««i>n  l  l  I) 

IPREV-NSCf 1*11 
NlEVP-NPRITT 
INITIALP-IPREV 

C**»**»  CONFUTE  NSCIT21  . . 

NSCU21-N 


NSC! I31-P0ELTA7DELP 
1000  CONTINUE 

IT0TAL-IPRFV 

jmax-nlfv» 

WRITE! KTTX, 200 INSECT  ION, NLE VP, I  TOTAL 


X  L7VEL  IS  A’»'*yl6*«0B.  •  ) 

WPITE!KTTX,220) 

??0  FORMAT!*  0 ENTER  MAXIMUM  DEPTH  OF  THF  DATA:*1 
RFAOfKIN,*) ITNAX 
I *NSECT I 0N>1 

- rr*r*x - 

n«3*T 

NPPfI)-NLEVP*5 

NPRf I21-! !TM*X«5001/INlEVP-NPREV> 

NPR  f 1 3 1 »NPRFV  -  " 

NSC! II-IPRFV 

- RTCTT7TS3j - 

NS C!I3)*P DELTA  70ELP 
NSECTTON-NSECTTON*!  " 

RETURN 

FNO  - - — . 

FOEPF  SUBR  PT5B1  ♦♦♦♦♦♦♦*****♦*♦*♦*******♦♦***♦•*♦♦♦ 


c  •  ■  . . 

C  CONFUTES  COEFFIC  TENTS  FOR  A  LINEAR  TRANSFORMATION  X*AX*8 
C  FOR  POLYNOMIAL  OF  ORDER  N-l,  INPUT  ARRAY  C*  OUTPUT  0.' 

C 


liiAi 


dimension  cm ,.om 

DO  75  I-l.N  ■ 

»  •  1.0 


WNf  »  N  -  I 

TFIN8J ) 17*1?  1 5  '  - 

5  DO  10  J-1,N*M 
TP J  •  I  ♦  J 

R  •  »FLOATIIPJ-«il/FLOATI  J))4B4R 

to  5  -  5  ♦~R»cmm - 

i?  tmi  •  r  -  i 

TFIT81I15, 15*20  - 

15  R  -  1.0 

00  To  26  -  - 


I 


30700 

3OB00 


iTitTtin 


31000 
711  Off 
32300 
37400 
32500 


JflUU 

32700 

37800 

32900 

77'OffO 

33100 


iUiIi] 


33300 

17400 

33500 

36900 

37000 


37200 

37300 

37400 

37500 

37600 


20  p  •  A*4r*fi 
25  0U1  -  S*R 


®IN>  FCN  44444  PTSB1  *♦♦♦♦♦♦ 
FUNCTION  PIN)  - 


INTEGER*?  IA,N 
FA  -  "77777 
?F|NI1,2,2 

P  -  FLOAT!  (  HANOI  N,  mi  *32 761. 


RETURN 


lllIU 


-w 


37800 

37900 

37950 

37975 

38000 


38200 

78700 

38400 

36500 

38600 


38700 

38800 

38900 

39000 

79T00 

39200 


IfJUU 

39400 

30500 

39600 

79700 

39800 


ilO.Mil 


40000 

40100 

40200 

40700 

40400 


RETURN 

FNO 


C  4* 44*444 44444  44  4  4 44444444' 
SUBROUTINE  EDIT! JERRI 


C  EDIT  TE"B  ANO  SALINITY  IN  T€ CRT S  YTO  RTTaO tF  • 
C 

C  JAN  28  1975  N  FOE ON OFF 
C 


EQUIVALENCE  IPDIFF.AII _ 

IERR  *  0 

0VP8AX  «  -.1?  - - - 

00  10  H«N1,N? 


PIN)  -  SIN) 

10  THIN)  •  TIN1  -  ' 

15  DVN  -  0.0 

PN  -  0.0  “ 

THN  •  0.0 


XNOP  -  NOP 
00  20  N-N1,N2 

0V»  «  OVH  *  DYffff - 

P"  •  PN  ♦  PIN) 

THM  -  THH  ♦  TWfWI 
PTIN)  -  PIN) 


■  I  n 

20  CONTINUE 

DVN  -  OVN/XNOR  - - - 

P N  -  PN/ XNOP 

THH  -  TH87XNUR  - 

CALL  LSET 


K 

00  60  N-N1,N2 

TFl  ABS I PTT  N)-PWl-AT*7TT4ff»7ff*7fCr 
30  CORR  •  POLY10VIN) ,DVN,CP,N,OiO,6000.0) 

KERR  -I  - 

IFI ISSWI31 131,32,32 


i 

M 

* 


40451 

40452 

40453 

40454 


111  i‘ri-1 1.|| 


51  *  0ATAXIN^r,2) 

52  •  OATAKIN,^! 

53  •  OATAXIH*!^?! 
TH1  •  THFTATOTTO-’ 


f 


5‘tJ, • 


99 


40475 

40456 

TH?  » 
TH3  - 

THETA(DVfN) *  DA  TAX  fH v 1 ) *s?VRE ) 
THETA(07(4*ll*0ATAX(4*l*l),S3tPF) 

* 

40450 

071  • 

‘  DVAIPEi  fHlt'Sl) - 

40470 

072  • 

OVA (PE  *TH2»S2) 

•  4 

40475 

0V3  • 

07  A  (  PF  *TH3»53 1 

- 

40482 

0VP1 

-  (0Vl-07?)70ELP 

404^0“ 

0VP2 

-  ( 0V2-DV3 ) /DEL  P 

40500 

WRITE(KLISTt3100)D7(N),P(  N)  ,CORR«  OFLTAVOVPl  ,07P2 

40400 

37 

P(") 

*  C  0**K 

40700 

40 

!F(AB$(TT(N)-THN)-A3*Z2) 55  »50 v50 

40800 

50 

CORR 

-  p  HL Y (  0 V  (  N  )  ,0V M ,  CS  V.O  .OVSCFOO.  0)  ♦  THH 

40900 

KERR 

-  1 

41000 

IFf ISSB*3) )51„52*52  - 

'  " 

! 

41050 

51 

OFLTA 

-  CORR^TH(B) 

' 

"WT  0T1 

"SI"  * 

41052 

S?  - 

0ATAX(*,7 J 

4T05T 

41054 

41055 

41056 
TTWT 

41066 

41074 

41082 

4T09O 

41100 


51  -  OATAYf *♦!,?)  " 

TH1  -  THET4f0Vf  N-l ) * DATAX ( *-l *1 )*S1 »PP I 

TH?  -  THtT»TtmwiT01TAYmiTr^57TFFl - 

TH3  -  THETA(0V<M*l).DATAX(fUl*l)*S3»PF) 

"trrr  -  o  vat  p  f  ,;tht  vttt - 

0V?  -  OVA(*>F,TH2,S2) 

D73  «  0VA(PF,TH3,ST)  - — —  -  - 

DVP1  -  <OV1-OV?)/OFLP 

DVP?  -  (  077 -07 T)  TUELP - - 


4t?oo 


-*rr 

55 


HRTTE(KLIST,3100)D7(H)fTH<H)*CBRR.DELTAvD7Pl,D7P2 

—  THtWI  '  ■  COPP - 

IF(KERR.FQ.O.AND.IERR.EQ.O) THEN 

JFRR--?  -  - -  - . -  --  - 

FNDIF 

TERR  «  IFRR  ♦  l  - 

60  CONTINUE 

mirnnr  tto  vtovto' 

70  00  75  H«W1,N2 
Sf*)  .  PCH) 

TCI  -  TH(  *) 

75  P(")  «  07(5) 

JERR  •  JERR  ♦  1 
40  RETURN 


!l 


4?  TO  O 
42200 
42750 
42100 
47400 
42500 
47600 
42700 
47000 
42900 
47000 
43100 


1100 


F0RPAT(F7.l,F9.3,*  REPLACED  BY*  ,.F9.'3,# 
♦  *  5R.  70L.  CRTTOTENTST  A  B 0 70  -  *VE9.T.» 

ENO 


CHANCE  I5t*,F9.3, 
BEUW  -  ^iT9.7) 


♦  #%***♦*♦♦♦♦*♦♦♦♦*♦*♦ ******** 

- ^ji|Pn\JT INF  L5FT - 

******************************** 


LEAST  SQUARES  REGRESSION  SUBROUTINE  FOR  POTEN. 
PAR  6  1976  N  FOFONOFF 

INCLUDE  ,CO*,POTEN.FOR» 

1  NA  -  N*(N*1)7? 

L  -  1 

no  10  I » l *  NA 

io  coiir  ■  o.o 

00  12  I-l,N 
CPU)  •  0.0 
1?  CTI  I  )  -  0.0 
15  OO  20  T-l,N 


100 


**100 

BP  til  -  0.0 

***00 

8TII1  -  0.0 

**5oo 

20 

Tontinuf 

**6oo 

00  8  I »P1  *  *2 

**700 

x  -  ovm  -  ov* 

**800 

**900 

00  too  J»1*N 

TFf J-l)90,90,9* 

**000 

00 

BIJ1  -  1.0 

4^  YOTT 

XU  TU  TUO 

**200 

0* 

8  I  J)  -  X*MJ-«1> 

*1100 

100 

CONTJNUF 

***00 

***00 

JK  -  0 

x  -  pn  n  -  p*  — 

**600 

XT  -  TTUJ  -  T*T* 

**700 

on  8  J-1,N 

**800 

8P(J)  -  8P(JI  ♦  8 1 J ) *X 

**900 

8T(  J)  -  BT(  J)  4.  8f  Jf#*T 

*6000 

*6100 

105 

IF  f  L—1 1 105*  105  *18 

00  7  K  •  J  *  M  ’  '  - 

*6  200 

JK  -  JK  4  1 

*6  iOO 

7 

C0(  JK)  *  CQ  f  JK1*  ♦  8(  Jl)  ♦  6<KI 

*6*00 

8 

CONTINUE 

*6*00 

IFa-ni73,171vt7* 

*6600 

*6700 

171 

17* 

CALL  5BINVfC0,N.B,BA,PR*IFA!U 

00  200  P-1  ,N  - 

*6800 

SP  -  0.0 

*6^00 

5T  .  0.n 

*7000 

JP  -  *-N 

*7100 

00  17*5  J-l.N 

*7?00 

*7100 

17*0 

PF  IJ-P 117*0, 17*0* 17*2 

jw  -  jp  4  n  j  j  4  r  — 

*7*00 

CO  TO  174* 

*7*00 

17*7 

JP"  -"  JW  *V  1 

*7600 

17** 

SP  -  SP  4  C0IJP)*8P(J) 

*7700 

17** 

IT  •  5T  4  C0CJP7*FTTJ1 

*7800 

*7000 

200 

CPIPI  -  CPIMI  4  SP 

CT(P)  «  CTIP*  4  ST 

*8000 

c  copputf  rfsidual* 

B  |  H'i1! 

T7T 

Rr  -  0.0 

*8?00 

RT  -  0.0 

*moo 

no  18*  I«Pl,P2  - 

*8*00 

F P  «  0.0 

*8*00 

FT  -  0.0  - - 

*8600 

x  -  ovm  -  dvp 

*8800 

NJ  -  N—  J 4 1 

*8000 

FP  .  pp*x  4  CPTWII  -  -  ”• 

*0000 

180 

FT  •  FT*X  4  CTWJl 

*0100 

SP  -  Pm  -  FP 

*0200 

PT  f I )  •  SP 

M  LI  Wu 

5T  -  THI1I  -  FI 

*0*00 

TTU)  -  ST 

*0*00 

TFa-KSMMF*,181Yl87 - 

*0600 

181 

RP  «  <  SP-PP1 **2  4  *p 

*0700 

RT  -  IST-THP)»*Z  4  FT 

*0800 

18* 

CONTINUF 

*0000 

C  *  L  ♦  1 

*0000 

IFIL-KSWU*,1*,19* 

*0100 

10* 

XN  -  NOP  -  W 

*0200 

*0100 

21  -  SORTIRPAXH) 

r?  -  soRTfrmon 

^  50500 


50700 

“5000(1 

50000 

51000 

51100 


7lfUU 

51300 

“51500 

51500 


IF ! TSSWI 6 1  1 300*350»  350 

300  no  no  i»i«n 


8HI  ■  A85ICF!  1 1 1 /( Zl«$ORT( AB5 (C0< 1 1 1 1 1 1 
9Fm  -  AH5nmr))Tr77^30irriw5<cocnvu) 

310  CONTINUE 

HR  ITT  I  Kt  T 5 T , 3 1  TO ) PT,N , 50P",  T  O  m  ,K -1 ,51 
MRITFlKLI5T,31lOMBFf  K»,K«l,N) 


3100  F0f»PAT(F6.0,t4,I3,8F8.31 

3110  FHR5AT  f  ITJIVBFB  131 
FNO 


102 


Appendix  C. 

Program  Listings  for  PEPLT 


PE PLT/PEPLS! 
K  PR  ISW 


103 

SHORT  OOCUHENT  AT  ION 

JSW  KLIST  DESCRf PTI ON 


jcall  t*ble  subroutine*  ltst, 
plot ♦  Output  in  hap  eornat 
plot  is* 

AI*V«PL<NXll»A?*VPBL<NX?l*A3*VRPLIN)r3MA4PCtIREC*’l» 
B1PVRBL  (NY1 )  ♦B2*VR  BL  ( NV?  I  ♦B3*VRBL  C  WT3  )»84*CtIREC*>?) 

-  -  CHANGE  PARAMETERS  FOR  PLOT 

INITIALIZE  PAR-AHETERS  TOR  PLOT 

-  CHANCE  DATA  SELECTION  VARIABLES 
CHANCE  PLOT  PATRAHETERS 

-  call  avrgs  Subroutine*  horizontal 

AVERAGES.  FOR  fYET AILED  DOCUHENTA- 


950 

1000 

1100 

r?oo 


AVRGS  BRANCH. 

SET  ISSN  (SMITCHI  ARRAY 
RESTART  MAIN  PROGRAH 


TOO  RERlT/PfPLS:  BRANCH  3 — PARAMETERS — SHORT  DOCUMENTATION 
300  K  Rff  3  [SW3  DESCRIPTION 

5W  i  '  o  '  prTnT” Out- p a r are t er s  on  klistt  store  common  to 

soo  FILE  KPLCM.  RETURN  TO  PEPLS. 

600  1  INPUT  VARIABLE  SELECTORS  NX l  TO  NZ  3 

700  2  ENTER  A1  TO  kb 

TIOO  2  2  ENTER  Bl  TO 

1200  3  2  ENTER  Cl  TO  Cfe 


105 

100  PFPLf/AVRGS?  SHORT  DOCUME NT* T TON 


200 

KBR 

ISW 

JSW 

KLTST 

DESCRIPTION 

__j 

250 

4 

0 

- 

- 

SHORT  DOCUMENTATION 

300 

4 

1 

• 

0 

READ  FROM  DATA  FILES  VARIABLES  IN  COLUMNS 

400 

JSW  TO  KLIST 

• 

500 

2 

0 

• 

ZERO  COLUMNS  JSW  TO  KLTST 

j 

600 

3 

1 

- 

INITIALIZE  AND  INPUT  PARAMETERS 

700 

0 

- 

INPUT  PARAMETERS — NO  INITIALIZATION 

-  _ u 

HLfMifl 

4~ 

I 

T 

n 

000 

5 

• 

- 

AOO  COLUMN  JSW  VERTICALLY  FROM  THE  TOP 

1000 

6 

- 

LU 

PR  TNT  OUT  DATA  ARRAY  ON  UNIT  KLIST 

1100 

7 

0 

— 

CALL  NCAR  PLOT  PACKAGE  TO  PLOT  ONE  FRAME. 

1200 

OEFAULT  IS  COLUMN  JSW  AGAINST  PRESSURE. 

1300 

GENERAL  PLOTS* 

_ j 

X»R1*C(  I. JSWI»B24CII*NX?)»B3*PR 

1500 

Y-ftl*pp4A?*CU,NYl  >*A3*C(  l,NY2» 

1600 

MULTIPLE  PLOTS  ON  ONE  FRAME  ALLOWED 

1700 

KBR 

ISW 

JSW 

KLIST 

1800 

8 

- 

- 

COMPUTF  DYNAMIC  HEIGHT  ANO  POTENTIAL  ENERGY: 

1900 

ASSUMES  DVI  IN  COLUMN  1  (NVlll-lB)  ANO  OVE 

■RlllfliH 

Rl|  I,  M  M lil'l 

2100 

9 

0 

M 

INTEGRATE  COLUMNS  JSW  TO  KLIST  AS  A  FUNCTION 

2200 

OF  PRESSURE 

2300 

10 

0 

0 

SUBTRACT  RFFERNCE  LEVEL  VALUE  CIJREF.t)  FROM 

2400 

COtUMNS  JSW  TO  KLIST 

2500 

11 

1 

- 

INPUT  JC1 * CR 1  TO  JC4.CR4 

- 7600 

~  IT 

— r t 

- 

2700 

Cl  I» JC1J-CR14CII, JC1I*CR2*CI IVJC2I* 

2800 

CR34CI I,JC3)*CR4*C< JREFWJC4) 

a 

2900 

12 

- 

- 

RETURN  TO  PEPLS 

3000 

13 

— 

- 

INPUT  COLUMN  «»S  AND  CONSTANTS  TO  PERFORM 

3100 

THE  FOLLOWING  COLUMN  MULT IPL ICATION: 

• 

Cl  IREC»  1 1  -CUN'l*CI  1  RtC»T  >VCONZ*CI  IRVC*  JI  * 

3300 

C0N3  *C  I  IREC»  KJ  .  INPUT  ORDER?  I,J,K,'C0N1 

3400 

C0N2  * CONI'?  AN” INDEX  II, J, OR  KT  OE  “V ALOE  -1 

3500 

PREVENTS  THE  INCLUSION  OF  THE  ASSOCIATED 

1600 

ANO  FOLLOWING  COLON  NT  YT. 

4 

3700 

K8R 

ISM 

JSW 

KLTST 

TWOTT 

nr 

— 

- 

3900 

ACCESS  TO  THIS  8RANCH  QUERIES  WHAT  HORIZONTAL 

4000 

LEVEL  »  IS  DEST1EO 

4100 

15 

- 

- 

NOT  USED 

4200 

16 

1 

- 

INPUT  X,J  " 

— 

4300 

0 

- 

Cl IRECt J)-C( IPEC, J)»RX.  I  SHOULD  FOLLOW 

5300 

^ 1 1 6 %  1  IHHf or^TFl.Y  IN  f  ^COTIOKr  J 

4500 

17 

- 

- 

ERROR  SUMMATION:  VERTICAL  INTEGRATION  WITH 

3600 

ITTELTA  PI*»2  AS  THE  INCREMENT 

4  700 

18 

- 

- 

INPUT  DELTA  P  INTO  C<IRECt5) 

TBOO 

19 

- 

- 

EXCHANGE  TWO  COLUMNS  OE  C 

-  ■* 

4900 

20 

- 

- 

CHANGE  A  S INGLE  ELEMENT  OF  C 

21 

— 

— 

5100 

X-BAR  IN  CUREC»41  ANO  X*X-B AR  IN  CfIRECtSI. 

5700 

RESULT  IS  STORED  TN  CITRECtll. 

6300 

22 

1 

1 

CALCULATE  DYNAMIC  HEIGHT  AT  A  GIVEN  LEVEL 

5400 

RELATIVE  TO  PATS  SURE  CORRESPONDING  TU  JREF 

'  " 

5500 

AND  OUTPUT  IN  MAP  FORMATtALONO  WITH  VARIABLES 

— 3*00 — 

F^t)^  COLUMNS  ?  ANO  4  AT  TMAT  j  jlJtE  • 

5700  NVfll  HOST  BE  •  18  ANO  NVI2I-19. 


COOPFPIT.FOR  FILE:  OTMFNSI ON*COWffOH  TND  FOITIVIOKCF  PO*  PFPLT 
01 SPLAY  PROGRAM.  N.BRAY 

- PARAMETER  KCH«943 - 

PARAMETER  JDIM*1'00 
CWA»ACTFR*8  DOIT 
CHARACTER*  12  GNAMp 
OIMFNSION  CSTr6*'4),VR(T5) 

0!  ME  NS  ION  KHDGM80)  *KBUF(46) 


404 

478 


COUPON 

COMMON  KIN 

BF GINNING  OF  STORED  COMMON 


600 

TOO 

800 

iroo 

1300 


COMMON  NY1  *NX?,iNX3,NYl,NY2*NY3*Nn,N72*N77 
COMMON  ISH*  JSW*W,NWVTN7TO 

COMP ON  YM I N  *  YM AY  ,'YN  T N , YW1Y 
COMMON  A1  *  A  2  «  A  3  *'A4  *  A8  «  A6 


1800 

T600 

1700 

T80TJ 

1900 


COMMON  C1,C?*C3*'C4*C8,C6 
COMMON  01*02  *D3  ♦'04*08*  D6 
COMMON  XI*  ZL TO  *i7LC0  » DA  Y *  XPL *  YPL 

COMMON  WT  '  - 

COMMON  NV  (  6 ) «NX(6)*AV«  BV  *  CV 


2100 
77C0 
2300 
7400 
2800 
‘  7500 
2633 
7666 
2683 
7700 
2800 
~  2900 
3000 
7T00" 

3200 

7700 

3400 


3600 
— T6T6 
3632 
— 3670 
4100 
— 4700- 
4300 
-  7400 
4428 
“  4670 
4478 


COMMON  CRl,CR2wCR3,CR4 

COMMON  JMAX, JREF 

COMMON  NX4,NV4„N74 

COMMON  IV1,TV2,TX7,  JBOF,  JHDX, TOO  _ 

COMMON  JSHPt  6) „0AY1*DAY2*PMIN,PMAX 

~roMMnN  yfmn.yfpxtxnmnvxnm' 

COMMON  C(  100,61 
COMMON  IS SHI 161 

COMMON  PL  A  BH  101 , XL  ABL  (  10) vYLABL  1 110  1 

COMMON  LTYPF,MHOR, I CON , I SHP • I C AST ,X0 AY  * 1PR ,LPR 

"COMMON  XL  AT*  XLONO*  W6T*XLTO*TL'GO - 

COMMON  LBBLf  3>*LBVfl3>  ,N5Cf  601  *NPRf  60UNSECTI0N 

COMMON  PF,TO,SO,OVO 

COMMON  PT"*THF,"SIF*TJTF -  -  -  -  - 

COMMON  PMf  THM, SM*’fJVM 


COMMON  CPm,71<.CT<8>,Z2*Fi,F2*F7 

END  OF  STORED  COMMON  . . . . . 

COMMON  XDATI 108V7) ,YDATt 100*7) 
COMMON  TOEtP,OP-  - - - - 


COMMON/CHARACTER/  CN1ME ( 200) *DOC f 10 1 

FOUIVALENCE  C A1 ,CST) * < VR ,PF ) 
FOUTVALTNCE  fXHDG,LTYPFirnniUF^XTYPE ) 
FOUIVALENCE  ( KTTX*  KPLC  M ) 


I 


C  PFPLT  PPOG  ***♦♦*♦**♦  SFPT  24  1<JT7 


t 


J 


TOC 

200 

- TTTtT 

400 
'  ~50C 
600 
roc 
800 
— tot r 
1000 
rioo 

1150 

tttjc 

1300 
"T500 
1500 
VS  00 
1700 
1800 
1900 
"7000' 


C  PFPC5  SUBROUTINE  ♦♦♦♦♦♦*♦ 
SUBROUTINE  PERES' 

c 
c 

C 
C 


SEPT  24  1977  *♦♦♦♦ 


RROOPA"  TO  PLOT  POTFN  VARIABLES. 
JUNE  27  1976  N  FOEONOFF 
VHXTF17TTTTN--N0V  I9>W 

INCLUDE  »CO«PEPLT.FDR# 


40 


120 

T2TTP 


WRITE(KTTX,40> 

FOlfNATtlH  ,  *!-NTTIALl7FTrCTAy  VARUBtE  SELECtlON  PARAMETERS 
IYFS  OR  NO)?M 

TF f NOYES ( KlN,RTTX) .FO • 1 VCO  TCT06 
WR!TE(KTTX,1?00»KBR, ISW*JSW*KtlST,KTP*KOUT,KlN 

'  -FTTgmnT  ♦  RPEPUrtKBR  ,  IS*'. TSVt KL"TTT« K  rPirmmiOTfT'A, - 

♦  6X,3I3, 15, 13*14,13) 

KUST  -  6 
KOUT  -  8 
KTp  .  n 

REA0(KIN,*)K8R',  ISW,JSW,KLIST »KTP, KOUT ,K TN 

- rr<  K  gg-.T  T77  TKIWS7 - 

2100  rE(KBRU20,B0OV130 

2700  130  GO  TO  (200, 20,30, 400, 500,600,T00TKBR 

2300  C  *♦♦*♦♦*♦*♦♦♦*♦♦  +  ****♦♦♦**♦♦♦♦*♦***♦*♦****♦♦♦*■**♦♦**♦•♦*♦**•**♦♦**♦ 

7400  C  TNTTYALT7E  ~  " 

2500  C 


jnJO 

3400 
T50C  ~ 
3600 
"  3700 

3800 

Tuo 

a  FTN*~7TJ 

XBAX  -  100.0 

YNTN  -  0.0  ' 

Yf»AX  -  5000. 8 

A I  -  1.0  - 

B l  -  1.0 

LI  - 

I  .o 

4000 

01  - 

1.0 

4100 

07  * 

1.0 

4200 

07  - 

1.0 

4100 

04  - 

1.0 

4400 

05  « 

1.0 

4500 

“06  * 

T"#  0 

4600 

NX  1 

12 

4  TOO 

NX? 

0 

4800 

NX3 

0 

- 4900 

NY1 

19  ~  ‘ 

* 

5000 

NY? 

0 

5100 

NY3 

0 

5200 

N71 

25 

5300 

N77 

0 

-- 

5400 

N73 

0 

5500 

NV  - 

7  -----  -  ■  . 

-  - 

5600 

8W  - 

0 

- 8-T00 

■  ~KTP 

■"“Tt - 

5800 

CALL 

OAT  A I  -l  »iIEOE  ) 

5900 

GO  TO 

120  “  ----- 

6000 

C> 

***** 

****** 

5100 

700 

CALL 

table  —  •  — 

6200 

GO  TO 

120 

6  iuu 

~C 

Ttff 

•  R****1 

RTT  7JR ANCH  7--C5ANGF  OAT*  lELFCTTTJN  VARIABLES'  ****** 

6400 

20 

CALL 

OATAfO,lFOFI 

6500 

GO  TO 

120 

- 

6600 

C 

6700 

30 

OPEN  tUNI  T  «50  ,.Tf»5€«~*yE  PLSl.OtTTT' mr  P  F-*OLO»  ,9EITOTTNLTT 

-  -  — 

~6iW 

6900 


7100 

7200 

7300 

7400 

7500 


ri.iiM 


7700 

“7100 

7900 

0000 

8100 


8300 

8400 

8500 

8600 


IliU 


itlTi: 


00  3350  N-1,200 

REA0(50»  3325»lE'ND-3340 )  <D0C (  I  l-Vt ■1*9) 


W 

3350  CONTINUE 
3325  FORMAT  (  9A  8  I 
3330  FORMAT! IH  ,9«8) 

3340  CLOSE  (UNIT. 90) 

300  WRITE ( KTTX, 3000) 


J  t  i  t.llijLl 


.1*1 


KX  -  4 

R E  Af» f  K f  N,  ♦  ) KtR3,r3W3VKXVMV,MW'  ~  ~  ' 

IF(KBR3>30, 350*31 

31  IF  (ISW3-1) 32*32*74 

32  WRITE(KLIST,32001NXI,NX2,NX3,«fYl»NY2*NY3,NZl,4l72fW73 


iiii  i.r  mi: 


:  k>/i ;  k'A  i 


8700 

34 

Booo 

3400 

8900 

9000 

35 

9100 

37 

9200 

38 

9300 

fS*rYiTMi 

t 

9500 

350 

9700 

9800 

- 9900 

10000 


10200 

T030TJ 

10400 

10600 

10600 


10800 

tuooit 

11000 

TTTOO 

11200 


IIJUW 

11400 

~TT50(T 

11600 


IF  (  TSW3-1 )  34»33«'30 

READ!  K  !  N*  * )  NXWNX  2  *  N  X  3  *  NY  1 «  N  Y2  *’NY  3  *  Nf  1 4 NZ  2  «  Nf  3 
CO  TO  300  _ _ _ 

WRITE(KL IST,3400)KBR3,(CST( JCVK8R3t« JOT,  6) 


V1IC9 

TF(  ISW3 )  30,  38,,35 

READ*  X  I N , ♦ ) ( CS T ( JC ,K  8 R  3  )  ,  JO  1  *6 1 

GO  TO  300  _ _ 

KRP3  -  KBR3  ♦  1 
IF ( KBR3—KX ) 34*34, 350 


WR ITE ( 10  *  1 )KPLCN 
CO  TO  120  - 


C  ♦♦♦♦♦• ♦♦♦A V  E RAGES  >4 
400  CALL  AVRG5 


LTiMJ 


rimmi 


11800 


k  kivv 
12000 
I?  TOO 
12200 


C5  ♦*♦♦♦♦  SET  IS5W  SWITCHES 

5000  FORMAT ( 2 ( 1H  , X, 16T  4, / )  ,  *  ENTE*  K,ISSW!K>M 
REAOtRTN,  ♦HX^ISSWfTn  ,7(-r7T6I  ~ 

GO  TO  120 


MAIN 

600  RETURN 

C  *♦♦*♦*♦***  EXIT  PROGRAM  ♦♦♦**♦♦♦»♦*♦  « *  *  *  *  WVWmWW 
700  WRITE ( KTTX*7000 ) 

TFIN0YF5  ( KTNVNTTX) . NE .11 GO  TO  120 - 

STOP 


f wwu  runn» i t in  rn 

O*  PPPLS*  SHORT  DOCUMENTATION — BRANCH  0  ♦♦♦*♦**♦♦♦♦#*♦ 
BOU  OPEN!  UNI  T^SO  t.MAME p  EELS'.  OPC «  ,'TTP  E  -  * BCTT*  ,  R 

DO  850  N-1,200 

READ(5P*B?5»EN'P»BT7TnJ0C(  I  TTT^TTTf - - 

WRITE(KTTX,830I(D0C( I),I-T,9> 


CLOSE (UN IT -50) 

furmitptab  ) 

FORMAT! IH  ,9*8) 
GO  TO  T70 
ENO 


12400 


X  TVTnTT  SUBP50 G  * ♦  *  **99**  * MV ***M MMM 
SUBROUTINE  AVRGS 


400 
'^50  (T 
600 


C 

C  TOR  WJRT70NTAU  IVEPACFS  CDPPUTELT3T1ND  PLOT.’ 
C  JUMP  27  1976  N  FOFONOFF 

c  ^rrTFRsrmv— uuvT'rBo;  nvbrat - 

OfMFNSION  0(6) 


kidif  vim 


1000 

1050 - 

1100 

— 1200  X 
1  TOO 


1  TVU  K, 

1450 

1500  X 


DIMENSION  XY*<41 
CHAR  ACTER*  12  0UTNA4F 
TNCLUOE  *C04PEPfcT.f OR* 

EQUIVALENCE  <D1VD) 


CHARACTER*^  1DSTN 


1600  10  GO T0 <100*200. 300*  400*  500.600. 700*800* 900*1 000. 1100 .1200.1 300. 

1700  *  l400,ISS0.1^00TT700,T8W."190r072r000.2100^2ZOO¥T500TTSir 

1800  C 


111111 22123 ll 


I  i.l»: 


:IT11ITTTTTTT1T1 


2000 

2700 

2300 

T40d 

2500 


100  CONTINUE 

DO  ID  l  K«1 *  NO  - - -  ~ 

RE AO <12. 101 l *END»101 2 lM.GNAHEt K  >  ,  NT 
T0I1  TrrRNATTI4,A17?F5.Tl 

GNANEfRI <  9tl?)-*.AVG* 


2625  GO  TP  1013 

7650  101?  CONTINUE  -  "  “ 

2675  NO  -  K-l 

?687  REWIND  1? 

2700  1015  CONTINUE 


i  ■ 

2900  IFOF  *  0 

3000  OPFN<  UN  T  T  -  RTP  ,  NA  «  F-G  N  A  PEINIT  THREAD  0  NLY  .TYPE  *  ffTLO1 

3100  ♦  *  UNFQRN A  TTEO * »PRR*168 ) 

3700  1 02  CALL  OATACl.TXDFI  - ~ -  - - 

3300  IF< IE0FI165, 105.105 

— y^OO - TD5 - WT  ~  -  T  .'O - 

3500  IF<ISSW<15)U10,'115,115 

1600  I ID  WT  -  WGT  - 

3700  115  00  160  I ■ JS W.KL1 ST 

— 3900 170  XT  -  wBLTNvrrrr  - - -  - 

4000  CdREC.I  J-C(t*FX.T>  ♦  0(1 )  *WT*fAV*XT-*T‘BY«’CV*XT)  * 

5TUD - f — 7RTXTNXTTTT) - 

4200  160  CONTINUE 

— 5300  DO  TO  102  -  -  - - 

4350  165  TF< TSW.EO. 22)00  TO  800 

-5155"  GO  TO  170  - - -  - 

4382  168  WRTTE<KTTX,4)*FRR0R  READING*  .GNANE  (N3T1- 


9TVV  1  (  V  v  1719  1  I  HUT 

4600  WRITE<lO*l)KPLC* 

"4700  130  GO  TO  1500  - - -  - 

4800  C  ♦♦♦♦*•*♦*♦♦  «2  ZERO  TA8LF  SET  PARAMETERS  *****9**** 

-4900  200  00  210  1-JSN.RtTST  - 

5000  00  210  J-1,100 


7  OLD*  *FORM> 


4200  IF(  ISW.EO. 2  7.  A<ND.LE  ILE.EQ.OIGO  TO  170 

5100  GO  TO  1400  - 

5400  C  ♦*♦♦•♦•♦♦♦♦♦  *3  SE?  PARAMETERS  ♦*♦♦♦♦♦♦♦»#♦♦*♦♦♦♦♦#♦♦**♦ 

5500  TOO  ITT  TSUI  370*  370  7310-  - - - 


in 


•5600 

5700 


5900 

6000 

6100 

6700 

6500 


6500 

6600 

6700 

6000 

6900 


u«m«i 


7100 

"7700 

7300 

7*00 

7500 


r  ouu 

7700 

7000 

7900 

8000 

8100 


8300 

T5U0 

8500 

8600 

8700 


*1 


9500 

9600 

9700 

9000 

9900 


10100 

T0700 

10300 

T0500 

10500 


IUOUU 

10700 

T0600 

10900 

TTOOO 

11100 


11300 
11500 
11500 
IT  600 


>10  -  l 
NVC11  -  51 


NVm  -  86 
NV 161  -  87 
NVC51  -  63 
NVC61  -  -I 
00  317  1-1,6 


JC1  -  l 
CR1  -  1.0 
JC 7  -  ? 
C»?  —1.0 
JC3  -  3 


JC5  -  * 
CR*  -  1.0 
AV-1. 
BV-0. 

cv-o. 


JEFF  -  50 

370  TF  I  JSW)  350, 375V775 

375  WRITECKTTX,3?tt01N0,1NV!Kl,K-lW6l,  JREfUJOAX 
RFAOC  KIN  ,  *  1  NO  ,'  f  NT  IX  T  ♦  K -1,61 ,  JNEF  ,  JWAR 
330  WRrTE(KTTX,3300)AV,BV,CV, (NXCTI^I-1,61 


fvmrAn  i;t(«yn  o  m  i 


3*0  HR  ITE  C  KTTX,  3*001  A1',A7,  A3, 81, 62*63 

R  F AO  <  RI N , ♦ 1 A  lyl7.1T,  8 1 ,87183 
WRITE  CRT TV, 3500 1 NX 1 , NX?,NVl,NY? 

RFAOmNi* )  NTT  /NX?  ,NYT,NY7 - - 

775  GO  TO  1500 


3300  FflRHATIlH  ,  «  AV^B  V  ,  CV  ,  NX  f  6  >  •  ,  f  *3F6. 3 , 61  3  V 

3500  F0RW1TCIH  .'PLOT  R8WA8FTFR5: - 11 -  A? - 

♦  B1  B7  B3«,/rl:6X,6F9.3> 

3500  FOROATCTH  ,f  NY  1  "NR7  NTT - NY7r,  7»5T5T 

C  ♦♦♦♦♦♦♦♦♦♦♦  •  *  AVER AGP  TABLE 


IT (CCJ, 6 11*05,615, 505 

505  TO  510  T-JSH,RLTST  - - 

910  ru*n  -  c  <  j,i  >/cc  j,m 

515  JNAR  -  J  -  T  - 

IFCISSWC?!  1*70  ,’150  0,1500 


■  l  v  •  »  u  ,  »▼  »  t  /  t  jit, 

575  F0R8ATC1H  ,5!T3,7R1) 

GO  TO  1500  - - 

C  ***************  • 6  A 00  COLUBN  JSH  ♦*♦♦♦♦♦**9* 

500  TO  510  J-7,  JWAY  - 

510  CCJ, JSH 1  «  CCJ-IVJSW)  ♦  C( J, JSH1 


*6  L  TST  TABLE 


HRfTEfKLIlT  ,60551  NO,  JREF,  JBARft  NV(K  1  .IC-T ,61 
HR  TTE  art  T3T  ,6056  ITT,  BV,C7t(NTF3  1,1-1,61 
00  610  J- l , J WAR 


HR  I  TEC  XL  I  ST, 60501 J,KP,(C4 J,K)VK-1,61 

610  CONTINUE  - - - 

GO  TO  1500 

6000  FORMAT ITH  VITR51  — - 


$ 


-VI 


't 

m 


□  rrroo 

11800 


i  nwu 
12000 

12200 

T73TJ0 

12400 


uitlt] 


12600 

T7T00 

12800 

17300 

13000 


l  3  L\J\J 

13200 

13300 

13400 

1T500 

13600 


13800 
13000 
14000 
'14300 
14200 
— 14300 
14400 
14300 

14600 

TT700 

14800 


i  i-»uu 
16000 
36100 
16200 
T57UO 
16400 


16600 

T6T00 

16800 

16000 

16000 


161DU 

16200 

16100 

16400 

16500 

16600 


•  V  '  W  V 

16800 

16300 

17000 

13100 

17200 


1  I  JUW 

17400 

17600 

17600 


6050  FORMAT ( I3,T5»7X»6FiO.  41 . 

6066  FORMAT < 1H  *314*618) 


ni'-jo  rutrin, 

C  ***********  # 7  PLOT  TABLE 
700  CONTINUE  ” 

JM1N-1 

’  NTURV-l  - - 

HRITE(KTTX,'7500) 


jiiiiLtUfUL'MiLi'niyi: 


IFfNCURV.CT.6)NCURV»6 
MR ITE  f  KTTX,701O)PLA8L 
IF  fNOYESCKIN.KTTX) .EQ.1 1THEN 

RFA0(K!N,7020)RXABt  —  . 

CALL  STRIPfPLABL) 


c*u  i  r 

MR  ITEIKTTX, 70301 
IF (NOYES (K  CN, KTTX I  . F 0 . -1 ) THEN 
MR  ITEIKTTX, •7040)X4fN,XM AX *YMIN*YM AX 
R  F  A  n  (  K I N  ,  vyXMf  Ifi  AT,  Y«  IN*  Yfmf 
CALL  AGSETE!  6KX/MI7I.  ,  XMIN) 

TaCI 

CALL  AGSETE ( 6HY/MIN. ,YMIN) 

CALL  AGSETE  <frHY/NAX.,¥HATO~' 

FN^IF 

MR  !TF  (KTTX,7060)XLA8X  . 

IF  f NOYES CKTN,KTTX).E0.1 ) THEN 

~  R  E  A 0  f KTN ,  70201  XLTBt - 

CALL  STR I P(XLABL) 

FNOIF 

MRITEfKTTX,  7060) YL  ABL 
TF  f  NOY F  S  t  K  T  N  , KTTX1  .TO TT  TTHFN 
RFA0IKIN,702O)YLABL 


LIAR 


FNOTF 

5TT  UP  PLOT 
CALL 
CALL 


LABFL  - - 

AGSETF(11HLABEL2NAHE.,1HT) 
ACSETF II 7HL TNF/NUMBTTR.  ,V85 ) 


PARAMETERS 


E7MXT  PLOT 


I  %  ff  Nk  % 


CALL  AGSETEI 4UR0M. ,2. ) 
CALL  AG3ETFr6HFRABF.  ,TT 
READ  DATA  INTO  PLOT  ARRAYS 
00  710  K-1,NCURV 

MRITE(KTTX,7070) JSM 


IW‘ 


TEITSSW(S) ,E0.— 11THFN 
HR  ITF  nCTTY,70B0) 

RE AO I K T N, 7090)1 CHARI K) 

FNOTF  -  - 

I  JH-0 


I JM» I  JM*1 

PR*FLtTA  TT  JPRTT^HFR  ,'NSEC1 
XOAT( I JH,*)-B1*C< J,JSH) 
YOATI  T  Jff«KT«a T*WT  r~t7+t 
710  CONTINUE 


♦  82*C 1 J,NX21  ♦  B3*PR 
T  JfWTl  )  *  S3*CT7,NT7t 


CALL  AGSETE (8HY/0RDFR. ,1* ) 

TL5F  - 

CALL  AGSETE  f  8HY/0R0E  R. ,0. ) 
FNOTF  ~  - 


3  ■  ■}  i  .  .  «r  .1.  m  '  I, 


113 


17800 

17900 

18000 

18100 

moo 

18100 

TBT600 

18500 


CALL  ANOT  AT  <  XL  AftL  *YL  AftL,0»0, 0,0. 1 

CALL  F7NXYIX0AT,YDAT, JDIN,NCURV,T  JN,PtABLl 

~mX—AC'CC TP  ( 15H5EC0nD'ASY?uSF* .  r?Y*76  1 - 

?Ff  O.GT.XVNfl  l.AN0.0.LT.XYNf2)lTHEN 
CALL  LTNFNCARI0.,XYWm,0 .,XYNT4 n 
FNO  IF 

TFUSSWI5).F0.-il)THFN 
00  720  1-1 , NCURV 


- 1 8  BOO 

18700 

18800 

18900 

19000 

19100 

720 

- CALL  "PniNTSUTOAT'i  l,rr.YDAT  11,71,1  JW,  ICNrRTITyTn - 

CONTINUE 

ENOIF  '  ~~  . 

CALL  AGSE  TF  f 6HX/N ON. « l , E36 I 

CALL  ACS ETF ( 6HX/NAX. , 1  .E36 1 

CALL  AGSETFf6HV/NIN. ,l.E361 

— 

- 

19300 

CALL  FRANE 

19506 

GO  TO  1500 

19500 

C  FORNATS 

19600 

7000 

FOR  NAT  f  1H  »,TNPUT  *  OF  CURVES  lN  THIS  PtOTfHAX  IS  61? 

19700 

►  AND  INOEX  OF  FIRST  POtNT*»l 

- 19800 

70TTT 

F09"8'A"T  E1H  *  ’CHANGE  PLOT  LABEL?  "DOmCS  FL  I5SVV/7AH 

7TTJT5T 

19900 

7020 

FORNAT  f 1 0A6 1 

2O0OO 

7030 

FORNATCIH  ,  *  USE  OFF AULT  AXIS  P  Aft  A NET  E  ITS  ?  *  1 

20100 

706  0 

FOPNATIlH  *  *  CURRENT  VALUES  OF  XNIN, XNAX ,7NIN, YNAXY * , 7, 

6E10.31 

20200 

7050 

FORNATIIH  , *  CHANGE  X-AXIS  LABEL?  OLD  LABEL 

,10161 

20300 

7060 

EORNAT ( 1H  ,'CHANCF  Y-AXIS  LABEL!  OLO  tATBEL  IS*N2V4H 

, 10A6 1 

20500 

7080 

FORNATflH  ,  ’INPUT  IDENTIFYING  CHARACTERS  1 

20600 

7090 

F0RNATIA1) 

20700 

C  ***********  *a  COMPUTE  OH  AND  PE  ♦**♦♦*♦♦******♦*** 

’28800 

800 

PPR  -  0.0 

¥■* 

20900 

OEL A  -  Cfl,lf-Cfi,2) 

p  i 

DF LB-'- "DEL A 

•  ■ 

21100 

00  820  J-1,JNAX 

21200 

805 

PR  -  JPRU,NPR  wNSFCTIDNl 

21300 

OELP  -  PR— PPR 

21  VO  0 

OFLA  -  Cf  J,l)-*Cf  J,77 

21500 

DHX  -  0, 5* f OEL A ♦DELB  1 P0ELP 

—J miuniinj J.iiii :  m  i  i  in.  — —— 

21700 

IFfJ-11815  *  818V815 

7IOT0 

810 

C  11,11  -  DHX 

21900 

Cfl, 2)  -  PEX 

72900 

GO  TO  817 

22100 

815 

CfJ,il  -  CfJ-1,1!  ♦  OHX 

C  t  J  »  2 1  •  CfJ-lV2l  ♦  PER 

22300 

81  7 

DELB  -  DEL A 

27500 

820 

PPR  -  PR  . .  '  '  .  '  '  - 

~  . .  ■ 

22500 

rFfISW.EQ.221 THEN 

77575 

RLTST-7  '  -  '  ' 

22550 

GO  TO  1000 

* 

77575 

FNDTF 

22600 

TFI  TSSWf 21  1825,1500  ,1500 

% 

77700 

825 

WRTTFI5,5?51X5ir?TSir,  JSWtKLTST 

- - - - — ~~ 

■  ■ 

22800 

GO  TO  1500 

72900 

— 

23000 

900 

DO  950  I- JSWyXLTST 

75100 

PPR  ■  0,0 

23200 

CPR  -  Cfl, 11 

77500 

DO  960  J-I,J*»AX  - - - 

- — -  — 

23*00  910  PR  -  JPRf J,NPR,HSFCTI0N1 

77500  TFT  TSW.FfT.T77  00  TO  91 7 


1 

114 

. “~a 

27600 

23700 

— mm - o- 

pex  -  o.5*(c( j,ri*T:PRT*n»p~PFii 

GO  TO  918 

23900 

7407a 

24100 

-24700- 

24300 


£ 

24500 

24600 

24700 

74600 

24900 


25100 

75200 

25300 

27400 

25450 


25700 

25800 

25900 

76000 

26100 

-26700 

26300 

76400 

26500 

76600 

26700 


26900 

77000 

27100 

77700 

27300 


C  f  TVV 

27500 

77600 

27700 

27000 

27900 


C  ouuu 
28100 
70700 
28300 
70400 
28500 


t 


28700 

70000 

28900 

29000 

29100 


29300 
5*  29400 

$  29500 

^  79600 


910  IFI J-l 1930,920,730 

720  Cfl,Tl  •  PFX  "  - 

GO  to  940 

737  COO  -  TT  JVn  - - - —  — 

c(j*d  *  ccj-*wii  ♦  pex 


950  CONTINUE 

TF  ( I55VT(  2  J 1 975,1500,1500  - - 

975  WRITE(4,425IKRR, ISW, JSW.KLTST 

GO  TO  1500  -  - 

C  ***********  «10  SUBTRACT  REFERENCE  VALUF  *♦♦**♦*♦♦* 


iVlUV  UU  IU  l"  JJR15U 

C»EF  -  CfJREF,!) 

00  1040  J*l,JMAX  . . 

1010  C  f  J  *  I)  •  CREF  -  C(  J*  II 

1040  CONTIWF  - - -  - - - 

TF ( TSW.FO.22l  GO  TO  1425 


IF  USSWI  2  )  110  75,1500,  1500 

1075  WR I  TF  f4,4?5)KRR^I5W,  J  S  W  , XL  TSTT "  - 

GO  to  1500 

C  ************  *ii  A00  COLUWNS  ♦♦***9449444*47****# . . 

1100  IF f JSMI 11 10,1120,1110 

— mo 7RTTFTKTTT ,  ITT  5 1 T CT,  CR 1 ,  JT7,Cr2  ,JC  TVER  3  *JC  4  ,C  *  4 - 

RE  AO (KIN,*) JCI,CR1,4C2,CR2, JC34CR3* JC4VCR4 
GO  to  1500  -  - 

1115  FORMAT ( 1H  ,*  JCl,CRl,  JC?*CR2*  JC7*CR3,  J<C4,CR4»*A,4*T4,E12.4U 
1170  00  1125  J-l,  JMAT  -  -  -  - 

1125  C(J,JC1)  -  CR1*C( J,JCl  l*CR2#C(i,JC2»*CR3*C( J,JC3) 


IF (  TSSWf  2  I  11150,1500,1500 

1150  WRTTE(4,425IKB*,1W,  J5H,Ktn7 - 

WRITE (4, 11151 JCI,CR1, JC2,CR2*JC3*CR3, JC4,CR4 
GO  TO  1500  -  - 

C  ********* 


C  *******************  ll3  MULTIPLY  UP  T03  COLUMNS  ♦*#*♦♦**♦*♦*♦*•** 

non  i— i  -  -  -  -  - . -  - 

j--i 

CONl-l 


C0N3-1 

WRITE  (KTTX,  13101  -  -  - 

1310  FORMAT! 1H  , • INPUT  COLUMN  NUMBERS  UP  TO  3  VALUES*  AND  CORRESPOND! 

X  MULTIPLICATIVE  C0NTTANTST1 - - - 

REAO(KTN,*l  I  ,‘J **K , CONI ,CON? , COM3 

- rrry  1390  trec»i  *.jmax - 

IFfI.LF.0IG0  TO  1500 

A-ct  rRFC,n#cmrr  -  - 

IFfJ.LF.OIGO  TO  1380 

8*Cf  I7FC, J)*C0N2  - 

TEfK.LF.OIGO  TO  1381 


luetic 

GO  TO  1385 

1780  CONTINUE  — - 

8*1. 

13WT  CONTINUE  -  - 


•a 

♦ 

>1 


79700 

29800 


Mi/Mi] 


30000 

T0T00 

30200 

tutou 

30400 


30600 


30800 

~3UB75 

30850 


JUBD  C 

30875 
'  “  30887 
30900 
” "3T0UO 
31050 


3  1VOO 

31082 

7TTO0 

31206 

- TT2TZ 

31218 


31236 

31747 

31248 

31254 

31260 


JICUO 

31272 

"31771 

31284 

TT79U 

32183 


32300 
32400 
3  2500 

- 32550 

32552 


jtjii 

32577 

776U0 

32800 

T7900 

32912 


32936 


1 

Ml 


32958 

3Z979 

33000 


33150 

33173 

33200 

T4TOO 


1389 


C0-1. 

C!  IREC,I  l-A*1»*C0 


lLIiMKiI 


L'Hl; 


1395 


ITT  IT 


IF!  I SSW!2)  11395,1500,1500 
MRTTET  4,4  25) KNIT,  T 3W ,J 5 W VXL  TST 
MRITF  ( 4 , 1 3 1 1 »  I’,  J  ,  K  ,  C ON  1  ,C0N2  ,-C0N3 

GO  TO  1500 


1400 


1402 


142 

3475 


KT0  -  60 

JRFC  -  1  -  - - 

WRITE!KTTX,1404)OUTNAME 

IFTNOYFS I  KTV,  KTTX7.F0  .  I TTHFN - 

READ1KTN, 140?) OUT NAME 


t  NIJ  1  r 

FORMA  T !  A  l  2 ) 

MR I TF  f  K  T  TX , 1 403 1  JR  PCI  , JRFC7 - 

REA0CKIN,4) JRPC1, JREC2 


OPEN! UNIT  -KTO*NAHE«OUTNAME  »TYPP*#  NEW*  1 


KLIST-6 

GO  TO  100  - - - 

IF!  IREC.NE. JRECI1G0  TO  102 
LFTLF  -  0 

I0STN(lS2)  -GNAME’I  NST  1(1:2) 


iiiAtLikaimLUACki 


1430 


1 4917 


1421 


1422 


1403 

1404 


TE!  ISW.NE.22)  60  TO  1430 
00  145  KREC-JREC1, J9FCT 
VR1  «  C!KREC,7S«l*l,E-3 
VR2-C(KRPC,3) 
VR3-C(KREC,4) 


VR5-C(KPEC,6) 

KP-JpR<KREC,NPR,NSfCTIOST - -  - 

MRITE(KT0,1421)  I  DSTN  »KP,XL AT  t'XL’ONG,  VR1  4VR2 , VR3,VR4, VR5 

CUNTTNUF  - — - - 

TE ! NST«EQ, NO  ILF  ILF -1 


•  u 

VR1  -  VRBLfNVmi 

VR2  -  VRBLTWTZTl  - - - - - 

VR 3  -  VRBL!NY!3fl 

CONTINUE  -  - - - - - 

KP- JPR ! KREC ,NPR ,  NSEC  T ION) 


iu  j  *  fi  i  irnrnu 

I0STN!3:5)-GNAME!NSTM6*8I 

WRTTFIKT  0,14  71 TTD3TN  ,5F,XL  AT,TLONG  ,Y4r,rR7,YR9  ,5R5VVR5  — 

iM  utw(Ev 
KLIST-6 


bU  UJ  £IIU 

ENOIF 

TF fTSW.FO , 22, ANP«LFTT  E  »  F0« 1 ) HRTTE I R  TO, TVZZl - 

FORMAT!/) 

CtOTFTUNlT-irrOT  - —  - 

TF  C  LF ILF ) 200,102, 200 


»  i  1  4n  9  nt  If  uw  *  r 

FORMAT! 1H  ,»LFVfL  NUMBERS  ARE  •,2131 

FOR51TI1N  ,*TNPUT  NF)T  OUTPUT  FILE  NANETTE  OR  NOIT  "STITT 
BRANCH  0— SHORT  OOCUHENTAT TON  ♦♦♦♦♦♦♦♦♦♦ 


1550  GU  TO  1500 


116 


-  —  l  jij.i  Ijji  W30th  ‘  J 


347( HJ  C  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

34300  1400  WRITE!KTTX,1505)KBR, ISW, JSW,KL1ST 


wwo - rooo - fprmatiih  ^RVRTrstrB-RisnjswvinsT  *,/,'5x,m,i53 — 

34400  KLIST  -  6 

74500  RFAD1KTN,  ♦)K9R*T’3TJ*  J5W,KLT3T -  *  ~ 

34700  IFIKBP.E0.0)  CO  TO  ?300 

34806  FF (K8R.GT • 771GO  TO  1500  - - 

34900  GO  TO  10 

J56UU  C  ♦♦♦♦♦♦♦♦♦♦♦♦♦»»»5  lT^~TAKr~TlTRFCt  TTTT**¥V¥WFT*¥T¥W 

35100  16 Off  FF( JSW)I610,16?0, 1610 

35700  1610  WRITFfKTTy* 16113  ~  - - 

35300  RE  A0It<IN,4>X,  J 

15400  GO  TO  1500  '  - - -  - 

35500  16 70  00  1630  IREC-l,JMAX 


- 34600 - C  1 TRFC ?  J  T *~A HTCTTR F C~,"J 1  T»»7 

357QQ  1630  CONTINUE 

~3510O  FFfTSSW!  2)  1 1650*T500»  1500  - 

35900  1650  WRITE!4,4?5))C4R,!5W,JSW, KLIST 

3 &mnr  •  '  '  WR  TTF  !  4  , 1651 6*  «1  - 

36100  1651  FORMAT  ! 1H  ,E17.6, 7X, l 7 I 


"16700 

36300 

16400 

36500 

16600 

36700 


- 00101500 

1611  FORMAT  !  1 H  ,  *  INPUT  EXPONENT ,  COLUMN  *  1 
t  BRANCH  17— SUWWITTIJN  OF  ERRORS  0¥ER  P 

1700  GO  TO  900 

1765  FORMAT!  1H  ,2!T7,?YTF  - 

C  *****+****+********  •#!«  —  DELP  INTO  C!IREC*5)  ♦*♦**♦***♦♦**♦*4** 

— rttyty ppp«q.o - 


♦♦♦♦♦♦♦♦♦♦♦♦♦ 


- 36800 

36900 
— 37000 

37100 

1T700 

37300 


00  1810  J«l,JMAX 

P R  «  JPR  !  J  ,  NPR  *iN5FCT  I0N1 
OELP-PR-PPR 

PPR.PTT 

C!  J,5)«0ELP 


17400  mo  CONTTNU? - 

37500  IF!  ISSm?)U875, 1500,  1500 

37600  IR?5  WRTTE(4,4?5)KW,T5H»  JS¥,XI:T3T  - 

37700  GO  TO  1500 

“37800“  C  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  »1 9—  EXCHANGE  COLUMNS  ♦♦*♦♦#♦♦*♦♦*♦♦**♦* 
37900  1900  WRITE!KTTX,1910) 

18000  T9TO  E OP *67 T 1H  T INPUT  COLUMN"  NUWBUTS"  TO  OF  "EXCHANGED*  1 - 

38100  READ<KIN,*H  ,  J 

39700  00  1970  IREC«1*JWAX  -  “  "  - 

38300  C  f REC  -  C(IRECVI) 

31400  CUREC*!)  -  C !  f  RFC  *  J I  - 

38500  CtIRECtJ)  -  CIREC 

"18600 - PT70HNT  T"NUE - 

38700  IEITSSW!?)  >1925*1 500* 1500 

11800  T9?5  H»TTFT4,4?5JKlR,VTWfJ5irflLT3T  - -  -  - 

38900  WRITE!  4, l 755) I. J 

19000  GO  TO  1500 

39100  C  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  CHANGE  SINGLE  ELEMENT  OF 


C  —  «70 


19700 - 7000 - WRtTgfKTTX,  70101 - 

39300  7010  FORMAT ! 1H  , ‘INPUT  COLUMN,  ROW,’  NEW  VALUE* ) 

39400  REAoinNY*ir»“JVTCHtr  - - - 

39500  C!J,D  •  XCMG 

T9600  rF!15SWT  7 1)7075 ,0400 ,1500  - - 

39700  7075  W» T TE ! 4, 4 75 1K8R ,1 SW, JS W*KL FS T 


39800 - WR  HT  14 , 7030)  T,~J,  XCHO - 

39900  7030  EO»MAT!lH  ,7! r?,7X) ,F17.6) 

40000  GO  TO  1500 

40100  C  ♦♦♦♦♦♦♦♦♦♦COMPUTE  ST0.DFV.1X1  IN 
40700  7100  TEttS5lf!7»l?tl0V7ir0,?l?0 


COL  1  FOR  X-8AR,X#X*BAW  IN  COL  4*3  ♦♦' 
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"50TOO  7110  WWlTFH*A?*>)K9rVISV*  JSIT.KUTTT 

40400  2170  00  2130  !REC  «  tVJHAY 

7  -  C 1 1 R  EC  •  4  !•  “ 

40600  C(IREC,1)  «  S'OR'TIN04<C(IREC,3>-Z*n/f  NO-i  II 

~~50Too  21  to  comtinot 

40800  GO  TO  1500 

40500  C  •*♦**♦*♦*♦*•♦♦  *i22--CJH.nJtATF  DH  FOR  EACH"  STITT W-«*OTTW  TNNXF' FORMAT' 
41000  7200  GO  TO  1400 

— 41  TOO - CT*************  wnfNTH  0 — SHORT  OngCUfTENTAriUN-i-AVKGS  »»»»»»»»»» - i 

41700  2300  0f>ENlUNIT-50,<NANE*»AVRGS.D0C#*:TYI>F-,0t!T,»RFA0Of«LYI 

*1100  00  7WN- 1,700  ~ 

41400  RE AD< 50, 2325, END-2355 H DOC  til  VI -1*9) 

CT500  WRTTE  FKTTY  ,  2  3  TOTTODCTT  TYT-TilT 

41600  2350  CONTINUE 

4ttoo - rm — P0ftH/m«>A8) - 

41800  7330  E0R8ATC1H  ,9A8) 

moo  7355  CtflSF  TUN  IT-50)  "  ~ 

47000  GO  TO  1500 

47T00 -  ~  f  NO -  - 


1 


J 


too 

200 

—jtnr~ 
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TTCT  77  1977  ♦  ♦*** 


C  TA5LT  SUBPROG  PFPLT ***** 
SUBROUTINE  TABLE 


.  1 
> 


t 


t  . 


2900 
1000 
II 00 
1200 
T100 
3400 
3500 
3600 


5700 


3B00 

4J00 

4700 

4100 

4400 

4100 

4600 

"4T100 

4900 

4000 

1100 


9200 


1300 
1400 
1500 
1600 
1700 
5  BOO' 
1900 
6000 
6100 
6700 


400 

C 

w 

- 100 

r 

TO  COWRUTF  AND  PLOT  PTJTFV  V A RTA5LTSV . 

600 

c 

JUNE  2  7  1976  N  FOFWOFF 

l 

- rmr 

C 

VAT  VFPSTON — NOT  1910  N.' BRAT  “  '  . 

BOO 

c 

900 

TDCLUOF  ' ' CONPFFCT.F OR* 

910 

DIMENSION  V AR 1 100*9 ) 

'  1000 

c 

'  .  "  . . 

1100 

C 

PROO 

1700 

TF I ISW-1 1600, 11,10  - 

1300 

10  NOSTN  -  0 

- MOO — 

NO*  1 

1100 

NX1-1 

~w" 

1600 

NX7-64 

1700 

NY1-2 

-  T600~ 

NY7-12 

1900 

AMI. 

- 2000 — 

87*7 

2100 

B  M I . 

7700 

B?*-« 001 

2600 

X20IM-3. 

2  TOO 

Y20IM-10 on. 

2720 

XMIN--UO. 

- 7710 — 

"TWAT»=40 “  ’•  -  ' 

7  760 

Y*IN--109. 

~  7710 

Y«AX--10 

2600 

11  WR  ITEMT  TX,  1101)  N0,PMIN,PMAX»X20M,Y?OIM 

1101 


FfTRMTdH  ,*60.'  STATnjNSTTNIJiTWTN,  Plir.T^DI  W,T?DT*"T 


%T  4,477.0 


RF4DIK!N,*)NOyP«!NiP«AX,X2D!6.r20Ii 
,T*TN«T - 


WRITE(KTTX,7000) 

9F8DCKTN,*) JWTO 
WRITEIKTTX,  7010)PLABL 
if  (noyes  <  kin,  xrmn  .to.iithfw 
RFAOIXIN*  7O20)PLABL 


FNOTF 

WRTTF(XTTX,704O)X1TN,XW*X;TNTN7Y*1X 

REA0<KIN,*)X1TN,XBAX,YNIN,Y*AX 

call  aGSFTriifTx/wrw.  txrpn - 

CALL  AGSETF ( 6HX7NAX. , XMAX ) 


CALL  AGSETFf 6HY7BAX. ,Y*AX) 
WR ITFfKTTX,  7OT0TXLABL 
IF I NOYES (KrN.KTTXI.EQ.l) THEN 
RE  AO I TIN ,  70201  XL'ABU 
CALL  STRIP(XLA«L> 

FNOTF 


WRITE  I XTTX, 70601 YLABL 

TF  I NOYFS IRTH,KTTX)  *TO  .1 1THF1 -  -  - 

RFA0(XIN*702O>YLABL 
CALL  STRIPIYLTBL3 
ENOIF 

~  -CALL  AGSE  I  FIIHFEAWE. ,2. > - 

RF  TURN 

*****  PLOT  RH-RTTrO  F0RNAT5  **»4n4»4*4R4»*44F4*4*4R****4* 

FORNATIIH  **INPUT  INOFX  OF  FIRST  POINTS*) 

FOONATTtH  ,  *  CWAWCr  RLOT  L ABTLT  OLOLXBET  TS : *, /  , 40 


I  Ait 


7000 

7010 


□ 


,T014) 
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6300  7020  FORMAT (10441 

6500  7040  F  ORMA  T ( 1 H  ♦•CURRENT 


*600 

FORM  AT  ( lH 

6700 

7060 

FORMAT! IH 

~  7100 

600 

conttnuf 

7200 

C 

7700 

18 

DO  180  J-4, 

7400 

DO  180  1-1, 

7500 

7800 

!  80 

TTTTTTT  “  0. 
DO  101  K»l, 

T^mr 

* ’CHANGE 


6 

100 


VALUES 

x-asts 

Y-AXIS 


OF  XMIN,XMAX»'YMIN.,YNAX*  **/*4F10.3) 
L'ABFLT 
LA8FL? 


TO" 

OLD 


CAttU 

LABEL 


7TT»T7v4Tr 


rypn 

10A4I 


IS***/ *  AH 


NO 

~79D7}  RFA0I12, 1010,FND-106JM,C^AWmnTMl 

8000  1010  FORMAT! I  4  *  A 1 2 ,  F5 . ? ) 

8100  GNAMF TK  )  1 9 : 1 71  »"*  «~A  VG  * 

8200  101  CONTINUE 

~  *75U - RT6 - T  F7W70T.  TT-T1W-*T^I - 

8300  RFHINO  12 

~!TA00  00  DO  175  JST«I,N0 

8500  IE0F-0 

TTCOO  OPFNIUNTT-KTP.TUMF-t 

8700  ♦  'UNFORMATTED*  *«ERR*1 75 ) 

■  '8  BOO - 


’OLD^VFORTf- 


-75 - artt  "trer*  r  1  mom 

IF(  IE OF. EO. -1100  TO  111 

X  -  Al*VRBL(NyiWA?*VRBLnnr7T>HT*VR0iLTTfim4A4*rriirEr,n 

Y  -  Bl*VRBL!NYll4R?*VRBL<NY?|4B3*VRBL(NY3)*B4*C«fR1FC*?l 
7  -  CT*7RBL  1571 J>C7  A7RBLT 
XPR-A1*VRBL(NX1I 
_ —  YRR-BIWRULTimr 
WT  -  1.0 
VARfTRFC  *  1  )*PF 
VARIIRFC*2)-XPL 
VATr(TRFC»31  "YTT 
00  950  M-1,6 


8900 
9000 

9100 
_9?00 
9300 
~95U0 
9500 

<7509 

9518 
95  Z  7 
9536 

’  9545 - YARirRFC*M»3  I  »'VRBL  TNV f  Mf | - 

9554  950  CONTINUE 

9578  IFITYSWriOl .FQ.I1TWFN  - 

9581  WRITE  (KOUT,  1621 H  VAR  !  I  PEC  «  K  )  *  Jt-1*  9 ) 

9587  FN0IF  - 

9590  1421  F09MAT19F8.3I 


9595 

9600 

9700 

9800 

9900 

10000 

142? 

940 

980 

985 

FORMAT!/) 

IF  !  TSSWI15)  >980 *'985*985 

WT  -  WGT  - - 

C! IREC»4)-X 

C  f  TREC*  5  I  «Y 

C  1  TREC  »6 ) *7 

— 

10 100 

10200 

110 

WR ITE !<L IST,*)PY,X,Yf7 

TD300 

1100 

FnPMAT!?I4,A?TTr,T7.0*3F10^7'  “ ' 

“  '  - - -  -  — 

10400 

113 

GO  TO  95 

0500 

10600 

111 

CONTINUE  - 

IF!  ISSWI10I.E0.-1.AND. JST.EO.NOI 

10700 - 

11300 

IT*  00 

- AOR  ITF(KOUT,  1427) - 

IF!  ISSW!6).E0.-UG0  TO  175 

171  I  JfT  -  0 

11500 

00  172  K-JMIN,IRFC 

11600 

I J9-T JW41 

_  - - -  _ 

- - - 

. .  -  - 

11700 

XOAT ( I JM  *  1 )  *  C ! K  *4 ) 

11900  17?  CONTINUE 

12000  TF!  J5T.0T.il  TWF19 

12100  CALL  AGSE  TF 1 11HBACR6R0UN0. *4# I 

17700  FNOTF  - 


120 

T7300  CALL  ERSTPTOflJATdtll  *  YDAT (1*1)) 

12400  CALL  ANOTAT»XLABLtVLABLtO*OtO,0.» 

12*500 - C*rr-E7"XT  <  XPAT  ,TnA7t  JDIW«ltUW»E 

12*550  IF  I  IS$VM5).EQ«— 1 I  GO  TO  175 

~T7500  YPR2“YPR  ♦F?*Y70TT5  - 


12700  XPR2»XPR-A?*X?0I* 

T7mu  XPR3-XPRVA2YX?T)T* 

12900  CALL  LINENCARl  XPR,YPR*XPR*YPR2» 

rJW0  CALL  LTNFNCAR FXP1T7,Y PT? , XPR TTYTHTT" 

13100  175  CONT INUF 

ITI50  ■"  lFCT55W(6ir.T0.-nCTT  TO 

13200  CALL  FRAME 

IT300  CALL  AGSETFf  EiHYTinTJi  * ITFJfe) 

13400  CALL  AGSETF ( 5HX/NAX*  »1.E36) 

- TTYTTO - CALL  A  CTETF  ( 6HT  /7TT1T.  , 1 TFJH1 - 

13600  CALL  AGSETF ( 6HY/MAX, *l«E36 1 

TTTOO  CALL  AGSETFf  IlWRICTOinjUND . , l *  ) 

13900  178  IFUSSWflOl.EO.-l I CUJSF  (UNIT-ROUT) 

13900  "  RETURN 

14000  FNO 


C  DATA  SuPR  PEPlT  *****  SEPT  15  1977  ♦♦*** 
SU9R0UTINE  D  A  T  A  <  MS  W  *  I  F  OF  l 


[  T  IITflTIIITIIITII  IIIIIIIHI3IIITII11I11I' 


C  PROGRAM  TO  RFAD  AMO1  SELECT  POTEN  DAT#* 
C  JUNF  2  7  1976  N  F OF OMOFF 
C  VAT  VFRSION--NOV  TWO.  M.BRAY 
INCLUDE  'COUPE PIT. FOR* 


900 

1000 

1100 

T70U 

1300 


MW  -  1 

rF  CMSW11«?D*?0T) 
CONTINUE 


2700 

7100 

2900 

3000 

3400 


3600 

3700 

3900 

3900 

4000 


4100 

4200 

4300 

4400 

4600 

4700 


IFLAG-0 

RETURN 

20  CONTINUE 

172  WRITEIKTTX,lf3)OAYl,DAY2 


REA0<KTN,*)0AYl,0AY2 

174  MR  T  TE IK TTY, 175 1  VENN ,YE*X,YNHN,YNHX 

175  F0RHATI1H  * 7HE-N  LTM*4F7.2I _ 

RFAD(KTN.*)YFMVfYFMY*YNWNF5rWIY 

WRITFIKTTX,l?7IZLT0»ZLG0* JDO VPN1N* PM1X 


177  F0RWATI9H  HRIOTM:  »  ?  f  X  ,  F  B  .  3  )  ,  X  ,'4H  J  DO  f  ,¥ 
RFAO<KIN,*)ZLTO,ZLGO*  JDO,PNINf<P'*AX 
RETURN 


200  TEITFLAC.EO.nGU  TO  217 
RFAOIKTp  * END ■ 2  90 I KHDG 


4975 

5000 

5100 

5200 

5300 


■5^00 

5500 

5600 

5650 

5700 

5900 


7011) 

5900 

6000 

6100 

6700 


251  TF ( VOAY-OAYl 1 790,252 1  252 
752  IFIDAY2-X0AY»290,254,254 

264  IF (XL0NG-XE9N)290»256*256 
256  IF  I XFPX-XLONG I260f258»259 


260  FF  f  XNMX-XL  AT  I 79  0*26 7» 267 

762  YPL  -  -111.1 7*(XLnNG-7LG07*C05rnXLlT>7VrOT/ll4.5«r7> 
♦  ♦  JD04FL0ATI  It:ON» 

YPL  •  Hl.l2*fYU*rr-2LT01  -  - 

DAY  «  XD AY 


.  l|l_VJ.|l|  I  1,(1] 


2615 

7620 

2625 

763X1 


TFIISSMI l 3)12620, 263*263 
TEITREC  -1 1 2630 *267532630 
MR  ITE (KOUTIKHOG 
WRTTEIROUTIXOO^ 


I 


PFPLT  5UBR  *♦♦*♦  PFPLT5UB  SFPT  15  1977* 
VAX  VERSION  NOV  1980.  N.BRAV. 


IL'UtAilL'HJlia'AAi.mij 


TO  EVALUATE  POLYNONI AL  OF  ORDER  N-l  tfITH  COEFF  CPftl, 
JAN  28  1976  N.  FOFONOFF 
OINENSION  CPU) 


-i<l 


800 

850 

900 

950 

1000 


l  v7U 

1100 

1150 

1200 

1250 

1300 


1600 

~T650 

1500 

7550 

1600 


1D7U 

1700 

1750 

1800 

1550 


n 


2300 

7150 

2600 

“7550 

2500 


2600 

~7650 

2700 

7750 

2800 


2900 

7950 

3000 

1050 


V  «  VO 

TFl V— VWN1I *  2*2 
IV-  V5N 


'  M  1'M'AI 


3  V  -  VNX 
6  POLY  »  0.0 
X  -  V  -  0V8 

no  m  r*T»x 

NT  »  N  -  I  ♦  1 


RETURN 

FNO  -  -  ' 

DPDV  FCN  A*********************************** 

FUNCTION  D P TTY  f VOYUYWVCP, N ,VHNf VHY1 - - 

*************** *********************************** 


TO  COMPUTE  DERIVATIVE  OF  POLYNOMIAL 
JAN  28  1976  N.  FOEONOFF 
DIMENSION  CPfl) 


V  «  VO 

TF(V-V5N?I,2,2  - 

1  V  -  VHN 

7  TFIV5T-V)3,6,6  - - - 


1900 

3 

V  -  VNX 

4 

WWT  *  N  -  1 

2000 

X  -  V  -  OVN 

7050 

OPDV  -  0.0 

2100 

00  20  I  -1*N81 

7T50 

N5I  -  N  -  T 

2200 

20 

OPOV  -  DPDV*X 

PND 

8 NO  FCN  t********************** 999*99999 
FUNCTION  BN0Uv781N,f8AX) 

*******  9*  *******  ****  ***  ***  ******  ***********  *  * 


5N0  -  Z  ‘  " 

TFI7-75IN) 10*20t'20 
lO  8 NO  «  7NT7T 
RETURN 


30  5N0  -  ZMAX 
60r  9FTURN 
ENO 

JFF  FCN  *****  FT58T  ******* 


I 


124 


1 


3TU0 

3150 

y?oo 

3250 

TTOO 

3350 

~3400 

3450 


FUNCTION  JPR( !IFC»NP**M51 


C  GENERATES  PRESSURES  CORRESPONDING  TO  IREC. 

C  - 

C  OCT  28  1975  N  FOFONOFF 

OINFNSI ON  NPR ( 1 ) 


- 3500 

3550 
~  3000 

3650 
T70TT 
3750 

~c - 

no  ioo  j-i*ns 

IFfTREC.LT.TIPRTJTlTHFK 

JPR.NPR(NS*JIMIREC-NPR<2*N5* jm 

RFTURN 

FNOIF 

i  ' 

* 

•4 

-  i 

3800 

3850 

7900 

3950 

4WU 

100  CONTINUE- 

JPR-NPR f  2*NS1 M IREC-NPR  <3*NS  >  » 

RFTUPN 

END 

^*•4  ******************** t******************************** 

\-s 

4050 

SUBROUTINE  STRIPtAI 

'■) 

4150 
“  5700 

4250 
6700 
4350 

C 

C  STRIPS  TRAILING  BLANKS  AND  PUTS  A  t  AT  THE 
C  IN  ARRAY  A  FOR  CONEORNANCE  KITH  NEAR  PLOT 

C  N.RPAY  17NOV80 

C 

f  NO-  OF  CHARACTERS 
PACKAGE  LABELS 

r- 

4400 

4450 

tJfHf  NSION  ft  f  l"l 

C 

4500 

B-’  *  "  . . . 

-i 

: 

4550 

no  ioo  j-i»io 

4600 

K  *  l"0— J  ♦  1"  - 

4650 

If  <A(»0-B»200*'100*200 

9750 

5800 

4850 

5900 

4950 


5050 
5T0TT 
5150 
'  5700 
5250 
— TTOTT 
5350 
5400 
5450 
5500 
5550 


5650 

570U 

5750 

5500 

5850 


7VOU 

5950 

7>mro 

6050 

6100 


700  NCH«K*l 

TFINCM.GT.IOINCW-IO 
A (NCH !•* t  * 
RETURN 
FNO 


ItAHUL 


0?POV  FN--SECOND  DER  TVATTVE  OF  POLYNONIAL 
FUNCT I  ON  07PnVTV0,DVNfCP  .wvNN.'vion 


nr^rmni  rTUFBVDrr 

0THFN5 ION  CPfll 
V  -  VO 

T  F I V-V8N) 1*2*2 


i  »  “  inn 

7  TEt  VHX— V I  3*  4**4 

3  "V  ■  vhx 

4  NN1  -  N  -  2 
X  «  V  -  0T8 
07P0V  •  0.0 


* 


ril 


N*I  •  N  -  I 

1TNI2  -  NET  -  1  - - 

20  02POV  «  02P0V4X  ♦  FLOATINNIRNHIZMCPf HMMI 
return  -  - -  - 


6150  FNO 

6200  C  SEAMATEP  PROPERTIES 
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r.>xi 


6300 

6350 

6400 

6450 

6500 


0171) 

6600 

6650 

6700 

6750 

6800 


6500 

6550 

7000 

7050 

7100 


7200 

7750 

7300 

7350 

7400 


LVlil 


7500 

7550 

7600 

7650 

7700 


7600 

7850 


C  SCO  ♦♦♦•♦ 

FUNCTION  SGOIS) 

C 

C  STOWl-O  XNUOSEN  - 

C  FEB  15  1976  N.  FOFONOFF  _ 


SCO  -  C  (6.  76786136E-6PS-4.8249614F-4)PS'*0.81487657T»*S 
X  -0.0934458632  ~~ 

RETURN 

FNO 

C  SCT  FCN  ♦♦*♦*♦♦♦♦*  _ 


I  »  3  * 


lik'd 


irl 


C  SIOEl-T  KNUDS EN 

C  FEB  15  1976  N  FQFONOFE 

C 

SC  -  SGOIS) 


LVI.llIltI.il 


X  ♦4.53168426)*r>/<T*67.26IM  (  <  t  l‘i667E-8PT-8. 164E-7  F*T 
X  ♦1.0O3F-5  »4n#SC»(  f-1  .0843 E -6*T*9.  8T85E-5>  *T-4.786TE-3>47 
X  ♦l.OlPSG 

RFTURN  . 

FNO 


v.  ruUI »  iu  pi  ur  i  i  *  i  r  run  )r*n«icn  cuvou 

C  *♦*♦♦♦**♦*♦*♦♦**♦♦♦****♦♦♦♦♦#♦*•♦•**♦♦♦♦♦*♦♦♦*♦♦•♦♦♦♦** 

REAL  FUNCTION  E0S80I PI *T*S1 
C  *•**♦♦♦♦**♦♦**#♦♦*♦*♦♦**♦♦♦•♦****♦♦•♦♦♦♦♦♦♦♦ 

C  EOUATf ON  OF  STATE  FOR  SEAMATFR  PROFITJED  FEDOTS  TOOT 
C  RFEERFNCFS 


l»M 


fMg  rrxvDi.n 


C  JPfJTS  NINTH  REPORT  197Bt'TENTH  REPORT  1980 
C  UNITS* 


7900 

7950 

8000 

C 

c 

c 

PRFSSURF 

INPUT  PRESSURE 
TF8PERATURE 

P 

PT 

T 

BARS 

DE  C  TEARS 

DEG  CELSIUS  flPTS-681 

c 

5 alInTTV 

S 

NSU  I'lPSS— 78 y 

8100 

c 

DENSITY 

RHO 

KG/N943 

¥150 

c 

SPEC.  VOL. 

ETJS80 

WRrj/RC 

8200 

8750 

8300 


BJ1U 

8400 

5450 

8500 

8550 

8600 


.r.Vlij 


8700 

8750 

8800 

5550 

8900 


9000 

9050 

9100 

“9150 


C  CHECK  VALUE*  E0S89  -  9.435561E-4  NP*37KG  FOR  S  -  40  5SV, 
C  T  V  40  OFC  C*  P  -  TODD  B APT. 


C  NOOTFIFO  TO  TAKE  OB  TNPOT  PRESSURE  «•  AND  OUTPDT  IN  C5R*S7GW  28N0V80 

c  N.wny  -  -  ■  . "  - - 

REAL  P1»P,  T«S*RH0*SR*R1*R2«R3JLR4 
REAL  A,P*C,0,E*AT,Bl*AN,BNfR,KOi'RW 
C  FOUIV 


LlkH  17.1 44x15: 


mi  rmn*.rn.n.uwiiCT  * 


EQUIVALENCE  I A««'A,RL*R01 * (KW,KOV* 1 
C  rrTWFf  T  PRESSURE  TO  TART  IWO  SQUARE  ROOT’ SI 
P  •  P 1  ♦  •  1 

SR  •  SORTIABSmi  - 

C  CONFUTE  OENSITY  PURE  MATER  AT  ATM  PRESSURE 


I wa* ii i . k I sQ.i  a. b  *\iii 2 i.rui  j~i  nil 


X-9. 095290E  -3  )  ♦  T*Ui  793952E-2 1  *Tx999.  842594 

C  SEAWATER  OENSTTT  AT5  EWETT. - 

R 2  -  <«  I5.3875E-9*T-8.2467E-71PT4  7.6438E-^5>PT-4i0599F-3)RT 
X46.24493E-1  - 
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9700  R3  -  (-i.6556F-6*T*l .02?7F-6>*T-5.?2586E-3 

9250  R5  -  5.8315E-5 

9300 - ♦  rirsr  *  V7T*s-*~wt - 

9350  C  SPECIFIC  VOLUME  AT  ATNOSPHER IC  PRESSURE 
9500  ALPHA  -  I «  E  ♦  3 7WH0 

9550  EOS 50  -  ALPHA 

95tro  tecp.eo.o.otreturn 

9550  C  CPHPUTF  COMPRESSION  TERMS 

"9300 - F-.-rT.~I  A97E-1  OAT  V?  .08 1  ST-  8T  *T  -9  V  93  A  *  F^T 

9650  BW  -  I5.2787E-8*T-6.12293E-6)5T58.50935E-9 

9700  B  ■  PW  ♦  EPS 

9750  C 

■9W00  T»  «  I.91075E-5 

9850  C  -  <-1.6078E-6FT-1.0981E-5)FT*2.2838E**-3 


9950 

10000 

10050 

lOTOO 

10150 

-19200" 

10750 

10300 

10350 

10500 

10550 

"T0500- 

10550 

10600 

10650 

T0700 

10750 

"imroo 

10850 

19900 

10950 

11000 

11050 


X*3. 739908 

A  .  f05S*  ♦  C1*S  ♦  AW 

PI  -  1-5. 3009F-5*T,1.6583E-?)*T*7. 955**7 

*1  -  (  (-6.1670E-55T«-1.09987F-?»-*T-0.603559)5T*56.6766 

nr  «  r  r  t  -  5 .  ip  *?mrr=5n-  rrrjou  atte-^? i  *' r  - 1 . 3 77T usm - 

X*  l  58.570* ) *T* 19657. 71 
KO  «  ( R1*SR  ♦  A 1  1  * S  ♦  KW 


K  .  fP*P  ♦  A  I  •  P  KO 
ALPHA  -  AIPHAM1.0  -  P/KI 

- ptiiPo  *tm  - 

EETU»N 

FNO 

V350P  FCN  OCT  7  1980 

PEAL  FUNCTION  V350P!FI) 

~ vp FC'TnT~~Ynuijn:~ re r * 9771*9 1  for  s  -  35  nsu  rmr=7ir? 

TF NPrR A TUP*  0  0*0  CFLSIUS  IIPTS-68)  AND  PRESSURE  IN  DECIBARS. 
FOUATION  DERIVED  FRO"  ETTS8  0 

CHECK  VALUE:  V350P  -  9.337531E-5  M5R3/KG  FOR  7  -  1000  BARS. 
"ODTETED  TO  ACCFPT  TNFT/T  PRFTSURT  TN  OB  AND  OUTPUT  5P. VOL  T9 
CM553/G"  78  NOV  80.  N  PRAY. 


11150 

11200 

11250 

11100 

11350 

11500 

11550 

11900 

11550 

11600 

11650 


ALPHA  -  9.  7?662’E-5R<1.0-P/I21582.27»I3.35951»5.032EJ-55P1*P1I 

ALPHA  •  1 .  E  ♦3RTSVPHA  - 

V350P  -  ALPHA 

RFTUPN 

END 

rw  fc'N  'TRFT*— orr  ~7~TW0  TMT» - 

REAL  FUNCTION  OEPTHIP l ,L  AT) 


OF  PTH  IN  PETFRS  FRO"  PRESSURE  IN  DECIBARS  OS TNG 
SAWDERS  AND  F0E0MUFF*5  "ETHDU. 

DEEP-SEA  RES.,  1976 ,f?3 , 109-1 1 1 . 


11750 
11800 
11850 
11900 
11950 
12  ODD 
12050 
12100 
12150 
I?200 


REAL  LAT 
P  -  Rl*.l 

V  -  SINILAT/57. 29578) 


GR  •  9. 780318511. 0 ♦  C 5.  2708€-3*2. 36E-5PJT ) **)  ♦  1.097E-55P 
DEPTH  -  f  f (-1 . 87F-1 l*P ♦  779E- 7 ) *R- 2 . 297 7F-1 1 **♦ 97 .76591 *P 

DEPTH  •  0EPTH76R 
RETURN 


I 
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12250 

12100 

T7T5TT 

12400 

12450 

12500 

12550 

12600 

T755TT 

12700 

12750 

12600 

12650 

12000 

TTOTT 

12954 
12958 
12962 
1 7966 
12970 


PND 

ATG  PCN  ♦♦*♦*♦♦*** 

- rUNCTTW  A  T OTPVTTTi - 

************************* 

AO  TAP  AT  I C  TEMPERATURE  G» AO ! ENT.  JMtrOPN  1973. 

OS  -  5  -  35.0 

- ITG  «  -«Tr=2;T687F-16*r4i'.B6  761-141*1=47670*8-133*? 

XM(?.7759E  — 1?*T— 1.1151E— 10l*DS4|f— 5.4481E  — 14*T 
X*8.733F-l2>  ♦T-6.'T796P-10  >*T*l."87r4lF^81T*P’ 
X*(-4.2393E-8*T>l.8932E-6»*DS 

7*  ( (  6 . 6228 F—l 0*7—6. 836E— 8  I  *f  ♦  8.  5258E-6)  *T*  3.  58038-5 
®ETURN 

- PlfO - 

OVA  PCN  ♦  *♦♦♦  PTSfll  **♦*♦ 

FUNCTION  OVA(P«TyS) 

4* ********* ** *****  ***** ** 

SPFCtFIC  VOLUME  AN0MAL7 


12978 

12982 

12986 

17V90 

11000 


11100 
rn  so 
11200 
1T750 
11100 


OVA  -  SVAN(P,Ty5»SPV| 

RFTURN 

FNO 


SVAN  FCN 


(IL'I* 


IIUr»  iTBUir, 

************************* 

SPFCTFrr  VOLUME  AN0WIUVnF5 
FEB  T5  1976  N  FOFOflOFF 
V  -  F0S80 ( P »T *3 1 
SVAN  -  1.0E5*<V  -  V350PIPM 


t 


4 

i 


11500 

13550 

13600 

13650 

11700 


l  jnu 

11800 

13850 

11900 

TT95U 

14000 


FNO 

THFTA  FCN  *************** 

FUNCTION  THETAtP0,T0,S,PF» 

******************  1  *f*r*  *****  *  *** - 


ir  L OC « l  ruirni  i  rl  icnrc««ium  ■  i 

FOURTH-ORDER  RUNGE-KUTTA  integration  using  STEPS  OF  100  08 
PR  LFSS.  IRALSTON-WILF  VOL  1  PTT5VWT61 


OCT  12  1975  N.  F0FTW7FF 


■\  14050 - 

P  14100 

i;  14150 

|j|  14200 

I  14250 

f  14300 

- p-Tr-ptr 

T  «  TO 

H  -  PF  -  P 

N  -  ABS ( H ) 71000. 0  ♦  1.0 

H  -  H/FLDAT IN) 

00  10  I  *  1  *  N 

-  -  -  -  • 

fK  *  H^#T6  I  F<S) 

14400 

T  •  T  ♦  0.5*XK 

14450 

0  **  XK 

14500 

P  -  P  ♦  0.54H 

14550 

XX  •  hratgip»t„ti 

.  • 

14600 

T  -  T  ♦  0.29289122MXK-Q) 

.^•A.I.V  111  1  < 


LUtfLLU1! 


14700 

14750 

14800 

14850 


XK  *  H*ATGtR*T,S) 

T  -  T  ♦  1 . 707106 78T •  f VIC-01  - 

0  -  1.41421356 ?*XK  -  4 • l 2 13203 44PQ 
P  ■  P  ♦  0.5*M 


128 


woo 

14950 
1500U 
15050 
T5100 
15150 
157U0 
15750 
15700 
15350 
15 400 
15450 
15500 
15550 
15600 
15650 
15700 

15750 
15500 
15650 
15900 
15950 
16000 
16050 
16100 
16150 
15700 
16750 
16300 
16350 
16400 
16450 
16500 
16550 
16600 
16650 
16700 
16750 
"1550  O’ 
16850 
16900 
16950 
17000 
17050 
17100 
17150 
17700 
17750 
17100 
17350 
1741)0 
17450 
17500 
17550 
17600 
17650 
17  700 
17750 
17800 
17850 
17900 


XK  -  H*ATG(P.T*SI 

T  -  T  ♦  <XK-?.ff*Q)/6.0  _ 

THETA  -  T 

RETURN 

PNO 

T65  FCN  ****♦♦*•♦♦ 

RJNCTIQN  T68t  T  )• _ 

TO  CONVERT  T-48  TO  T— 68  TEMPERATURE  SCALE 
FEBf  15  1976  N  FOFONOPF 

T68  -  T  -  4.4E-6*T*< 100.0-T) 

RETURN 


T48  FCN  ♦♦♦♦*♦**♦* 

FUNCTION  T46(T) 

************************* 

TO  CONVERT  T-68  to  T-48  TEHPER1TTWF  SCALF 
EE8  T5  1976  N  EOEONOFF 

T48  -  T  ♦  4.4F-«-6*T*f  100.0-T) 

RETURN 

FNO 

OVTTT  FCN  ♦*♦*♦♦♦•♦* 

FUNCTION  OVDT(P*T»St 

- - 

OF R IV AT l VE  OE  SPECTFTC  VOL.  WITH  TEMPER ATURERIE5 
EFBT  70  1976  N  FOFONOFT' 

N  *  0.75 

OVOT  -  (5.0E4/H14|EOS8OfP,T*H,S»-E0S8f)IP,T-H,SI  » 
- ^turn - 

FNO 

OVOd  FCN  **********  - 

FUNCTION  OVOS  C  p»'T»S  I 

OFRTVAT I VF  OF  SPECIFIC  VOL.  WITH  $ ALTNI TY*IF5 

FT  5  70  19  iT6  ~N~TnFONTTFV - 

H  -  0.5 

OVDS  -  f5.0F4/H)*fEOS80|P,T,S*H»-EOS80f  P,’T,S-H3  I 

RETURN 

FNO 

WfjF  TCN  *****  *  **** 

FUNCTION  OVOP{P«T,S) 

****************  ********** 

AO  I AM AT  I C  DERIVATIVE  OF  SPEC.  VOL.  WITH  PRESSURE*IE5 
FEW  70  1976  N  FOFONOFT  - 


OVOP  -<5.0E4/H>.ME05  80<P»H*T,S>-E!)S80(P-W*T*SH 

X  ♦  ATGfP,TfS>*WroTfPtT,5>  - 

®E  TURN 

FNO  - 

OKOT  FCN  •*♦♦♦*•*♦• 

AOTARATTC  C0MPRFSSIB1LITV  TEMP  DTRTVITIVE 
EE 8  70  1976  N  FOEOMOFF 


129 


1 7950  M  -  1.0 

1*000  0K0T  -  I0.5/H»P!0VDPfP,T*H,S)  -  DVDPf P,T-N,S») 

TffOW  '  RETURN . 

1*100  PNO 

18150  C  0*05  FCN  •*♦♦♦♦*♦♦♦ 
l*?00  FUNCTION  OKDSIf>,T,SI 

1*250  C 

1*100  C  *DfA*ATIC  CO"P«gSSl*n.ITV  SALINITY  DERIVATIVE, 
18*00  C 

18*50  H  -  2.0 

1*500  OKOS  -  (0.5/H»*(0VDP<P  *T,S4.H)  -  DVDP  I  P 

18550  RETURN 


1*600  END 

Twspro — r  t#  rw  * ***  ***  *?* 


1  *  TOO  FUNCTION  S  AL  !  P  ,<T  ,0  I 

1875 0  C  *****♦*♦*♦**♦****♦******♦* 

18*00  C  COMPUTE  SALINITY  GIVEN  PRE  SSURE  , TE  MPFRATUR  E  ANO  SPECIFIC 
15850  C  VOLUME  ANOMALY! 10**5*DELTA1 
18<»00  C  FFB  16  1976  N  FOFONOFF 


- r*950 

C“ 

19000 

K  -  0 

19050 

SAL  •  15.0 

19100 

10  S  «  SAL 

I9I50 

SAL  •  S  ♦  (D-Sf 4N(P,T,S,V>)/DV0STF,TVS1 

19200 

R  -  K*1 

1 9750 

TF  ( K-50)  20,10,10 

19300 

20  TF( ABSISAL-S) -0.0005110,10, 10 

19150 

10  RETURN 

— ■ 

19*00 

FNO 

19*50 

C 

NOYFS  FUNCTION  *♦**  DEC  1  1979  ♦  ****♦“ 

rj 

19500 

FUNCTION  NOYES! KIN, KTTXI 

'  J 

19550 

c 

“■*♦»**•**♦*  ®p*  *****••*•• 

•1 

19600 

c 

RF  TURNS  l  FOP  YES  -1  FQ®  NO 

19650 

NOYES  -  0 

19700 

1  RPA0fKIN,10ILB 

19750 

lO  FORMAT ( A? I  - 

19800 

c 

TF"!  LB.FO.  ?HYF  iNdYFS^l 

19900 

rF!L*.F0.2HN0IN0YE5— 1 

19950 

TF ( NOYFS ) 10,20,10 

-  - 

20000 

c 

F  R  RflP 

20050 

20  WRITF!KTTX,IOOI  -  - 

.  . 

20100 

100  FORMAT ! *  t  YES  0*  NO?  *) 

20150 

00  TO  1 

20200 

c 

20250 

10  PETUON  - -  ' 

‘ 

20100 

ENO 

20150 

c 

CTflSO  F  ILF  **♦♦♦♦♦  JULY  15  1977  «*******♦*♦*♦*♦*♦* 

20*00 

c 

YHSAL  FCN  ♦**♦♦*♦*  JULY  6  1977  ******♦♦♦*♦**♦♦♦*♦♦ 

4 

1 

TOT5T5 - FUNCTION  THSAH  KTNiYI - 

70500  C  ♦••♦♦♦♦♦*♦**♦******♦*♦ 

20550  C 

20600  C  T«KP*  UP  TO  25  CUBIC  SPLINES  TO  GENERATE  A  SALINITY  FRO" 

20650  C  ®OTFNTTAL  TEMPER ATURF  REFERRED  TO  THP  SURFICTY. INPUT  UATA 
20700  C  CONSISTS  OF  LOWER  SPLTNF  BOUNDARY  FOLLOWED  BY  FOUR  COEFFICIENTS, 

20750 - C  COEFFICIENTS  ARF  FRO"  TWF  ~FTT  UF~IRBT~  AND’  BRAT  119811  TO - 

20800  C  ISELIN  ANO  WORTHINGTON  METCALF  THETA-SAL  DATA. 

20850  C  - 


20900  OIMFNSION  CIS, 25) 

20950  C  OAT* 
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2rooo 

21050 


miy 


21150 

2T200 

21250 

7T300 

21350 


MCliltl 


21450 

71500 

21550 

21500 

21650 


21750 

21000 

21850 

7T900 

21950 


e.  cmw 

22050 

22100 

22150 

27200 

22250 


£.  C  JUU 

22350 

27400 

22450 

72500 

22550 


DATA  C/0. 00, 34. 738063, 0.0, 0.0, 0.0* 

♦  0. 50, 34. 7  3  805  3  v.'t  07290,.  5848497-02, -.253429E-02, 


ritnlniiu  i  ■  w  wta  i#i,i  %  ifki.rj  u  i  gjj  \ 


♦  1.50, 34. 85029  7*. 127  785.. 529320E-01, -.  135  3  79* 

*T.  75, 34. «83436,Tl  28868  ,-i 48582  87-01 ,  -.  WW, 

♦2. 00, 34. 91 05 87*. 802 174E-01,-. 146093*. 228920, 

♦2.  25,34.925087  r.'50091FF-01  T.25T484E  -0 1  ,-;7fc73B2F-0  1  ,r 
♦2.50,34.938790 578 544E— 01,  .  552526E  — 02  »,— •  359945E— 01 , 


1 1  airciM  cuttaivflKiiviiijSiii 


♦3.00,34. 9645  75, ..‘360969  E-01*  —  .  495364E— 01 * • 509274E— 01 , 
*T.  20, 34.9702  20, .2  73 936F-0I,-. 189292 F-0 1  »V5806 837-01 , 
♦3.40,34.9744  06,.217901E— 01,.157868E— 01  ,«•  4T9  730E— 02  * 
*1.60,34.979434  2  86805  E-OI,.  1859757-01,-.  794T72F-01, 

♦3.80,34.985679  /• 3 25 895 E-01, •  102958E-02,— .279688E— 01 ,' 


111  Wf.T'VKW’f  A  *T  *i  4  ■ 


♦5.00,  35.01238,  .'175223E-01*.  357759E— 02 , •  114377E— 02, 
*7.00.35.070  89,.-4*v5579F-01 ,'.T04186F-01  ,.3655927-05, 
♦  ro.  00, 35. 301 74,. -10«423,.  1051 727^-01, -.763343E-03, 
*rr.  on  ,35 . 7  01 05 ,7.T50916  ,“i  36 4T90'F=07,  .31 08057 -04 
♦16,00,36. 18748, .1 73 643,. 3929267— 02, —. 689732 E— 024 


,3n.-n  f^3,t  tow  r  nt-ui,-« 
♦21.00,36. 90401 18V0 .0,0.0, 0.0415*0.0/ 


HAT  A  KNOTS/22/ 
250  X  -  0.0 


U'J  7  1  .  1 

nr  -  cn.ll  -  T 
TFn)T>305, 320,370 
305  x  «  -or 

310  C0NTTNUF 
320  0  -  X 


♦ 


4 


22650 

22700 

22750 

22800 

22350 


C  c  TVU 

22950 
23000 
23050 
23  TOO 
23150 


C  KUU 

23250 

23300 

23350 

73400 

23450 


c  J  /uv 

23550 

21600 

23650 

21700 

23750 


TF(  10)325,325,330 
375  TO  -  1 
0  -  0.0 

310  TH5AL  •<  f  C  <5  ,TU1*TT*C  <  5 
PETUPN 


[  0 1 1  W*CT7VTD1” 


V88U  FUNCTTONS  PFPLT  ************* F************************* 

FUNCTION  VRBLfHVRI 


[Jill 
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25300  F3  -  -.050968*F24F’2/F0 

- 26350 - F5~«  OPDVl 5F iUVNVCP  ,N, FT,F2 T 
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132 


2  7100 
27150 


F5  *  -50.96fl*GK»P0F**2 
585  NVR*  -  NV»  -  57 


irii  mmi  m  uiunAjfiri  n.>j 


[.l.fl.VFlill 


'! 

ME 


29 950  VRP!  -  71*Z1 

30000  770  Y  -  VHP I  ♦  ?*PBF*r&F 

10050  VRBL  -  .5#VPPI*P5Yf6*Y 

30100  PFTUPM 


30150 

78  VRP  I  -  3 • *  7 1 *71 /NOP 

30700 

GO  TO  770 

30250 

7<>  VP P I  -  16. 

30300 

GO  TO  770 

30350 

80  VRP  I  -  71*7 1  /  NOP 

30400 

800  Y  -  . 4* VP P 1*1 VRPT  ♦  2* PDF ) 

30450 

VRBL  •  GK *  *  7 * Y 

30400 

PFTU»N 

- 3TT55TJ - 

8 1  wr  ■  T6. 

30600 

GO  TO  800 

10650 

87  VR8L  -  -1 

30700 

RETURN 

30740 

83  VRBL  -  -l 

30800 

RETURN 

30rt*>b 

84^ W5IT  -  ~D VTTP  (  P  F  ,  THE,  W1 

30900 

RETURN 

30950 

85  VRBL  -  OVOS(PF*TMP*SHFI 

31000 

RETURN 

31050 

86  VRBL  -  OVOP(PF,THF,SHF)*POF 

31100 

RETURN 

jii5u - nr vrbl  ■  ~=T5WTTFyypr^ry 

31200  RETURN 

31750  88  VRBL  -  -1 

31300  RETURN 

31T50  89  VRBL  -  -1 

31400  RETURN 

“TTT'TO - 9  0' VRBL  «  -1  ~ 

31400  RETURN 

3V550  01  VRBL  -  (5F-F3)**2 

31600  RFTURN 

31650  02  7  -  -THP*PDF 

31700  VP BL  -  1/7 


31800 

93 

VRBL  - 

l/( E0*F01 

TT550 

BFTURN 

31900 

94 

VR8L  - 

THI ♦ THI 

11940 

RFTUBN 

37000 

95 

TH8  5  - 

CfTREC«4t 

lijam 

32100 

RETURN 

371 5TJ 

96 

TH8"  « 

C(IBFC»41 

32200 

VRBL  - 

-( THF-TNN5)/(THP*P0F 1 

32T50 

RETURN 

32300 

97 

VRBL  • 

THI 

37350 - 

R*.  turn 

37400 

98 

FO  -  1. 

O/EO 

37450 

VRBL  » 

-  .  5  *  F  6  *D  2P0VE  OVF  t  H  VW,CF  *.WWF3  *T  21  *  (  P  D  F  *E01  **7 

37500 

RETURN 

34400 

200T7 

FNO 

34450  C  KOAY  ECN  ♦  ♦*♦♦  K0AV5  JULY  6  1977  **♦*♦♦ 

14  WO - ru qrf  TON  rt  0  AT  ( TO  JfflTTWl - 
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34650  C  USF5  LAST  4  OICTTS  OF  JULIAN  DAY.  A00  244000(7  TO  GET 
34700  C  FULL  JULTAN  OAY. 
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35100  20  *  -  TWO  *9 
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T  -4. 240F-1  3 T*TG *12. 13r=T5*T-l,24064T-T21»Ti-6 ^616=143*7 
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